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SUMMARY 


An  attitude  control  system  to  regulate  the  temperature  of  a 
manned  lifting  spacecraft  during  reentry  into  the  Earth's  atmosphere  is  proposed. 
Its  use  prevents  the  peak  skin  temperature  that  is  experienced  during  the 
reentry  from  rising  moderately  beyond  that  which  would  occur  during  an  equili¬ 
brium  glide  of  the  same  vehicle. 

The  effects  of  Earth  rotation  and  oblateness  upon  the  performance 
of  the  attitude  control  system  were  found  to  be  moderate  and  predictable.  The 
maximum  temperature  increment  associated  with  them  was  found  to  be  only  100  F° 
for  the  worst  set  of  initial  conditions.  The  cross  range  shift  of  the  foot¬ 
print  due  to  rotation  was  found  to  be  within  70  miles  of  the  value  that  would 
occur  for  the  corresponding  orbit  in  vacuum.  Oblateness  could  generally  be 
accounted  for  by  using  the  effective  initial  glide  angle  relative  to  the  Earth's 
surface  rather  than  the  geocentric  initial  value  relative  to  the  central 
coordinate  system. 

The  results  of  density  variations  in  the  Earth's  atmosphere  were 
not  serious.  Large  increases  in  the  maximum  skin  temperature  occurred  only  when 
extremely  large  spatially  random  density  disturbances  were  encountered  by  the 
vehicle . 
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I.  INTRODUCTION 

1.1  Historical  Background 

When  the  problems  of  manned  space  flight  were  being  solved  in 
the  late  1950 's  and  the  early  1960's,  the  state  of  the  art  with  regard  to  pro¬ 
tecting  reentry  nose  cones  was  oriented  toward  an  ablative  type  of  material 
since  the  heat  pulse  for  ballistic  missiles  is  characterized  by  intense  heating 
rates  which  normally  last  for  quite  short  Intervals  of  time.  Thus,  the  occurr¬ 
ence  of  manned  space  flight  so  soon  after  the  successful  implementation  of  the 
ballistic  missile  program  produced  technological  solutions  based  primarily,  on 
concepts  developed  for  the  ICBM  system.  As  a  consequence,  the  heat  shield  con¬ 
cept  which  is  presently  associated  with  the  manned  capsule  configuration  seemed 
both  a  natural  and  logical  extension  of  the  techniques  which  were  and  still  are 
generally  accepted  practice. 

However,  the  capsule  shape  is  able  to  attain  only  a  very  low  ratio 
of  lift  to  drag  (L/D).  Therefore,  only  very  limited  manoeuvring  capability  is 
available  to  the  pilot.  He  is  unable  to  land  the  spacecraft  outside  a  relatively 
small  area  once  the  retro  rockets  have  been  fired  and  the  deboost  phase  has  been 
completed.  As  a  result,  only  a  small  portion  of  the  Earth's  surface  can  be 
reached  during  any  given  twenty- four  hour  period.  This  is  because  each 
successive  orbital  track  is  displaced  by  about  22.5  degrees  as  a  result  of  the 
Earth's  rotation.  Therefore,  during  that  time  period  there  would  be  only  about 
sixteen  passes  over  the  whole  Earth  each  yielding  only  a  thin  strip  of  coverage. 

As  a  result  of  the  above  restrictions,  a  large  pickup  and  re¬ 
covery  force  must  be  maintained  in  the  principal  landing  zone  during  the  entire 
course  of  an  exercise  in  case  of  an  emergency.  This  procedure  is  unsati sf actory 
for  many  reasons.  First,  and  perhaps  the  most  obvious,  is  the  sheer  cost.  How¬ 
ever,  and  this  reason  is  probably  of  greater  importance,  weather  conditions  may 
dictate  the  choice  of  landing  area  and  a  sudden  storm  might  require  the  movement 
of  the  pickup  force  to  a  more  suitable  position  within  an  impossibly  short  time 
interval.  Other  difficulties  are  also  apparent  upon  even  a  preliminary  analysis. 
The  efficient  solution  to  the  landing  problem  appears  to  require  a  different 
shape  of  reentry  vehicle.  The  basic  requirement  is  greater  manoeuvrability, 
which  requires  a  higher  L/D.  Thus,  the  capsule  configuration  must  be  abandoned 
in  favour  of  a  winged  vehicle  or  lifting  body. 

Another  obvious  requirement  for  space  flight  for  the  foreseeable 
future  will  be  the  minimization  of  space  vehicle  weight.  It  is  not  necessary 
to  point  out  in  any  detail  the  well-known  correlation  between  the  cost  of  the 
first  stage  booster  and  the  mass  of  the  final  payload.  Therefore,  the  first 
step  is  to  design  a  vehicle  with  a  moderate  L/D  which  possesses  no  excess  pro¬ 
tection  to  enable  it  to  reenter  the  Earth's  atmosphere  and  land  on  the  Earth's 
surface.  Considering  the  minimal  degree  of  protection  available,  the  basic  flight 
path  upon  which  the  preliminary  design  is  based  should  be  a  simple  trajectory 
for  which  the  problems  of  protecting  the  spacecraft  and  its  occupants  are  as 
moderate  as  possible.  However,  when  the  choice  for  the  basic  path  and  the  method 
of  protection  is  being  made,  consideration  must  be  given  to  incorporating  provision 
for  departures  from  the  basic;. path.  The  departures  are  necessary  to  allow  the 
vehicle  to  manoeuvre  to  landing  sites  other  than  the  basic  terminal  position. 

Before  considering  the  choice  of  the  vehicle  and  the  basic  path, 

5wo  general  areas  should  be  reviewed.  The  concepts  involved  in  protecting  the 
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vehicle  against  extreme  heating  rates  will  be  presented  first.  Then,  the  manoeuvr¬ 
ing  capability  or  'landing  footprint'  available  to  a  vehicle  will  be  discussed. 

1.2  .Types  of  -Protection  Against  Heating 

There  are  two  basic  methods  of  withstanding  reentry  heating. 
Ablatiqn  materials  have  successfully  been  used  on  both  ICBM  nose  cones  and  the 
current  series  of  manned  capsules  in  addition  to  other  classes  of  vehicles. 

Very  broadly  speaking,  the  basic  principle  is  to  use  a  material  which  has  a 
large  heat  of  sublimation  in  order  to  absorb  as  large  a  total  heat  load  as 
possible  for  a  given  mass.  The  actual  heating  rates  experienced  during  the  re¬ 
entry  are-  usually  unimportant  for  such  a  design.  However,  the  maximum  heat 
load  is  very  strictly  defined.  When  exceeded,  the  heat  shield  will  fail  and 
the  spacecraft  will  disintegrate. 

On  the  other  hand,  a  radiative  material  can  withstand  much  higher 
total  heat  loads  as  long  as  the  maximum  heating  rate  does  not  exceed  the  capa¬ 
bility  of  the  material.  The  permissible  heating  rate  is  primarily  a  function 
of  the  temperature  versus  time  capability  of  the  material  (see  Fig.  l),  its 
emissivity,  and  (for  higher  than  average  heating  rates  of  a  short  duration)  its 
heat  capacity.  If  the  material  is  thin  enough,  its  heat  capacity  may  be  neglect¬ 
ed  and  the  temperature  then  assumes  the  radiation  -  equilibrium  value.  In  this 
situation,  the  rate  at  which  heat  is  absorbed  locally  by  the  skin  of  the  space¬ 
craft  equals  the  rate  at  which  heat  is  being  radiated  from  its  surface.  Changes 
in  the  heating  rate  are  then  quickly  reflected  by  a  change  in  the  skin  temperature. 
A  higher  rate  requires  a  higher  temperature  to  radiate  more  heat  and  vice  versa. 

In  general*, ,  the  position  at  which  the  heating  rate  is  most  severe  is  the  nose 
stagnation  point  of  the  vehicle,  since  the  nose  cap  has  the  smallest  effective 
radius  of  curvature .  Therefore,  the  nose  stagnation  point  will  have  the  highest 
temperature  on  the  spacecraft.  Thus,  if  a  suitable  material  can  be  found  for 
the  nose  stagnation  point,  it  can  safely  be  assumed  that  the  rest  of  the  vehicle 
can  also  be  adequately  protected.  Consequently,  unless  explicitly  stated,  the 
radiatipn-equilibrium  nose  stagnation  point  temperature  will  be  the  value  to 
which  reference  is  being  made  whenever  temperature  is  referred  to  in  this  paper. 

•  Two  other  points  which  should  be  mentioned  in  comparing  abl  ative 
and  radiative  cooling.  The  materials  which  employ  the  latter  method  must  always 
be  replaced  after  use.  On  the  other  hand,  while  recoating  may  be  necessary  for 

a  severe  reentry,  it  should  be  possible  to  reuse  in  a  conventional  manner  a 

vehicle  which  employs  radiation  cooling.  It  should  be  noted,  however,  that 
ablative  materials  are  able  to  handle  very  large  heating  rates  (until  the  maxi¬ 
mum  heat  load  is  reached)  that  are  usually  far  larger  than  the  heating  rates 

which  the  maximum  temperature  of  a  radiation  cooled  material  can  support. 

1.3  Footprint  Terminology 

In  order  to  initiate  a  return  to  the  Earth's  surface,  an  orbiting 
spacecraft  will  fire  its  retro  rockets.  Once  this  deboost  phase  is  over,  the 
vehicle  becomes  a  hypersonic  glider  which  is  capable  of  influencing  the  position 
of  its  landing  site  only  by  means  of  aerodynamic  force  controlled  by  the  attitude 
of  the  vehicle.  A  representative  manoeuvre  boundary  or  footprint  from  Ref.  1  is 
shown  in  Fig.  2.  Note  that  the  initial  position  is  at  the  altitude  at  which 
aerodynamic  forces  begin  to  affect  the  path  of  the  reentry  vehicle.  Throughout 
this  investigation,  the  computations  were  started  when  the  perpendicular  dis¬ 
tance  to  the  Earth's  surface  was  350,000  feet.  As  a  consequence,  the  reentry 
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values  of  velocity  and  glide  angle  at  the  initial  position  were  treated  inde¬ 
pendently  of  the  orbital  altitude  at  which  retro- fire  takes  place.  The  actual 
flight  path  is  effectively  separated  for  computational  purposes  from  the  point 
and  conditions  of  retro-fire.  The  only  restriction  is  that  the  deboost  cannot 
occur  at  an  altitude  that  is  greater  than  the  maximum  height  which  is  compatible 
with  conditions  at  the  initial  position  • . 

In  Fig.  3,  additional  terminology  is  given  regarding  the  pro¬ 
jection  of  the  footprint  onto  an  imaginary  sphere  which  is  used  in  place  of 
the  Earth's  surface.  Geocentric  spherical  trigonometry  is  used  to  define  the 
geographical  position  over  the  imaginary  sphere.  Since  the  footprint  boundary 
is  given  for  an  altitude  of  about  100,000  feet,  the  usual  conversion  problems 
of  positioning  the  vehicle  in  terms  of  its  geodetic  latitude  and  longitude  and 
of  determining  the  distance  travelled  over  an  oblate  Earth  can  be  avoided. 

The  actual  method  of  presenting  the  results  requires  some  explana¬ 
tion.  Assume  for  the  moment  that  the  initial  position  of  the  vehicle  is  in  the 
Earth's  equatorial  plane  and  that  the  initial  azimuth  is  either  due  East  or  West. 
Note  that  the  initial  glide  angle  is  quite  immaterial  to  the  range  definition 
(although  it  is  of  great  importance  to  the  pilot  since  too  large  a  value  would 
cause  the  spacecraft  to  disintegrate  through  overheating) .  In  this  situation, 
the  equatorial  plane,  the  great  circle  projection,  the  orbital  track  projection, 
and  the  footprint  midline  all  coincide.  The  'down  range'  is  the  absolute  value 
of  the  longitude  change  from  the  initial  position  (by  definition,  down  range  is 
positive).  The  cross  range  is  the  final  latitude  (cross  range  may  be  positive 
or  negative)  if  the  initial  azimuth  was  due  East  and  the  negative  of  the  final 
latitude  if  the  initial  azimuth  was  due  West.  Note  that  the  definition  must  be 
consistant  with  the  reentry  trajectory  track  over  the  Earth's  surface,  not 
with  the  Earth's  co-ordinate  system.  A  third  range  parameter,  the  azimuth  or 
heading  angle  change,  may  also  be  simply  defined  in  this  situation  as  the  final 
azimuth  minus  (or  plus  in  the  case  of  a  due  West  reentry)  ninety  degrees. 

At  the  initial  position  (which  has  been  arbitrarily  defined  to 
be  the  point  over  the  Earth's  surface  when  the  vehicle's  altitude  is  350,000 
feet),  the  initial  values  of  speed,  altitude,  glide  angle,  azimuth,  latitude, 
and  longitude  define  the  path  which  the  vehicle  will  travel  over  the  Earth' s 
surface.  For  a  non-rotating,  spherical  Earth,  the  final  pattern  or  shape  of 
the  projection  of  the  path  on  the  Earth's  surface  will  not  be  affected  by  the 
initial  values  of  the  last  three  parameters.  The  resultant  equivalent  down 
range,  cross  range  and  azimuth  change  must  be  the  same  as  for  the  case  of  the 
due  East  equatorial  reentry  with  the  same  initial  values  of  speed,  altitude, 
and  glide  angle.  For  this  situation,  the  great  circle  projection,  the  orbital 
track  projection,  and  the  footprint  midline  will  still  coincide  with  each 
other,  but  they  will  no  longer  coincide  with  the  equatorial  plane.  However, 
for  a  rotating, oblate  Earth,  all  three  will  tend  to  diverge.  In  this  situation, 
the  great  circle  projection  remains  in  the  same  relative  position.  Consequently, 
all  range  parameters  are  measured  in  equivalent  degrees  over  the  imaginary 
sphere  relative  to  the  great  circle  projection  through  the  initial  position. 
Figure  3  illustrates  a  typical  footprint  over  a  rotating,  oblate  Earth.  The 
subroutine  which  computes  these  range  values  is  listed  in  Appendix  C  under 
the  name  of  SUBROUTINE  KJUGE. 


3 


1.4  Choice  of  Vehicle  and  Basic  Path 

The  type  of  vehicle  chosen  can  be  seen  in  Fig.  4.  It  satisfies 
the,  ^requirement  of  moderate  lift  to  drag  ratio  since  (L/D)max  can  range  from 
1.3'  to,  1.7  for  this,  type  of  configuration. 

With  regard  to  the  choice  of  the  basic  path,  the  manoeuvrability 
requirement  is  first  used  to  determine  the  method  of  preventing  the  spacecraft 
from  burning  up  during  the  reentry.  When  the  vehicle  must  reach  landing  zones 
corresponding  to  iarge  down  range  distances,  prolonged  flight  times  must  result. 
Consequently  the  total  heat  input  to  an  ablative  material  would  be  extremely 
large  and  might  become  prohibitive.  However,  in  a  design  which  radiates  most 
of  the  .heat  away  from  the  vehicle,  the  total  heat  load  is  not  an  important  con¬ 
sideration.  Therefore,  a  wide  variation  in  the  time  interval  over  which  the 
heating  occurs  will  have  a  negligible  effect  on  the  vehicle  design  as  long  as 
the  peak  temperature  experienced  during  the  flight  is  independent  of  the  point 
at  which  the  spacecraft  lands.  Consequently,  radiation  cooling  offers  the  only 
practical  solution  for  protecting  the  vehicle  against  the  extreme  heating  which 
is  encountered  during  the  reentry.  This  method  then  requires  a  trajectory  that 
satisfies  the  requirement  on  the  peak  temperature. 

An  equilibrium  glide  was  chosen  as  the  basic  path.  This  tra-, 
jectory  is  characterized  (see  Appendix  A  for  a  more  comprehensive  discussion) 
by  the  absence  of  any  long  period  or  'phugoid'  oscillations  along  the  trajectory. 
This  situation  results  from  the  equality  between  the  lift  force  and  the  differ¬ 
ence  between  gravitational  and  centrifugal  forces.  There  were  a  number  of  im¬ 
portant  reasons  for  the  choice  of  the  equilibrium  glide.  First,  it  satisfies 
the  requirement  that  it  be  a  simple  trajectory  since  it  is  a  fixed- attitude 
reentry.  Second,  it  offers  the  lowest  possible  value  for  the  peak  temperature 
which  is  experienced  by  the  vehicle  during  its  reentry  without  the  introduction 
of  a  sophisticated  optimal  programming  procedure.  Third,  the  equilibrium  glide 
also  satisfies  the  requirement  that  departures  from  this  basic  path  may  easily 
be  made  in  order  to  provide  the  manoeuvring  capability  necessary  to  land  at 
sites  other  than  the  basic  terminal  position. 

The  reader  should  also  note  the  simple  analytical  result  for  the 
peak  temperature  which  is  experienced  during  an  equilibrium  glide  (see  Appendix 
A).  The  ease  of  obtaining  this  comparison  standard  and  its  simplicity  of  de¬ 
finition  will  greatly  outweigh  the  benefits  of  obtaining  any  reduction  in  the 
standard  peak  temperature  by  complex  optimal  procedures. 


1.5  CONTAC  System 

The  problem,  then,  is  to  conceive  of  as  simple  a  procedure  as 
possible  which  will  allow  the  next  generation  of  reentry  vehicles  the  freedom 
to  manoeuvre  to  their  desired  landing  zones.  Ideally,  implementation  of  this 
procedure  will  not  require  any  additional  protection  for  the  pilot  or  the  vehicle. 

As  a  solution  to  this  problem,  the  author  proposes  a  CONstant 
Temperature  Attitude  Control  (CONTAC)  system  which  will  limit  the  peak  tempera¬ 
ture  which  is  experienced  by  the  vehicle  during  a  reentry  to  a  value  which  is 
only  moderately  (by  20F°  or  30F°)  greater  than  the  peak  temperature  in  a  stand¬ 
ard  equilibrium  glide.  As  a  consequence,  the  CONTAC..  system  will  operate  within 
the  constraint  that  a  negligible  weight  penalty  is  incurred  relative  to  an 
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equilibrium  glide,  which  it  is  assumed  the  vehicle  can  safely  fly.  A  further 
constraint  will  be  that  neither  predictive  nor  memory  logic  (with  a  minor 
exception)  will  be  used  by  the  CONTAC  system.  This  requirement  means  that  the 
CONTAC  system  must  produce  the  solution  on  a  real  time  basis  when  operating  in 
the  actual  reentry  environment. 

It  will  be  shown  that  the  CONTAC  system  gives  the  pilot  the  ex¬ 
tensive  manoeuvrability  which  he  requires  in  order  to  land  the  spacecraft  at 
any  point  under  or  near  the  orbital  envelope  (Fig.  5)*  When  coupled  with  the 
proper  range  controller,  the  complete  system  will  provide  automatic  guidance 
to  the  final  destination  and  prevent  any  dangerous  temperature  peaks  which  might 
destroy  the  vehicle.  In  addition,  no  undue  waiting  periods  will  be  required. 

For  most  landing  areas  within  the  orbital  envelope,  the  delay  before  retro- 
fire  will  be  less  than  twelve  hours  and  no  point  will  require  more  than  twenty- 
four  hours.  For  orbits  of  small  inclination  to  the  equatorial  plane  the  waiting 
times  would  be  usually  less  than  one  orbital  period.  Furthermore,  the  peak 
temperature  will  always  be  very  near  to  the  peak  value  encountered  during  an 
equilibrium  glide  except  in  those  cases  for  which  the  atmospheric  density  varia¬ 
tions  about  the  mean  are  both  large  and  random.  Under  these  extreme  conditions, 
it  will  be  impossible  to  prevent  a  moderate  rise  in  the  peak  temperature  whenever 
the  CONTAC  system  encounters  a  sequence  of  density  variations  for  which  no  further 
correction  by  means  of  attitude  control  is  available  (i.e.  when  the  optimum 
attitude  -  optimum  angle  of  attack  and  zero  roll  angle  -  has  already  been 
assumed) .  Since  adverse  conditions  which  occur  randomly  are  impossible  to 
predict,  a  statistical  result  must  be  given  and  vehicle  protection  provided  for 
the  worst  sequence  of  events  which  might  occur  with  a  certain  design  probability. 


1.6  Comparison  with  Other  Systems 

As  detailed  in  Ref.  1,  a  number  of  authors  have  studied  the  pro¬ 
blem  of  guiding  a  spacecraft  to  its  final  destination.  However,  this  investi¬ 
gation  is  not  primarily  concerned  with  the  actual  guidance,  but  rather  with 
the  means  by  which  the  heating  problems  may  be  minimized.  In  addition,  the 
effects  of  Earth  rotation  and  oblateness  upon  both  the  heating  and  the  guidance 
problem  will  be  ascertained.  Finally,  a  test  to  determine  the  effect  of  random 
variations  in  atmospheric  density  will  be  presented. 

With  regard  to  the  heating  problem,  a  solution  for  minimizing 
the  total  heat  load  can  be  found  in  Refs.  2  and  3*  However,  since  these 
solutions  must  operate  within  the  inherent  limitation  imposed  by  the  maximum 
allowable  heat  load,  the  manoeuvring  capability  is  restricted  to  landing  sites 
which  are  relatively  close  to  the  initial  position. 

The  only  other  solution  for  minimizing  peak  temperatures  that 
has  been  found  in  the  open  literature  is  a  Temperature  Rate  Flight  Control 
System  which  has  been  developed  by  Stalony-Dobrzanski  (Ref.  4).  The  results 
show  that  the  TRFCS  is  also  able  to  prevent  the  peak  temperature  which  occurs 
during  a  reentry  from  becoming  greater  than  the  value  which  is  experienced  dur¬ 
ing  an  equilibrium  glide.  (The  use  of  the  equilibrium  glide  as  the  basic  path 
for  the  TRFCS  is  an  additional  reason  for  its  choice  as  the  standard  of  com¬ 
parison  in  the  present  investigation.)  However,  in  order  to  implement  the  TRFCS, 
the  derivative  of  the  temperature  is  required  along  with  an  elaborate  analogue 
system  to  analyze  the  rate  of  temperature  change.  The  attitude  of  the  vehicle 
is  adjusted  using  either  the  angle  of  attack  or  the  bank  angle  as  the  primary 
control  variable. 
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In  contrast,  the  CONTAC  system  uses  a  simple  digital  program  to 
analyze  directly  the  actual  temperature.  In  addition,  the  attitude  is  controlled 
via  pitch  and  roll  channels  simultaneously. 

With  regard  to  the  effects  on  the  Earth's  rotation  and  oblateness, 
there  are  no  papers  available  which  analyze  the  resulting  changes  in  the  foot¬ 
print  when  constraints  upon  the  peak  temperature  (or  for  that  matter  maximum 
heat  load)  are  present.  There  have  been  a  few  studies  which  analyzed  the  effects 
of  the.  Earth's  rotation,  but  no  constraints  in  connection  with  heating  problems 
were  incorporated.  * 

T  Due  to  the  lack  of  significant  statistical  data,  the  influence 

density  variations  in  the  Earth's  atmosphere  is  still  extremely  difficult  to 
analyse.  A  number  of  average  curves  of  density  versus  altitude  are  available 
for  both  summer  and  winter  for  latitudes  from  the  equator  to  75°N.  Based  on 
this  type  of  information,  a  number  of  papers  have  appeared  which  deal  with  the 
effect  of  a  constant  percentage  change  in  the  density  upon  the  guidance  problem. 
However,  even  if  heating  effects  were  included  in  these  investigations,  the 
results  cannot  take  account  of  the  fact  that  the  atmosphere  is  a  dynamic  system 
whose  state  can  only  be  described  in  a  statistical  manner.  With  regard  to  the 
effects  of  random  density  variations  about  an  average  curve  of  density  versus 
altitude,  the  author  is  not  aware  of  any  other  studies  which  deal  with  the  heat¬ 
ing  problem  either  from  an  optimal  point  of  view  (such  as  those  represented  by 
Refs.  2  and  3)  or  from>a system  approach  (such  as  proposed  by  Stalony-Dobrzanski) . 
Indeed,  from  a  detailed  examination  of  the  way  in  which  optimal  procedures  are 
applied,  the  inclusion  of  random  density  variations  would  seem  to  rule  out  this 
concept.  The  iteration  procedure  used  depends  on  accurate  external  environ¬ 
mental  specifications.  If  the  density  is  a  random  quantity,  it  would  not  be 
possible  under  the  presently  used  procedures  to  arrive  at  a  final  solution. 

Only  a  real  time  approach  inherent  in  both  the  TRFCS  and  the  CONTAC  system 
can  attempt  to  provide  a  solution.  The  present  investigation  will  include  the 
results  of  the  attempts  by  the  CONTAC  system  to  deal  with  random  density  varia¬ 
tions  as  well  as  constant  percentage  changes  in  the  density. 

1.7  Trajectory  Phases 

In  order  to  simplify  the  description  and  analysis  of  the  various 
aspects  of  the  CONTAC  system,  the  reentry  trajectory  has  been  divided  into  five 
rather  distinct  phases.  Each  phase  is  described  rather  briefly  at  this  point 
in  order  to  acquaint  the  reader  with  the  overall  concept.  What  the  CONTAC 
system  is  capable  of  doing  is  also  explained.  In  addition,  problems  which  it 
cannot  solve  are  mentioned  and  some  suggestions  are  made  regarding  the  solution 
of  these  difficulties. 

It  is  assumed  that  the  vehicle  is  under  the  primary  control  of  a 
navigation  and  guidance  system,  that  monitors  and  adjusts  the  two  primary 
attitude  angles,  a  and  q>  ,  in  order  to  home  on  the  desired  landing  site.  At 
certain  short  but  critical  periods  of  the  flight,  the  CONTAC  system,  as  a  result 
of  monitoring  vehicle  temperature,  automatically  cuts  in  and  overrides  the 
navigation  system.  This  situation  occurs  only  during  the  portions  of  Phases  III 
and  IV  of  the  reentry  when .the  CONTAC  system  is  in  active  operation.  There  are, 
of  course,  additional  constraints  on  the  attitude  angles.  These  will  be  pre¬ 
sented  during  the  discussion  of  each  phase. 
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The  five  phases  of  the  trajectory  are  shown  in  Fig.  6.  Although 
no  such  sharply  defined  divisions  will  occur  during  an  actual  flight,  it  is 
simpler  to  divide  the  reentry  into  the  five  parts  shown  for  the  purpose  of 
explaining  what  takes  place  within  each  phase. 

The  € OUT AC  system  has  been  designed  for  space  vehicles  entering 
the  Earth's  atmosphere  from  altitudes  up  to  five  hundred  miles.  It  is  conceiv¬ 
able  that  some  application  might  be  found  for  reentries  from  positions  higher 
than  this  altitude .  However,  the  primary  requirement  is  that  the  maximum  con¬ 
vective  heating  rate  first  experienced  during  the  reentry  (Fig.  6)  must  not  pro¬ 
duce  a  temperature  greater  than  the  peak  value  which  would  be  experienced  during 
an  equilibrium  glide.  If  this  constraint  is  satisfied,  the  COIfTAC  system  is 
capable  of  preventing  all  subsequent  maximum  temperatures  from  becoming  very 
much* greater  than  the  peak  experienced  during  an  equilibrium  glide.  But,  if  the 
deboost  altitude  is  above  five  hundred  miles,  the  spacecraft  must  reenter  the 
atmosphere  at  a  speed  which  is  substantially  greater  than  the  circular  value. 

In  order  to  prevent  skipout,  the  reentry  glide  angle  must  also  be  increased.  As 
a  result,  the  first  maximum  experienced  for  the  convective  heating  rate  exceeds 
the  peak  for  an  equilibrium  glide.  Unless  this  single  pulse  can  be  tolerated, 
the  vehicle  fails. 

Therefore,  the  CONTAC  system  is  intended  to  be  used  for  reentries 
from  low  orbital  altitudes ,  up  to  about  five  hundred  miles,  with  vehicles  that 
use  radiative  cooling  as  the  primary  protection  against  failure  due  to  heating. 
Examples  of  this  type  of  operation  include  shuttle  service  for  space  stations, 
low  orbit  reconnaissance,  and  low  orbit  satellite  service  and  repair. 

1.8  Phase  I  -  Pre-Atmospheric  Flight 

It  should  be  noted  at  this  point  that  although  the  CONTAC  system 
may  go  into  operation  as  soon  as  the  retro  rockets  have  been  fired  (or  some 
other  method  is  used  to  push  the  spacecraft  out  of  its  vacuum  orbit),  the  angle 
of  attack  of  the  vehicle  is  quite  immaterial  for  a  period  of  between  ten  and 
thirty  minutes,  depending  on  the  height  o*’  orbital  ejection.  This  delay  is  due 
to  the  fact  that  aerodynamic  forces  do  not  become  appreciable  until  an  altitude 
of  about  350,000  feet  is  reached  (Fig. 7)-  Thus  the  pilot  will  have  ample  time 
to  jettison  any  service  modules  or  extraneous  externals  and  to  perform  the  tasks 
associated  with  the  alignment. of  the  spacecraft  prior  to  the  buildup  of  atmos¬ 
pheric  density. 


One  point  of  importance  which  is  associated  with  this  phase  is 
the  deboost  attitude  alignment  and  velocity  increment.  Both  values  must  be 
carefully  computed  and  equally  carefully  implemented  since  a  large  error  in 
either  could  cause  excessive  heating  rates  leading  to  the  failure  of  the  space¬ 
craft.  Figure  8  shows  the  effect  of  errors  in  retro-fire  attitude  alignment  and 
velocity  increment.  However,  past  experience  in  the  Gemini  program  indicates 
that  the  require  accuracies  for  attitude  alignment  and  velocity  increment  will 
not  be  difficult  to  achieve.  Interpretation  of  the  results  indicates  that  the 
Gemini  retro-fire  system  experienced  maximum  errors  of  one  degree  in  attitude 
alignment  and  three  feet  per  second  in  velocity  increment  (Ref.  5 » •  A  simul¬ 
taneous  maximum  error  in  both  would  yield  a  maximum  resultant  error  at  the  end 
of  Phase  I  at  350,000  feet  of  0.1  degrees  in  flight  path  angle,  10  feet  per 
second  in  speed  and  2  degrees  in  down  range.  Cross  range  might  also  be  affected 
but  to  a  negligible  extent;  the  maximum  error  could  amount  to  only  0.2  degrees. 
Note  that  all  of  the  above  nominal  error  values  are  maximum  and  would  be  very 
unlikely  to  occur  simultaneously. 
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Since  experience  has  proven  that  a  highly  accurate  calculation 
of  the  alignment  attitude  and  velocity  increment  values  is  possible,  there  should 
be  no  difficulty  in  this  area.  However,  it  may  be  necessary  to  change  the  method 
by  which  the  calculation  is  made.  Consequently,  the  allowable  margin  of  error 
is  described  and  discussed  in  Section  VI  where  the  effects  of  incorrect  reentry 
conditions  at  the  beginning  of  Phase  II  .-are  considered. 

1.9  Phase  II  -  Overall  Down  Range  Control 

For  the  purposes  of  this  investigation,  the  reentry  begins  at 
350,000  feet.  This  is  due  to  the  fact  that  for  all  practical  purposes,  the 
atmosphere  has  no  effect  on  the  flight  path  until  the  spacecraft  reaches  this 
altitude.  Since  the  pilot  has  no  control  over  the  trajectory  between  the  time 
he. fires  the  retro  rockets  and  the  time  he  descends  to  this  height,  it  is  best 
to  separate  the  problems  associated  with  achieving  the  correct  reentry  conditions 
at  350,000  feet  (by  the  proper  attitude  alignment  and  retro  velocity  increment 
at  orbital  altitudes)  from  the  problems  associated  with  controlling  the  vehicle's 
down  range  (by  the  proper  angle  of  attack  once  lift  and  drag  become  appreciable) . 
In  fact,  in  order  to  simplify  the  problem  even  further,  the  portion  of  the  flight 
within  the  atmosphere  has  been  sub-divided  into  four  parts.  Phases  II,  III,  IV. 
and  V  of  the  reentry. 

During  Phase  II,  the  CONTAC  system  must  be  in  operation;  however, 
it  exercises  no  actual  control  except  to  signal  the  end  of  Phase  II  at  the  point 
when  the  first  qm«y  is  reached  (Fig.  9).  At  this  position,  the  CONTAC  system 
switches  from  passive  to  active  operation  and  Phase  III  begins.  Therefore,  during 
Phase  II,  the  pilot  (or  the  automatic  guidance  system)  has  the  function  of  devis¬ 
ing:’.  the  angle  of  attack  program  that  will  eliminate  down  range  errors  which 
may  initially  be  present  or  which  may  occur  during  Phase  II  due  to  the  Earth's 
departure  from  a  circular,  non-rotating  sphere  with  a  standard  time  invariant 
atmosphere.  Since  Phase  II  lasts  for  about  four  minutes,  there  is  ample  time 
for  any  down  range  corrections  which  are  required  in  order  to  ensure  that  the 
vehicle  lands  on  target. 

To  some  extent,  it  will  be  possible  to  compensate,  by  adjustment 
of  the  initial  retro-fire,  for  down  range  errors  which  it  is  anticipated  will 
occur  during  Phase  II  due.to  some  of  the  above  non- ideal  conditions.  As  a  first 
approximation,  a  constant  angle  of  attack  can  be  used  during  the  calculation  for 
Phase  II.  In  this  situation,  down  range  errors  which  will  occur  due  to  rota¬ 
tion  and  oblateness  can  be  evaluated  rather  easily  for  this  portion  of  the 
flight.  The  same  calculation  procedure  can  also  be  used  to  find  the  down  range 
errors  which  will  occur  during  Phase  I.  Thus,  at  the  deb'oost  position,  it  should 
be  feasible  to  make  the  necessary  corrections  for  down  range  errors  which  will 
occur  during  Phases  I  and  II  due  to  Earth  rotation  and  oblateness.  Note  that 
the  time  of  retro-fire  may  also  be  varied  in  addition  to  the  attitude  align¬ 
ment  and  velocity  increment.  This  additional  parameter  will  probably  provide 
most  of  the  necessary  correction. 

Errors  in  down  range  which  occur  during  Phase  II  due  to  variations 
in  the  atmospheric  density  are  much  more  difficult  to  evaluate.  Since  it  is 
unlikely  that  an  accurate  prediction  of  the  density  can.be  made,  an  estimate 
based  on  the  available  statistical  data  should  be  made.  It  may  then  be  possible 
to  compensate  for  the  anticipated  down  range  error  by  adjusting  the  retro-fire 
parameters. 


8 


mii.hiuiim  iim  . . . . . . -  - - ■„  - - — .  .  IIMJI  jr..— m  n  m  T~i.."m;'‘HOTa 


The  only  restriction  on  the  angle  of  attack  during  Phase  II  (angle 
of  roll  is  zero)  is  that  it  must  not  be  so  low  as  to  cause  the  first  qm»y  to 
exceed  the  maximum  allowable  value  nor  so  high  as  to  endanger  the  rear  portion 
of  the  vehicle  from  excessive  heating  rates.  In  general,  the  optimum  value 
lies  very  near  the  angle  of  attack  which  yields  the  maximum  lift  coefficient. 
Larger  values  for  the  angle  of  attack  will  reduce  the  down  range  and  smaller 
values  will  increase  it.  For  the  vehicle  studied,  the  optimum  angle  of  attack 
was  55  degrees.  The  largest  angle  of  attack  used  during  Phase  II  was  65  degrees 
and  the  smallest  value  was  35  degrees.  These  two  limits  yielded  a  down  range 
interval  of  about  75  degrees  or  about  5000  nautical  miles. 

It  should  be  noted  that  almost  all  of  the  incremental  down  range 
achieved  by  modulating  the  angle  of  attack  during  Phase  II  is  not  apparent  until 
the  end  of  Phase  III.  This  is  because  the  primary  effect  of  varying  the  entry 
angle  of  attack  is  a  variation  in  the  vehicle’s  speed  at  the  end  of  Phase  II 
without  any  major  difference  in  altitude,  glide  angle,  or  down  range  distance 
(Fig.  9)*  A  comprehensive  discussion  and  analysis  of  this  point  is  given  in 
Sec.  6.2. 


1.10  Phase  III  -  Equilibrium  Glide  Attainment 

The  most  critical  portion  of  the  reentry  is  Phase  III,  the  various 
sections  of  which  are  shoun  in  Fig. 10.  During  this  part  of  the  flight,  the 
CORTAC  system  must  remove  energy  from  the  vehicle  at  a  rate  designed  to  place 
the  vehicle  on  or  very  near  to  an  equilibrium  glide  trajectory.  However,  the 
rate  of  energy  dissipation  must  also  be  flexible  enough  to  allow  the  pilot  (or 
the  automatic  guidance  system)  to  control  the  vehicle’s  range  in  order  to  counter¬ 
act  all  departures  from  non- ideal  conditions  as  well  as  any  errors  present  at 
the  start  of  Phase  III.  Failure  Of  the  CONIAC-guidance  combination  to  cope 
with  adverse  conditions  will  result  in  the  vehicle  over-heating  during  the 
-  flight  or  landing  at  the  wrong  location,  either  of  which  could  be  disastrous. 

The  actual  details  of  the  COKTAC  system  are  given  in  Section  II. 

At  this  point,  it  is  sufficient  to  state  that  the  phugoid  or  long  period 
oscillation  can  be  completely  eliminated  before  the  peak  temperature  is  ex¬ 
perienced  and  as  a  result,  that  value  is  very  nearly  minimized. 

For  most  reentries,  the  attitude  which  is  adopted  subsequent  to 
the  first  q^y  is  held  for  between  two  and  four  minutes.  During  this  period, 
the  pilot  again  has  an  opportunity  to  continue  correcting  any  down  range  errors 
as  well  as  beginning  to  reduce  any  cross  range  discrepancy.  The  additional 
incremental  down  range  available  during  Phase  III  varies  from  about  1500  miles 
when  the  entry  angle  of  attack  is  65  degrees  up  to  about  5000  miles  when  the 
entry  angle  of  attack  is  35  degrees.  It  should  be  repeated  that  practically 
all  of  the  down  range  increment  available  due  to  Phase  II  attitude  control 
actually  occurs  during  Phase  III  due  to  the  variation  in  speed  at  the  beginning 
of  Phase  III.  The  pilot  or  the  guidance  system  oust  be  aware  of  and  be  able 
to  compensate  for  this  speed  difference. 

During  Phase  III,  a  number  of  pairs  of  values  are  used  for  the 
upper  and  lower  limits  on  the  attitude  angles  of  the  vehicle.  These  pairs  of 
values  are  inputs  to  the  CORIAC  system  and  may  be  regarded  as  constants  within 
this  system.  During  an  actual  or  simulated  reentry  in  which  a  pilot  or  range 
controller  is  available,  the  attitude  limits  will  be  changed  to  provide  the 
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necessary  range  control  in  order  to  land  the  spacecraft  at  the  proper  point 
on  the  Earth’s  surface.  It.  is  a  feature  of  the  CONTAC  system  that  it  is 
unaffected  by  the  values  for  these  limits  as  long  as  the  upper  and  lower  values 
of.  certain  key  pairs  do  not  coincide.  As  long  as  this  requirement  is  satisfied, 
the* CONTAC  system  can -damp  .the  long  period  oscillations  from  the  trajectory  and 
place  the  vehicle  on  the  equilibrium  glide  before  the  peak  temperature  occurs. 

In  very  simple  terms,  the  pilot  will  select  the  upper  and  lower 
limits  for  the  angles  of  attack  and  roll  which  best  fit  the  choice  of  the  land¬ 
ing-area  toward  which  he  desires  to  proceed,  although  certain  practical  con¬ 
siderations  will  restrict  the  pilot's  freedom  of  choice  (the  upper  and  lower 
limits  must  not  coincide  for  both  pitch  and  roll) .  The  CONTAC  system  then  auto¬ 
matically  adjusts  the  attitude  angles  (or  the  pilot  himself  if  the  associated 
control  hardware  should  fail,  but  not  the  temperature  readout  equipment)  to 
some  value  between  the  upper;. and  lower  limits  of  the  angles  of  pitch  and  roll. 

It  has  definitely  been  determined  that  the  CONTAC  system  will  hold 
the  temperature  at  very  nearly  a  constant  value  as  long  as  the  vehicle’s  attitude 
has  not  reached  either  an  upper  or  lo*er  limit.  A  single  skip  usually  occurs 
just  after  the  first  but  this  is  quite  intentional,  and  is  designed  to 

provide  down  range  control  capability  during  the  early  part  of  the  reentry. 
Indeed,  in  some  situations,  this  short  skip  is  completely  eliminated  in  order 
t ©substantially  reduce  the  down  range. 

Thus,  by  the  time  the  vehicle  descends  to  about  260,000  feet,  it 
has  been  placed  upon  an  equilibrium  gi.ide  corresponding  to  the  attitude  at  which 
the  spacecraft  is  currently  being  held.  The  lowest  value  for  the  peak  tempera¬ 
ture  is  achieved  by  using  the  optimum  attitude.  It  oecuvs  at  maximum  vertical 
force  (see  Appendix  A)  which  is  at  55  degrees  angle  of  attack  and  zero  roll 
angle.  This  optimum  attitude  is  then  held  by  the  CONIAC  system  quite  auto¬ 
matically  until  the  peak  temperature  is  attained.  At  that  point,  Fnase  IV 
begins  and  the  attitude  is  changed  to  provide  down  range  and  cross  range 
manoeuvring  capability.  In  practically  all  the  re-entries  studied.  Phase  IV 
began  at  an  altitude  of  about  245,000  feet.  However,  in  cases  where  the  density 
variations  were  quite  large,  Phase  IV  began  at  an  altitude  as  much  as  30,000 
feet  lower  than  the  above  value. 

1.11  Phase  IV  -  Attitude  Control  Release 

There  will,  of  course,  be  problems  associated  with  unpredictable 
increases  in  air  density  which  may  drive  the  temperature  during  this  phase  be¬ 
yond  its  previous  maximum.  In  such  a  situation,  it  might  be  expected  that  the 
CONTAC  system  would  have  to  revert  to  Phase  III  operation.  However,  the  reader 
may  recall  that  no  such  sharply  defined  division  occurs  during  the  actual  flight. 
Therefore,  since  the  CONTAC  system  has  neither  memory  nor  predictive  logic,  it 
is  impossible  for  any  such  division  to  exist  in  the  controlling  concept.  When 
in  operation,  the  CONTAC  system  is  attempting  to  keep  the  temperature  constant, 
and  if  an  increase  occurs,  whatever  the  cause,  the  CONTAC  system  will  change 
the  attitude  so  as  to  reduce  the  temperature.  If  the  attitude  is  already  at  its 
optimum  value  and  the  CONTAC  system  is  unable  to  prevent  a  further  increase  in 
the  temperature,  then  the  conclusion  to  be  drawn  is  quite  obvious,  i.e.  Phase  III 
had  not  yet  been  completed.  But  the  reversion  is  only  with  respect  to  the 
author's  arbitrary  division  of  the  flight.  The  CONTAC  system  is  not  affected. 


10 


In  fact,  the  pilot  (cr  the  automatic  guidance  system)  need  not 
in  principle  even  know  that  the  CONTAC  system  exists,  let  alone  hov  it  operates. 
Naturally,  he  would  wonder  why  he  is  unable  to  assume  a  particular  attitude  for 
range  control  during  certain  portions  of  the  reentry.  However,  the  CONTAC 
system  is  completely  automatic  so  that  it  can  cope  with  unpredictable  situations 
about  which  the  pilot  may  be  unaware. 

The  pilot  would,  of  course,  actually  be  informed  about  both  the 
existence  and  the  operation  of  the  CONTAC  system  and  he  cay,  therefore,  depend 
on  it  to  properly  reduce  the  attitude  from  the  optimum  position  at  the  beginning 
of  Phase  IV  to  the  desired  values  of  pitch  and  roll  sometime  before  the  end  of 
Phase  IV.  In  addition,  with  regard  to  where  Phase  III  ends  and  Phase  IV  begins, 
the  pilot  can  leave  this  quite  arbitrary  decision  entirely  to  the  built-in  logic 
of  the  CONTAC  system.  It  is  quite  immaterial,  as  far  as  he  is  concerned,  just 
where  the  peak  temperature  occurs  during  the  i light.  The  only  problem  is  to 
prevent  that  value  from  becoming  any  larger  than  necessary,  and  this  problem  the 
CONTAC  system  will  solve  all  by  itself. 

Therefore ,  during  Phase  IV,  the  CONTAC  system  will  gradually  re¬ 
duce  the  angles  of  pitch  and  roll  until  at  some  particular  point  in  the  flight, 
both  values  are  under  the  complete  control  of  the  pilot.  The  transfer  of  this 
control  may  be  revoked  at  any  time  after  it  is  first  cade;  however,  after  the 
temperature  has  fallen  several  hundred  degrees,  the  pilot  will  realize  that  the 
transfer  is  complete  and  that  the  CONTAC  systen  is  noilonger  in  operation. 

This  transfer  always  occurs  by  the  time  the  speed  and  height  have  dropped  to 
values  in  the  vicinity  of  12,000  feet  per  second  and  150,000  feet,  respectively. 
At  that  time,  a  procedure  to  increase  the  cross  range  will .have  come  into  effect. 
A  paper  by  Vagner  (Ref.  6)  discusses  the  variation  of  the  bank  angle  of  the 
spacecraft  in  relation  to  cross  range  control.  The  results  indicate  that  large 
angles  of  roll  should  be  used  at  the  beginning  of  the  reentry  when  the  speed  is 
high  and  that  small  angles  of  roll  are  best  near  the  end  of  the  flight  when  the 
speed  is  low.  Various  modes  of  dependence  between  the  angle  of  roll  and  the 
flight  parameters  speed,  altitude,  and  energy,  were  investigated  to  determine 
the  best  relationship.  The  calculations  showed  that  a  simple  linear  variation 
of  bank  angle  versus  speed  gives  the  best  overall  performance. 

1.12  Phase  V  -  Final  Range  Control  and  Touchdown 

As  this  paper  is  being  written,  the  procedures  for  the  final 
landing  cf  a  spacecraft  are  being  developed  in  principle  and  tested  in  practice. 
However,  an  analysis  of  this  phase  is  outside  the  scope  of  this  investigation. 

It  is  therefore,  of  little  value  to  speculate  upon  the  best  method  to  reach  the 
landing  site  and  perform  the  touchdown.  Several  papers  (among  them  Refs.  7,  8, 
and.  9)  have  appeared  which  deal  with  this  manoeuvre- so-  that?  some:,  ideas  are  2  already 
available  concerning  the  problems  which  say  be  encountered  and  the  methods  of 
solution. 


To  this  date,  all  of  the  flight  experience  in  the  supersonic 
regime  has  been  with  the  X-15  vehicle.  Thus,  it  may  be  quite  some  time  before 
the  best  solutions  to  the  various  problems  are  found.  However  a  basic  outline 
of  the  final  phase  has  emerged  from  this  test  program.  The  additional  local  down 
range  which  can  be  attained  is  about  200  miles  for  the  conditions  given  in  Fig. 

3  for  the  footprint  boundary.  But,  it  may  be  noted  that-  the  final  azimuth  is 
approximately  perpendicular  to  the  boundary  lice.  Consequently,  the  footprints 
which  are  presented  in  Sections  VI  ana  VII  are  a  conservative  estimate  of  the 
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vehicle's  capability.  One  point  of  considerable  importance  should  be  noted.  At 
the  end  of  Phase  IV,  the  valuesof  speed  aid  altitude  within  and  on  the  footprint 
boundary  are  independent  of  the  final  position.  Consequently,  the  same  pro¬ 
cedure  can  always  be  used  for  the  approach  during  Phase  V. 


II.  ATTITUDE  CONTROL  LOGIC 
2.1  Introduction 

This  section  discusses  the  detailed  concepts  behind  the  COMTAC 
system.  The  methods  by  which  the  system  may  be  implemented  on  the  spacecraft 
itself  will  also  be  suggested  although  the  actual  details  may  require  modifica¬ 
tions  in  order  to  conform  to  hardware  requirements. 

In  order  to  give  the  reader  some  insight  as  to  why  the  COMAC 
system  is  used,  a  typical  reentry  is  sh  m  in  Fig.  11a  in  which  no  attitude 
control  was  attempted  except  at  the  first  skip  point.  At  this  position,  the 
angle  of  attack  was  reduced  from  the  optimum  value  of  55  degrees  and  zero  roll 
angle  to  22.7  degrees  with  a  roll  angle  of  50  degrees  to  the  left  in  order  to 
obtain  greater  down  range  and  to  produce  some  cross  range.  However,  due  to 
the  reduction  in  the  angle  of  attack,  a  phugoid  oscillation  is  initiated  which 
results  in  a  large  temperature  oscillation  as  well.  Because  the  latter  is 
synchronized  almost  exactly  with  the  phugoid,  it  was  logical  to  -attempt  to  con¬ 
trol  both  the  phugoid  and  the  temperature  oscillations  by  monitoring  only  the 
temperature.  The  reason  for  the  use  of  the  temperature  as  the  controlling 
variable  is  that  this  value  can  be  measured  directly  on  board  the  spacecraft 
with  a  high  degree  of  accuracy.  On  the  other  hand,  the  altitude  is  a  very 
difficult  quantity  to  measure  when  the  vehicle  is  still  250,000  feet  above  the 
Earth's  surface.  In  addition,  the  temperature  oscillations  are  the  ones  which 
must  be  controlled.  It  is  quite  immaterial  whether  or  not  the  phugoid 
oscillation  is  controlled  and  damped  out  as  soon  as  possible  or  is  amplified 
and  made  larger  so  long  as  the  temperature  peaks  are  reduced.  Of  course,  the 
fact  that  the  elimination  of  the  temperature  oscillation  also  results  in  the 
elimination  of  the  phugoid  (Fig.  11b)  is  an  added  bonus  which  results  in  greater 
range  control  over  the  vehicle  by  the  pilot. 

It  may  be  of  interest  to  note  that  attitude  control  based  upon 
altitude  measurements  was  also  investigated  at  an  early  stage  in  the  development 
of  the  U0HTAC  system.  In  view  of  the  results,  which  were  inferior  to  the  C0KTAC 
concept,  and  the  difficulty  associated  with  the  measurement  of  the  altitude  dur¬ 
ing  an  actual  reentry,  all  effort  on  this  aspect  was  abandoned. 

A  discussion  of  the  basic  philosophy  behind  the  C0NTAC  system  will 
provide  a  greater  understanding  of  how  it  works.  The  main  purpose  is  to  dissi¬ 
pate  as  much  of  the  vehicle's  kinetic  energy  at  as  high  an  altitude  as  possible. 
In  order  to  achieve  this  result,  the  C0NTAC  system  is  designed  to  prevent  any 
decrease  in  the  temperature  after  the  first  maximum  is  reached  du3  *  g  the  initial 
pass  into  the  atmosphere.  In  order  to  prevent  the  decrease,  the  angle  of  attack 
is  reduced  which  increases  the  L/D  (actually  the  ratio  of  vertical  force  to  drag, 
but  at  constant  bank  angle,  L/D  is  just  as  appropriate)  and  thereby  extends  the 
glide  path  over  a  longer  distance  at  the  higher  altitudes.  However,  the  angle 
of  attack  eventually  reaches  its  minimum  value  and  at  that  point  the  temperature 
will  begin  to  fall.  However,  when  at  some  later  time  the  temperature  again 
starts  to  rise,  a  decision  must  be  made  as  to  when  is  the  best  time  to  increase 


12 


the  angle  of  attack  in  order  to  prevent  a  subsequent  temperature  peak  from 
exceeding  the  highest  temperature  experienced  prior  to  that  point  in  the  flight. 
With  no  prediction  scheme  available,  the  best  choice  of  the  time  to  increase  the 
angle  of  attack  is  immediately,  i.e.  as  soon  as  the  temperature  again  begins  to 
rise.  Keeping  the  temperature  constant  from  that  point  on  implies  a  gradual 
reduction  in  L/D  which  in  turn  may  result  in  lower  altitudes  sooner  than  is  de¬ 
sirable.  There  may  be  an  optimum  period  of  time  to  wait  before  initiating  any 
increase  in  the  angle  of  attack.  However  the  absence  of  both  memory  and  pre¬ 
dictive  logic  within  the  CONTAC  system  precludes  any  use  of  a  time  delay  interval. 

In  summary,  the  CONTAC  system  operates  on  the  principle  that  the 
best  place  to  lose  speed  is  at  as  high  an  altitude  as  possible.  Therefore,  the 
angle  of  attack  is  gradually  reduced  in  such  a  manner  as  to  prevent  the  tempera¬ 
ture  from  falling  after  it  achieves  a  local  peak  value.  That  action  results  in 
an  extension  of  the  glide  path  at  higher  altitudes  due  to  the  resulting  higher 
values  of  L/D.  On  the  other  hand,  it  was  found  best  (by  increasing  the  angle 
of  attack)  to  prevent  for  as  long  as  possible  any  subsequent  rise  in  temperature 
after  a  minimum  value  occurs.  The  result  is  a  system  which  always  attempts  to 
keep  the  temperature  constant  as  long  as  angle  of  attack  changes  are  possible, 
i.e.  as  long  as  the  angle  of  attack  was  at  neither  its  lower  nor  upper  limit 
following  a  maximum  or  minimum,  respectively,  in  the  equilibrium  nose  stagnation 
point  temperature. 

The  same  principles  are  also  used  to  control  the  bank  angle  in 
conjunction  with  the  angle  of  attack.  The  only  difference  is  that  the  "wings 
level"  position  is  substituted  for  the  upper  limit  since  the  maximum  value  for 
the  ratio  of  vertical  force  to  drag  occurs  for  zero  bank  angle.  On  the  other 
hand,  the  minimum  value  for  this  ratio  occurs  at  the  largest  angle  of  roll. 
Therefore,  the  "lower  limit"  is  the  largest  bank  angle  to  the  right  or  left  de¬ 
pending  on  which  direction  the  landing  zone  is  situated  with  respect  to  the 
orbital  track. 

2.2  Choice  of  Attitude  Limits 


There  are  certain  .practical  limitations  on  the  freedom  of  choice 
for  the  upper  and  lower  limits  for  both  the  angle  of  attack  and  the  bank  angle. 
Most  of  the  limitations  are  due  to  the  problems  associated  with  obtaining  as 
small  a  value  as  possible  for  the  peak  temperature  experienced  during  the  reentry 
of  the  spacecraft.  In  addition,  the  allowable  range  of  values  for  the  upper 
and  lower  limits  changes  with  each  phase  of  the  reentry.  Consequently,  the  best 
procedure  will  be  to  take  each  phase  in  turn  and  describe  the  reasons  for  the 
restrictions  during  that  part  of  the  trajectory. 

(a)  During  Phase  I,  there  are  no  restrictions  and  no  practicaT“limita- 
tions  since  the  atmosphere  is  not  effective  in  controlling  the  trajectory.  The 
only  point  at  which  the  attitude  must  be  controlled  is  at  the  beginning  of  Phase 
1  during  the  retro-fire.  At  that  time,  the  attitude,  control  must  be  quite  pre¬ 
cise,  but  there  is,  of  course,  no  restriction  on  the  choice  of  the  angles.  Near 
the  end  of  Phase  I,  however,  the  angle  of  attack  should  be  somewhere  near  the 
value  to  be  used  during  Phase  II  since,  otherwise,  stability  and  control  problems 
might  be  encountered. 

(b)  During  Fnase  II,  the  bank  angle  should  always  be  zero  in  order  to 
produce  the  maximum  possible  vertical  force  for  a  given  angle  of  attack.  It 
might  be  feasible  to  introduce  some  aerodynamic  siue  force  on  the  vehicle  at  this 
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time,  but  the  results  would  hardly  be  worth  the  added  complexity  necessary  in 
the  guidance  system.  This  is  because  only  large  angles  of  roll  would  really  be 
effective  in  producing  cross  range  at  this  high  a  speed  (Ref.  6),  but  the 
associated  reduction  in  vertical  force  could  not  be  tolerated. 

The  angle  of  attack  should  always  be  near  the  optimum  value  of 
55  degrees  which  gives  the  maximum  lift  coefficient.  In  the  presence  of  a  pre¬ 
determined  atmosphere  in  which  there  were  no  density  variations  about  the  mean 
curve  of  density  versus  altitude,  it  would  be  possible  to  determine  the  desired 
angle  of  attack  which  would  then  be  held  throughout  Phase  II.  However,  in 
general,  the  guidance  system  will  probably  be  making  small  changes  throughout 
Phase  II  in  order  to  compensate  for  predicted  errors  in  the  terminal  down  range 
position.  Since  the  CONTAC  system  does  not  go  into  active  operation  until  the 
end  of  Phase  II,  the  angle  of  attack  is  completely  under  control  of  the  guidance 
system  during  this  time. 

The  largest  angle  of  attack  which  was  used  herein  during  Fhasell 
was  65  degrees.  However,  values  as  high  as  75  degrees  would  probably  be  feasible 
if  the  spacecraft  could  withstand  the  resulting  higher  heating  rates  at  the 
rear  of  the  undersurface  of  the  vehicle.  The  smallest  value  used  was  35  degrees. 
When  still  smaller  values  were  tried,  the  first  maximum  temperature  was  higher 
than  the  peak  temperature  for  an  equilibrium  glide. 

The  designation  of  the  upper  limit  on  the  attitude  (no  lower  limit 
was  used)  was  the  entry  angle  of  attack  (which  varied  between  35  degrees  and  65 
degrees)  and  the  entry  angle  of  roll  (which  was  always  zero) . 

(c)  During  Phase  III,  the  CONTAC  system  is  in  active  operation  about 

half  the  time.  The  vehicle  is  initially  rolled  to  one  side  or  the  other  and  the 
nose  pitched  down  in  order  to  prevent  the  temperature  from  falling. 

The  lower  limits  are  designated  as  the  "skip  angle  of  attack" 
and  the  "skip  bank  angle".  The  reason  for  this  terminology  is  that  the  space¬ 
craft  usually  performs  a  nearly  constant  altitude  glide  at  the  beginning  of 
Phase  III,  involving  a  very  small  skip,  or  increase  in  altitude,  of  between  zero 
and  20,000  feet.  For  some  situations,  especially  when  the  entry  angle  of  attack 
is  large,  the  skipdoes  not  occur;  however,  the  designation  is  still  retained. 

Note  that  the  skip  angle  of  attack  and  the  skip  bank  angle  are  con¬ 
trolled  by  the  range  guidance  system.  These  attitude  angles  are  always  the 
actual  values  held  by  the  vehicle  whenever  the  temperature  is  falling.  The  skip 
angles  may,  of  course,  be  changed  at  any  time  during  Phase  III,  but  an  immediate 
change  in  predicted  down  range  will  not  result  if  the  skip  angles  are  not  in 
current  use  due  to  override  by  the  CONTAC  system. 

The  largest  value  used  for  the  skip  angle  of  attack  was  55  degrees; 
the  smallest  value  was  5  degrees.  The  smallest  value  might  require  adjustment  if 
its  use  results  in  upper  surface  heating  which  cannot  be  tolerated. 

The  "largest"  value  used  for  the  skip  angle  of  roll  was  zero  and 
the  "smallest"  value  w.  -  -90  degrees.  Notice  that  the  use  of  a  bank  angle  to  the 
left  preserves  the  algebraic  definition  of  largest  and  smallest.  If  roll  angles 
to  the  right  are  employed,  the  largest  would  still  be  zero  and  the  "smallest" 
would  then  become  +90  degrees.  The  key  point  to  keep  in  mind  is  that  the  largest 
roll  angle  must  always  occur  for  the  maximum  vertical  force  (which  obviously  forces 
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the  largest  roll  angle  to  be  zero)  and  the  "smallest"  bank  angle  for  the  minimum 
vertical  force. 


The  upper  limits  during  Fhase  III  are  simply  designated  the  "upper 
limits" .  The  values  are  always  the  ones  which  give  the  maximum  vertical  force 
which  is,  of  course,  obtained  at  the  angle  for  the  maximum  lift  coefficient. 
Therefore,  the  upper  limit  for  the  single  of  attack  is  55  degrees  and  the  upper 
limit  for  the  bank  angle  is  zero. 

(d)  Phase  IV  always  starts  with  the  attitude  at  the  upper  limits,  the 
same  as  the  upper  limits  for  Phase  III. 

The  lower  limits  for  Phase  IV  are  designated  "lower  limits". 

They  are  used  when  the  speed  falls  appreciably  and  the  CONTAC  system  goes  into 
passive  operation.  Within  the  limitations  imposed  by  the  stability  and  avail¬ 
able  aerodynamic  control  of  the  vehicle,  the  guidance  system  is  free  to  choose 
any  values  which  it  finds  necessary  to  control  the  range.  In  general,  the  lower 
limit  for  the  angle  of  attack  will  lie  in  the  vicinity  of  the  angle  for  (l»/D)max 
while  the  lower  limit  for  the  bank  angle  is  varied  in  the  way  which  will  result 
in  the  maximum  cross  range  (if  necessary) . 

(e)  During  Phase  V,  the  CONTAC  system  is  no  longer  in  operation.  All 
values  for  the  attitude  are  dictated  by  guidance  and  control  problems. 

2.3  CONTAC  Program  Flow  Diagram  (see  Fig.  12) 

The  basic  purpose  of  the  CONTAC  system  is  to  prevent  a  change  in 
the  temperature.  The  digital  implementation  of  the  basic  concept  is  quite  simple. 
The  program  to  simulate  and  test  the  concept  is,  of  course,  much  more  detailed 
only  because  of  other  requirements  such  as  program  error  checks  and  printouts  of 
intermediate  calculations.  In  general,  conflicts  often  occur  and  the  extra 
effort  to  take  account  of  the  numerous  possibilities  is  considerable. 

In  the  absence  of  these  other  details,  the  CONTAC  system..logic  is 
shown  in  Fig.  12.  There  must  be  some  memory  (but  no  predictive  requirement)  so 
that  the  system  can  remember  whether  it  is  testing  for  a  local  minimum  or  maximum, 
and  for  the  storage  of  Tlastmin  or  Tiast-g^,  respectively.  If  the  latter  is 
currently  being  experienced,  then  when  the  temperature  has  fallen  below  the 
current  local  maximum  by  a  reasonable  margin  (TfleadbandK  corrective  action  is 
taken.  The  determination  of  a  local  minimum  is  similar. 

However  the  simulation  of  the  CONTAC  system  by  a  digital  computer 
is  prone  to  difficulties  in  logic  flow  and  numerical  calculation  which  have  no 
counterpart  in  the  real  situation.  The  basic  cause  is  the  method  by  which  a 
digital  computer  must  perform  the  calculations.  As  opposed  to  a  continuous  value 
for  all  variables,  as  is  found  in  the  output  of  an  analogue  computer,  or 
analogue  instruments,  a  digital  computer  calculates  the  flight  path  only  at 
specific  points  along  that  path.  As  a  consequence,  the  deadband  test  shown  in 
Fig.  12  may  not- be  satisfied  at  one  point  in  the  calculation,  but  will  be 
satisfied  by  too  great  a  margin  at  the  next  point.  If  the  digital  calculation 
is  proceeding  in  steps  of  two  or  three  seconds,  the  time  lag  introduced  by  the 
digital  calculation  is  totally  unacceptable  (Fig.  13)  and  measures  must  be  found 
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to  correct  the  situation*.  The  simplest  method  would  be  to  reduce  the  time 
interval  (actually,  the  flight  path  distance)  between  the  points  in  the  calcu¬ 
lation.  However,  what  value  should  be  used?  A  step  size  small  enough  to 
satisfy  the  deadband  test  requirement  produces  an  unacceptably  large  increase 
in  computing  time  on  the  machine.  A  step  size  which  increases  computer  time 
by  only  a  reasonable  amount  will  not  always  be  able  to  provide  an  acceptably 
short  time  delay. 

The  solution  adopted  was  to  provide  the  calculation-  section  with 
backup  capability.  As  shown  in  Fig.  14,  new  points  are  calculated  at  shorter 
and  shorter  intervals  until  both  questions  (Fig.  lg)  are  satisfied.  At  that 
point,  the  distance  travelled  and  the  corresponding  time  delay  between  the 
actual  position  and  the  position  at  backup  point  No.  5  (Fig.  14)  is  acceptable 
and  the  commands  to  adjust  the  angle  of  attack  and  the  bank  angle  are  calculated 
and  implemented. 

The  backup  capability  allows  the  calculation  section  to  proceed 
with  step  sizes  which  Eire  as  lfiurge  as  calculation  errors  permit  and  thereby 
perform  the  calculation  as  fast  as  possible.  The  choice  of  the  smallest  step 
size  for  the  test  is  based  upon  practical  experience  with  the  runs  and  usually 
results  in  a  maximum  time  delay  of  between  ten  Eind  twenty  milliseconds.  During 
test  runs  to  determine  the  maximum  acceptable  value,  it  was  found  that  delays 
which  were  eui  order  of  magnitude  greater  produced  results  which  were  only 
negligibly  different.  However,  for  still  larger  time  delays,  a  discrepancy  was 
noted  sometimes. 

2.4  Vehicle  Implementation 

The  problems  of  isolating  the  control  logic  from  the  details 
of  the  calculation  will,  of  course,  disappear  once  the  system  is  placed  on 
board  the  spacecrEift  and  actually  put  into  operation.  As  shown  in  Fig.  15,  the 
calculation  section  cf  Tig.  12  will  be  replaced  by  the  vehicle  dynamics.  The 
vehicle  sensors  will  provide  the  forwEird  interface  to  the  control  hardware  and 
the  reaction  jets  will  be  the  rearward  interface.  In  the  case  of  a  digital 
control  system,  the  sampling  rates  for  the  einalogue  to  digital  convertor 
will  probably  be  on  the  order  of  ten  milliseconds  or  less  for  the  temperature. 
Other  variables  need  to  be  sensed  less  often  and  perhaps  only  when  needed. 

The  absence  of  inner  return  loops  in  the  vehicle  control  logic 
merely  indicates  that  no  signal  is  necessary  until  the  appropriate  test  is 
satisfied.  However,  since  the  calculation  will  be  in  real  time,  a  great  deal  of 
attention  must  be  given  to  the  way  in  which  the  available  computational  time  is 
divided.  Because  the  CONTAC  system  is  often  in  the  passive  mode,  i.e.  no  active 
attitude  control  is  being  performed,  time  shEiring  procedures  will  probably 
msike  the  bestuuse  of  the  available  computational  capacity.  The  range  controller 
will  undoubtably  use  the  major  share  of  the  computer's  time  for  predictive  cal¬ 
culations  since  in  view  of  the  complex  corrections  required  by  the  disturbing 
effects  of  the  EEirth's  rotation  and  oblateness  it  is  doubtful  whether  any 
simple  predictive  scheme  would  be  satisfactory  for  the  size  of  footprint  within 
which  the  CONTAC  system  is  able  to  operate.  However,  a  simple  interrupt  system 


*  It  is  not  so  much  the  length  of  the  time  lag  that  complicates  the  situation, 
but  the  fact  that  it  is  unpredictable  and  varies  from  zero  to  perhaps  five  or 
even  ten  seconds.  Since  no  real  system  will  possess  such  a  random  variation  in 
the  size  of  the  deadband,  results  computed  on  this  basis  would  be  highly  suspect. 


will  allow  for  the  monitoring  of  the  temperature  on  a  scheduled  basis  and  the 
subsequent  calculation  of  the  attitude  control  moments  when  required. 

The  control  scheme  proposed  here  requires  the  measurement  or 
computation  of  angles  of  attack  and  roll  and  their  derivatives.  No  specific 

means  are  recommended  for  obtaining  this  information.  However,  one  possibility 

is  to  compute  them  from  the  outputs  of  an  intertial  platform,  which  can  provide 
both  velocity  and  attitude  in  inertial  space.  Other  techniques  are  also  possible 
within  the  state  of  the  art.  When  the  CONTAC  system  commands  a  change  of  a  or 
<p,  it  is  assumed  that  pitch  or  roll  jets  are  turned  on,  producing  a  constant 

moment  with  zero  time  lag.  Both  |a|  and  |<p|  are  constrained  to  be  less  than 

certain  selected  maxima,  and  when  these  are  reached,  the  control  jets  are  turned 
off. 

An  extensive  effort  to  produce  a  thermocouple  which  is  capable  of 
measuring  the  temperature  at  the  nose  stagnation  point  of  the  vehicle  is  des¬ 
cribed  in  Ref.  4.  Actual  tests  were  conducted  which  showed  that  such  a  device 
was  able  to  measure  values  in  excess  of  3500°F  for  several  minutes.  Since  the 
peak  temperature  which  occurred  for  the  vehicle  used  in  this  investigation  was 
usually  below  3000*%,  there  will  not  be  any  difficulty  in  procuring  a  satisfactory 
measuring  instrument. 

Another  basic  problem  is  to  determine  the  actual  stagnation  point 
temperature.  The  position  of  the  thermocouple  is,  of  necessity,  fixed  ih‘  the 
nose  of  the  vehicle.  When  changes  are  made  in  the  angle  of  attack,  the  thermo¬ 
couple  will  no  longer  be  at  the  stagnation  point  and  some  method  must  be  found 
to  determine  the  actual  stagnation  point  temperature. 

The  solution  to  this  problem  is  also  given  in  Ref.  4.  The  re¬ 
lationship  between  the  measured  and  the  stagnation  point  temperatures  is  given  by 

^measured  =  ^stagnation  cos^^  (‘-'^stagnation) 

At  high  Mach  numbers  and  for  positions  up  to  4o  degrees  off  the  stagnation  point, 
this  relationship  is  in  excellent  agreement  with  both  theoretical  and  experi¬ 
mental  data  for  both  hemispheres  and  swept  cylinders.  When  the  CONTAC  system 
is  active,  the  maximum  variation  in  the  angle  of  attack  was  from  5  degrees  to 
55  degrees.  The  thermocouple  position  could  be  at  the  stagnation  point  when. the 
angle  of  attack  of  the  vehicle  is  at  30  degrees.  Consequently,  the  maximum 
offset  would  be  only  25  degrees.  For  a  stagnation  point  temperature  of  3000°F, 
the  temperature  measured  when  the  angle  of  attack  of  the  vehicle  is,  55  degrees 
is  only  4.2  percent  less  or  2875°F.  The  derivative  of  the  temperature  with 
respect  to  the  angle  of  attack  is  al-<=o  very  small.  The  corresponding  value 
would  be  9*7  F°  /  degree  change  in  angle  of  attack. 

In  general,  the  vehicle  will  experience  a  small,  but  variable 
degree  of  side  slip,  owing  to  the  presence  of  inadvertent  yawing  movements  from 
asymmetries  and  control  action.  Consequently,  two  extra  thermocouples  could 
be  provided  to  measure  the  sideslip  angle  and  to  control  reaction  jets  which 
would  achieve  a  null  angle.  The  suggested  positions  are  on  the  wing  leading 
edges;  offset  nose  positions  would  also  be  satisfactory. 

Note  that  regardless  of  the  position  used  for  the  thermocouples, 
several  extra  sensors  would  be  used  in  parallel  at  each  location  in  order  to 
guarantee  reliability. 


2.5  Vehicle  Attitude  Control 


As  part  of  any  simulation  of  a  complex  physical  situation,  a 
decision  must  be  made  as  to  just  how  accurately  the  mathematical  model  should 
be  abie  to  duplicate  the  real  situation.  In  the  case  of  the  attitude  control 
system,,  information  is  available  (Refs.  10  and  ll)  which  makes  possible  a  con¬ 
siderable  simplification  in  the  complexity  of  the  mathematical  model.  These 
references  indicate  that  the  short  period  mode  possesses  time  constants  orders 
of  magnitude  greater  than  the  characteristic  times  of  the  CONTAC  system.  As 
a  result,  most  of  the  terms  in  the  general  equations  describing  the  attitude 
of  the  vehicle  may  be  eliminated  -  in  particular,  both  the  damping  and  spring 
constants  are  deemed  negligible.  Consequently,  only  the  inertial  term  remains 
in  the  equations  for  rotation  about  any  axis.  It  is  also  assumed  that  both  the 
Earth' s  rotation  and  the  movement  of  the  vehicle  within  the  Earth's  reference 
system  will  have  a  negligible  effect  on  the  attitude  equations. 

•In  keeping  with  the  above  assumptions,  the  equations  of  motion 
for  the  attitude  become  very  simple.  Only  the  equation  for  the  angle  of  attack 
will  be  developed  since  the  bank  angle  equation  is  similar.  The  basic  equation  is 

Ia  a  =  M  (2.1) 

where  M  is  the  constant  pitching  moment  supplied  by  the  pitch  control  jets 
when  they  are  turned  on.  Thus 

O  -  kd  (2.2) 

where  is  zero  when  the  control  jets  are  off,  positive  when  set  for  nose-up 
change,  and  negative  for  nose-down  change. 


Assuming  that  the  angle  of  attack  was  decreasing  when  the  pitch  jets  were  first 
turned  on  and  that  the  angle  of  attack  should  now  begin  to  increase  (i.e.  is 
positive),  integration  gives 


<4 


(t-ts)<A  t 
(t-t  )>At 


(2.3a) 

(2.3b) 


where 

At  =  (dmax  -  &s)/*cc 


(2.4) 


A  final  integration  yields 

_  fas  +  M^s)  +  ^(t-t s)2  (t-ts)<At 

as  +  Aa  +  dmax  j^(t-ts)-AtJ  (t-ts)>At 

where 

=  (^max  +  ^s^max^s^ka 


(2.5a) 

(2.5b) 

(2.6) 


Figure  16  shows  the  various  relationships.  Since  equations 
(2.5a)  and  (2.5b)  are  simple  analytical  expressions  which  ere  a  function  only 
of  time,  they  are  easily  introduced  into  the  equations  of  motion  presented  in 
Section  V.  One  point  which  may  require  explanation  is  the  reason  for  the  limit 
placed  on  d.  On  the  one  hand,  as  large  a  value  as  possible  is  desirable  in 
order  to  be  able  to  respond  as  quickly  as  possible  to  rapid  changes  in  atmos- 


pheric  density  that  require  a  rapid  adjustment  of  the  angle  of  attack.  On  the 
other  hand,  if  the  inertia  of  the  vehicle  (with  its  negative  angular  velocity  - 
remember  that  as  was  negative)  is  to  be  quickly  overcome,  Eq.  (2.4)  shows  that 
a  large  &max  requires  a  large  angular  acceleration.  In  addition,  high  angular 
rates  lead  to  frequent  cycling  of  the  control  which  requires  a  high  fuel  ex¬ 
penditure  for  the  reaction  jets.  Consequently,  a  reasonable  compromise  must 
be  found  between  &max,  the  available  angular  acceleration,  and  the  available 
mass  of  fuel.  The  values  used  in  this  investigation  were  quite  conservative. 
Appendix  D  gives  a  sample  set  of  values  prior  to  the  actual  output. 

Slight  variations  in  sign  in  Eq'S'.  (2.3)  to  (2.6)  are  required 
when  is  negative.  Also,  two  sets  of  bank  angle  equations  are  required,  one 
for  positive  k<p  and  one  for  negative  kcp  .  The  bank  angle  equations  are  a  little 
more  complicated  to  use  since  the  optional  right  or  left  hand  turn:,  requires 
additional  logic  within  the  digital  computer  program.  For  a  left  hand  turn, 
all  the  a  symbols  need  only  be  replaced  by  <p  to  convert  equations  (2.3)  to  (2.6) 
into  bank  angle  equations  for  a  positive  k„  .  Of  course,  the  basic  equation 
for  the  bank  angle  .  v 

(p  =  kqp  (2.7) 

which  is  exactly  the  parallel  of  Eq.  (2.2)  and  applies  equally  to  all  situations. 


III.  VEHICLE  CHARACTERISTICS 
3-1  Introduction 


In  recent  years,  the  general  trend  in  lifting  vehicle  programs 
has  been  a  shift  away  from  configurations  with  relatively  high  lift  to  drag 
ratio  (L/D  >  2.5),  such  as  the  X-20  Dynasoar,  to  a  lifting  body  style  as  typified 
by  the  M2-F2,  the  HL-10,  and  the  SV-5  (Ref. 12)  wingless  reentry  vehicles.  The 
reasons  for  the  choice  seem  to  have  been  based  primarily  upon  the  advantages 
associated  with  much  lower  temperatures  resulting  from  the  elimination  of  small 
nose  cap  and  leading  edge  radii.  The  disadvantages  caused  by  a  much  smaller 
footprint  or  landing  area  capability  are  not  severe  enoqgi  to  cause  serious 
concern.  By  reason  of  the  Earth's  rotation  under  what  is  essentially  a  fixed 
inertial  orbit,  these  wingless  vehicles  are  still  able  to  land  at  any  point  on 
the  Earth  which  falls  underneath  or  at  least  near  the  orbital  envelope  (Fig.  5) 
without  having  to  wait  more  than  24  hours.  Naturally,  the  coverage  is  severely 
restricted  for  equatorial  or  near-equatorial  orbits;  but  even  the  Dynasoar 
could  not  claim  complete  coverage  for  equatorial  orbits  without  incurring 
extremely  high  temperatures. 

3.2  Vehicle  Geometry 

The  vehicle  chosen  for  this  study  is  quite  similar  to  current 
design  concepts.  Its  mass  and  wing  area  never  individually  enter  the  equations 
of  motion;  consequently,  they  have  been  left  in  the  form  of  the  wing  loading 
ratio.  This  parameter  was  then  chosen  from  about  the  middle  of  the  range  which 
is  presently  under  review.  The  value  of  40  lb/ft^  is,  therefore,  simply  re¬ 
presentative  . 


The  actual  vehicle  configuration  is  a  delta  wing  planform  with 
a  leading  edge  sweep  of  about  70  degrees.  Although  the  mass  and  wing  area  are 
not  needed  for  the  equations  of  motion,  nevertheless  some  estimate  of  the  nose 
cap  and  leading  edge  radii  is  desirable.  A  mass  of  about  8000  lb.  gives  a  wing 
area  of  about  200  sq.  ft.  which  would  then  lead  to  abase  length  of  about  17  ft. 
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However,  the  body  would  be  only  about  20  ft.  long  since  around  5  ft.  of  the  nose 
would  probably  be  cut  off  to  provide  a  nose  cap  radius  in  the  vicinity  of  2.5  ft. 
The  leading  edge  radius  would  also  lie  in  this  region. 

3.3  Vehicle  Aerodynamics 

One  of  the  most  vital  requirements  for  any  investigation  which 
seeks  to  predict  the  reentry  flight  path  of  a  spacecraft  with  a  reasonable 
degree  of  accuracy  is  a  proper  representation  of  aerodynamic  forces  throughout 
all  phases  of  the  trajectory.  However,  since  the  primary  object  of  this  study 
is  to  determine  the  effectiveness  of  the  COHTAC  system  under  varying  operating 
conditions,  only  that  portion  of  the  flight  where  the  heating  rates  are  extreme 
will  require  highly  accurate  models  for  the  aerodynamic  forces  which  are  exerted 
on  the  vehicle.  Portions  of  the  trajectory  that  are  relatively  unimportant  to 
the  primary  purpose  of  the  analysis  may  employ  a  model  which  is  able  to  give 
only  a  reasonable  approximation  for  the  lift  and  drag  coefficients. 

Such  is  the  case  with  the  aerodynamic  model  chosen  to  represent 
the  lift  and  drag  forces  on  the  vehicle  under  consideration.  Several  experi¬ 
mental  studies  (Refs.  13  to  15)  show  that  expressions  based  on  Newtonian  theory 
adequately  represent  the  sum  total  of  the  aerodynamic  forces  which  are  developed 
over  a  spacecraft  during  its  reentry  into  the  Earth’s  atmosphere  as  long  as  the 
Mach  number  is  greater  than  about  six  or  eight.  As  a  result,  a  number  of  authors 
(Refs.  3>  l6,  and  17)  have  used  either  an  exact  or  modified  form  ox-  the  flat 
plate,  hypersonic,  Newtonian  approximation.  Based  on  Ref.  l4,  the  best  esti¬ 
mate  for  the  pressure  coefficient  on  a  flat  plate  is 

Cp  =  1.75  sin2a  (3.1) 

The  value  assigned  to  the  minimum  drag  coefficient,  Cpmin,  is 
determined  by  the  vehicle  configuration,  primarily  by  the  degree  of  nose  cap 
and  leading  edge  bluntness.  The  nose  cap  radius  should  be  as  large  as  possible 
in  order  to  reduce  heating  (See  Sec.  3«^)>  however,  a  large  value  will  produce 
a  maximum  lift-drag  ratio  which  is  too  small  to  permit  the  required  manoeuvra¬ 
bility.  The  compromise  value  of  nose  radius  given  in  Sec.  3.2  yields  a  value 
for  the  minimum  drag  coefficient  of 

CDmin  =  0-0^25  (3-2) 

If  the  value  of  the  nose  cap  radius  which  has  been  suggested  does  not  produce 
the  above  value  for  Cp^  ,  an  adjustment  to  the  former  may  easily  be  made  without 
unduly  affecting  the  peak  temperature,  since  the  latter  varies  only  as  the  eighth 
root  of  the  nose  radius  (see  Eq.  3*7). 

The  use  of  the  preceding  values  for  Cp  and  Cp  .  gives  the  modi¬ 
fied  Newtonian  expressions  for  the  T*.xt  and  drag  coefficients,which  are 

CL  =  1.75  sin^a  cosa  (3.3) 

CD  =  0.0625  +  1.69  sin^a  (3.4) 

These  expressions  are  valid  for  0<O<90  degrees  and  are  quite  accurate  down  to 
Mach  numbers  of  six  or  seven.  For  Mach  numbers  smaller  than  five,  there  is 
every  indication  that  Eqs.  (3*3)  and  (3*4)  continue  to  be  quite  representative 
even  down  through  the  transonic  region.  References  18  to  20  indicate  that  for 


angles  of  attack  below  about  twenty  or  thirty  degrees,  the  lift  and  drag 
coefficients  continue  to  exhibit  the  same  general  pattern  while  the  L/D  slowly 
increases  with  decreasing  Mach  number.  Consequently,  Eqs.  (3-3)  and  (3*4)  may 
continue  in  use  down  through  the  high  subsonic  range  if  conservative  estimates 
of  manoeuvring  capability  are  acceptable. 

3.4  Equilibrium  Nose  Stagnation  Point  Temperature 

Since  the  primary  goal  of  this  investigation  is  to  determine  the 
applicability  of  the  COOTaC  system,  the  exact  numerical  magnitude  of  the  nose 
stagnation  point  convective  heating  rate  is  relatively  unimportant.  Indeed,  it 
is  even  questionable  whether  extremals  in  the  stagnation  point  temperature  (see 
Pig.  14  as  to  how  a  maximum  or  minimum  T  is  handled)  which  occur  as  a  result  of 
reversals’  in  the  gradient  of  the  heating  rate  must  be  accurately  predicted  by 
the  model  which  has  been  adopted.  On  the  other-  hand,  the  heating  rate  must  be 
simulated  in  a  way  which  is  sensitive  to  changes  in  atmospheric  and  flight  con¬ 
ditions,  especially  at  the  times  when  the  model  predicts  that  the  maximum  rate 
is  being  experienced.  In  other  words,  it  is  not  at  all  important  that  the  pre¬ 
dicted  heating  rate  reach  a  maximum  or  minimum  at  precisely  the  same  time  that 
the  actual  heating  rate  reaches  its  extremal  values  (in  fact,  it  is  very  dubious 
whether  any  model  which  simulates  the  heating  rate,  no  matter  how  complex  or 
sophisticated  it  might  be,  could  predict  the  exact  time  at  which  the  heating 
rate  derivative  changes  sign);  it  is  only  necessary  that  if  the  heating  rate 
is  to  remain  relatively  constant  as  a  result  of  attitude  control  for  the  simu¬ 
lated  model  that  this  condition  also  occur  in  the  real  situation  and  at  about 
the  same  time  in  the  flight.  This  rather  lax  requirement  is  due  to  the  adaptive 
nature  of  the  CONTAC  system  which  must  only  prevent  any  large  changes  in  the 
heating  rate  as  long  as  possible,  i.e.  as  long  as  the  attitude  angle  has  not 
reached  its  upper  or  lower  limit. 

The  particular  form  of  the  nose  stagnation  point  convective  heat¬ 
ing  rate  which  is  used  in  this  investigation  has  also  been  adopted  by  several 
other  authors  (Refs.  3,  17,  21,  and  22)  as  being  sufficiently  representative  of 
experimental  results.  The  basic  relationship  is 

n  =  17,600  (  V  P 

^  *JTn  \2b,000 )  \0. 002378 

Since  the  value  which  is  to  be  measured  on  board  the  vehicle  is  the  temperature, 
this  heating  rate  is  then  inserted  into  the  temperature  rate  equation  which  is 

-5£-  §  .q-EcrS1'  (3-6) 

s  dt 

The  radiation-equilibrium  temperature  is  defined  by  setting  the  left  hand  side 
of  Eq.  (3.6)  equal  to  zero.  When  an  actual  measurement  i£  being  made,  the  equili¬ 
brium  value  is  produced  by  using  such  a  small  value  of  m  c/S  at  the  temperature 
sensor  that  the  left  hand  side  becomes  negligible  and  an  excellent  approximation 
to  the  equilibrium  value  is  achieved  very  quickly.  The  data  in  Ref.  4  indicates 
that  there  will  be  a  relatively  negligible  time  lag  in  reaching  the  equilibrium 
temperature  since  the  time  constants  associated  with  the  cyclic  commands  pro¬ 
duced  by  the  CONTAC  system  are  of  the  order  of  five  to  ten  seconds. 

Then,  setting  to  zero  the  left  hand  side  of  Eq.  (3*6)  and  sub¬ 
stituting  Eq.  (3-5)  for  the  heating  rate,  qn,  gives  the  radiation-equilibrium  nose 


J'5 


(3.5) 
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stagnation  point  temperature  as 


/ 


The  only  terms  which  must  have  values  assigned  are  €  and  Rn.  Normally,  e  is 
between  0.5  and  C.8  while  Rn  is  about  2.5  feet  for  the  vehicle  under  considera¬ 
tion.  Thus,  the  products  of  &fRn  would  range  from  0.8  to  1.25  due  to  the  un¬ 
certainty  in  e  .  As  a  consequence,  a  "standard"  value  for  the  temperature  is 
used  for  which  a  value  of  unity  is  takerrfar  £\(5n.The  advantage  is  a  standard 
from  which  the  proper  result  for  the  temperature  may  be  computed  when  the  correct 
values  of  €  and  Rn  are  known  without  the  disadvantage  of  being  greatly  in  error. 
This  standard  stagnation  point  temperature  is  used  in  all  the  calculations  as  a 
measure  of  the  actual  temperature  which  would  result  for  a  device  (e.g.  a  sophis¬ 
ticated  thermocouple)-  placed  at  the  nose  stagnation  point  of  the  spacecraft. 

This  variable  will  henceforth  be  referred  to  simply  as  the  "temperature"  (as 
already  noted  in  Section  1.2).  Then,  in  view  of  the  above  assumptions,  the 
temperature  is  given  by 


T4  =  Q  V3*15  »Jp 


(3.8) 


where 


Q  =  17,600/cj26,0003*15  */0. 002378  (3.9) 


or 


Q  =  95397 


(3.10) 


IV.  GEOPHYSICAL  ENVIRONMENT 


4.1  Introduction 


One  of  the  most  important  goals  of  this  investigation  is  to 
study  the  effect  of  non- ideal  conditions  upon  the  peak  temperature,  primarily 
with  a  view  to  determining  the  extent  to  which  these  conditions  impair  the 
performance  of  the  CONTAC  system. 

There  are  two  distinct  types  of  non- ideal  conditions  which  are 
of  concern.  The  first  includes  all  the  effects  which  may  be  predicted  well 
in  advance  of  the  reentry.  These  effects  include  the  Earth's  rotation,  the 
Earth's  oblateness  (i.e.  the  distortions  of  both  shape  and  gravitational  po¬ 
tential),  and  the  geographical  variations  of  the  time  average  density  (which 
are  a  function  of  both  latitude  and  longitude  as  well  as  altitude).  The 
second  type  includes  those  effects  which  are  the  result  of  departures  from  the 
time  averaged  mean  values  of  the  density. 

A  third  source  of  disturbance  is  atmospheric  wind.  No  data  could 
be  found  showing  the  spatial  gradient  of  the  wind  at  the  altitudes  of  concern 
of  this  investigation  (the  same  problem  was  encountered  with  density  variations 
and  will  be  discussed  in  Sec.  4.6).  However,  local  wind  velocities  of  up  to 
500  feet  per  second  do  exist.  Actually,  almost  all  values  are  lower  than  250 
feet  per  second,  but  on  isolated  occasions  they  approach  the  higher  figure. 
Assuming  a  step  increase  of  500  feet  per  second  (from  -500  to  +500  feet  per 
second  seems  far  too  large  a  difference  to  support  even  if  both  conditions  did 
occur  within  a  reasonable  proximity) ,  the  percentage  increase  in  the  velocity 
at  a  speed  of  20,000  feet  per  second  would  amount  to  only  2g-  percent.  The 
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resulting  increase  in  temperature  would  be  slightly  less  (see  Eq.  3*8).  The 
percentage  increase  would  be  greater  at  lower  velocities,  but  by  the  time  the 
speed  falls  below  about  16,000  feet  per  second,  the  temperature  is  no  longer 
a  problem.  Since  it  would  be  the  changes  in  wind  speed  that  cause  the 
difficulties  (see  Sec.  8.4),  it  is  clear  that  w'.nds  should  not  present  a  problem 
for  the  COKTAC  system. 


4.2  Rotation  of  the  Earth 


There  is  a  relative  pseudowind  present  due  to  the  Earth's  rotation. 
The  West  to  East  speed  at  the  equator  amounts  to  about  1500  feet  per  second 
relative  to  inertial  space.  For  the  calculation  of  aerodynamic  heating  and 
forces,  this  relative  velocity  may  be  subtracted  from  the  initial  entry  speed 
for  equatorial  West  to  East  reentries  and  provides  e.  small  but  not  insigni¬ 
ficant  reduction  in  the  peak  temperature  experienced  during  the  flight. 

The  actual  rotational  motion  of  the  Earth  is  quite  complex. 

Parvin  (Reference  23)  has  done  an  excellent  summary  of  the  various  cotions 
and  lists  five  separate  parts  to  the  rotation  of  the  Earth. 

The  primary  motion  is,  of  course,  the  familiar  rotation  of  the 
Earth  about  its  axis  which  produces  the  normal  sunrise  and  sunset.  This  move¬ 
ment  takes  by  definition,  exactly  one  solar  day  to  complete.  However,  the 
rotational  speed  of  interest  is  with  respect  to  inertial  space  and  for  this 
figure  the  sideral  day  may  be  used.  The  length  of  time  taken  for  one  revolu¬ 
tion  in  inertial  space  is  23  hours,  56  minutes,  and  4.09  seconds.  This  works 
out  to  a  rotational  speed  of  0.004l7807  degrees  of  longitude  per  seconu. 

The  other  four  parts  of  the  rotation  will  be  mentioned  for 
interest,  but  are  all  negligible  with  regard  to  the  effect  they  have  on  the 
reentry  of  a  spacecraft  First,  there  is  the  precessional  rotation  about  the 
23-5  degree  inclination  to  the  ecliptic  pole  which  takes  about  25,800  years. 
Second,  other  bodies  in  the  solar  system  tend  to  cause  a  slight  irregularity 
in  the  precessional  rate  called  nutation;  one  component  with  less  than  ten 
seconds  of  arc  has  a  period  of  nineteen  years  while  the  other  component  with 
less  than  one-half  second  of  arc  has  a  period  of  only  two  weeks.  lfext,  the 
pole  itself  wanders,  which  would  tend  to  affect  the  alignment  cf  the  Earth’s 
axis;  however,  although  the  two  periods  are  relatively  short  -  twelve  and 
fourteen  months  -  the  angles  are  less  than  four  seconds  of  arc.  Finally, 
there  is  an  additional  revolution  due  to  the  spiral  motion  of  the  galaxy  which 
occurs  once  approximately  every  200  million  years.  A  better  and  much  more  com¬ 
plete  description  of  these  four  additional  motions  may  be  found  in  the  reference. 
It  is  quite  clear,  however,  that  they  are  too  small  to  have  any  effect  on  the 
reentry. 


4.3  Earth's  Oblate  Gravitational  Field 


During  the  past  ten  years,  the  use  of  artificial  satellites  to 
measure  the  potential  field  of  the  Earth  has  shown  the  inadequacy  of  represent¬ 
ing  the  gravitational  attraction  as  a  simple  central  force.  Using  the  notation 
of  this  paper,  the  general  representation  of  the  Earth's  potential  field  as 
given  by  Kozai  (Ref.  24)  is 
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rag 


are  spherical  hormonics.  Several  other  formulations  exist,  but  all  are  essentially 
the  same  in  that  they  use  an  infinite  series  of  spherical  harmonics  to  represent 
the  potential  field. 

The  gravitational  vector  is  divided  into  its  three  principal 
components  which  results  in  (see  Sec.  5*2  for  a  definition  of  the  fn,e,v} 
reference  system) 

g  =  n  +  e  gQ  -  v  gr  (4.3) 

and  requires  the  expressions  for  each  of  the  gravitational  components  which  are 
defined  by  the  spherical  gradient  of  the  potential  as 


dU 

gr  <5r 


1  dU 
r  3X 

1  dU 
r  cos  A  3e 


(4.4a) 

(4.4b) 


(4.4c) 


Substitution  of  U  into  equations  (4.4a),  (4.4b),  and  (4.4c)  and  expansion  of 
the  spherical’  harmonics  by  means  of  Eq.  (4.2)  is  tedious,  but  can  easily  be 
carried  out.  The  result  for  the  first  few  terms  in  the  series  is  shown  in 
order  to  provide  the  reader  with  some  idea  of  the  actual  results.  When  all 
the  substitutions  are  made,  Eq.  (4.3)  for  g  becomes 

g  =  v  ~  [l  +  (j)  j|  (3  sin2A-l)C2)04-  9  cos2A(C2>2cos20+  S2>2sin29)| 

+  |  2(5  sin3A-3  sinAJC-^o 

+  2(15  sin2A-3)cosA(C,  cos0+  S  .sin0) 

p 

+  60  sinAcos  A(C^  2cos20+  S^  2sin20) 

+  60  cos^A(C_  -cos30+  S_  _sin30)l  +  ....1 


n  ^3  sinAcosAC2  Q  -  6  sinAcosA(C2  2cos20+  S2  2sin20)J- 

+  (^)  |  |(15  sin2A-3)cosAC^  ^ 

o  2 

+  ^(ll  -  15  sin  A)sinA(C_  1cos0+  S,  ..sin0) 

d  Jy  1  J  J  J- 

p 

+  15(l  -  3  sin  A)cosA(C,  Ocos20+  S-,  _sin20) 

3>*-  3)*- 


+  45(sin  A-l)sinA(C_  ,cos30+  S_  _sin30) 
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(fr)  j^eosMSg  2cos20-  C2  2sin29^j' 

■  (15  sin2A-3)(S3jlcos0-  C3  jSino) 


+  30  sinAcosA(S2  2cos20-  C2  2sin20) 

+  45(l-sin2A)(S0  _cos30-  C_  _sin30) 
3y  3  JjJ 


}+ ....] 


(4.5) 

Although  the  forces  due  to  even  higher  ux’der  harmonics  will 
influence  the  orbit  of  an  artificial  satellite  to  some  extent,  the  presence  of 
terms  other  than  those  containing  C2  q  has  a  negligible  effect  upon  the  reentry 
of  a  spacecraft.  By  neglecting  all  iuch  terms,  the  complexity  of  the  gravita¬ 
tional  vector  is  considerably  reduced,  with  the  (e)  component  disappearing 
altogether.  The  result  is 


-  -  GM 

g  58  v  -p 
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Kr)! 


(3  sin2A-l)C2>0 


+  n  m 
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3  sinAcosAC, 


(4.6) 


The  three  remaining  constants  are  given  by 

GM  =  3.986032  x  1020  cm3/sec2 

or  {k‘7) 

GM  =  1.407654  x  10  ft3/ sec 

R  =  6.378196  x  108  cm 

or  0  (4.8) 

Ro  =  20,925,840  feet 

C  =  -1.08248  x  10"8  (4.9) 

dyU 

They  were  taken  as  the  best  average  of  all  the  values  suggested  by  Refs.  24  to 

26. 


(4.8) 


(4.9) 
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4.4  Oblate  Shape  of  the  Earth 

The  primary  departure  of  the  Earth  from  a  perfect  sphere  is  the 
flattening  of  the  poles.  Several  papers  by  Kaula  (Refs.  25  and  26  in  particular) 
have  also  indicated  the  additional  distances  by  which  the  Earth  departs  from 
its  ellipsoidal  shape,  which  in  turn  gives  rise  to  the  additional  coefficients 
in  Eq  .  (4.1).  The  variations  are  systematic  to  the  extent  that  they  indicate 
a  triaxial  tendency.  Wagner  (Refs.  27  to  29)  has  done  extensive  work  on  this 
subject  with,  however,  conflicting  results.  Consequently,  only  an  approximate 
(one  figure  accuracy)  average  was  adopted  with  the  major  axis  at  about  twenty 
degrees  longitude  west  and  a  minor  axis  about  210  feet  shorter  (fq  =  10"5) 
situated  at  seventy  degrees  longitude  east. 


In  order  to  calculate  the  perpendicular  height  above  a  triaxial 
ellipsoid,  it  is  sufficient  to  know  the  geocentric  geographical  position. 

The  three  variables  will  be  the  geocentric  radius,  the  geocentric  latitude, 
and  the  geocentric  longitude.  Several  different  methods  were  found  in  the 
literature  which  provided  suggestions  on  how  to  calculate  the  perpendicular 
height  above  an  oblate  spheroid.  While  most  were  unnecessarily  complex,  a 
very  simple  solution  is  given  by  Baker  (Ref.  30)  which  is  deemed  to  be  the 
best.  In  the  notation  of  this  paper,  the  suggested  equation  is 

y  -  r  -  Rq  [l  -  fp  sin2*  -  |  (fp)2  ^  sin22*J  (4.10) 

where  y  is  the  perpendicular  height  above  an  ellipsoidal  Earth. 


In  order  to  compensate  for  the  triaxial  tendency  of  the  Earth’s 
equator,  an  additional  term  of  the  form  R0  fq  sin20  (which  is  derived  in  the 
same  way  as  the  R0  fp  sin2*  term  for  the  polar  flattening)  is  required.  The 
higher  order  term  in  fq2  is,  of  course,  negligible  in  comparison.  This  additional 
term  must  then  be  multiplied  by  cos2*  in  order  to  reduce  the  triaxial  effect 
to  zero  at  the  poles.  A  biasing  constant  is  also  required  in  order  to  raise 
the  major  axis  and  lower  the  minor  axis  above  and  below,  respectively,  the 
Earth's  average  equatorial  radius,  R0.  This  biasing  constant(of  |)  in  equation 
(4.11)  is  required  due  to  the  definition  of  fq  as  given  in  the  nomenclature. 


In  addition,  sin22*  is  expanded  using  the  usual  trigonometric 
relationship  in  order  to  adapt  the  formula  for  digital  computation  (both  sin  * 
and  cos*  are  already  available,  but  sin2*  would  have  to  be  calculated).  The 
result  is 


y  =  r-  R  +  R  f  sin2*  +  (f  /r  +  f  )sin2*cos2* 
o  o  j)  r  c 


+  Ro  fq  {sin  (e  -  e2>2)  -  \  }cos2* 


(4.11) 


An  exact  procedure  which  required  several  iterations  to  yield 
the  result  and  which  was  not,  therefore,  useful  during  the  actual  computation 
of  trajectories  was  then  compared  with  the  answer^  from  Eq.  (4.11).  The  exact 
procedure  indicated  that  a  correction  factor  should  be  added  to  take  account  of 
higher  order  terms  which  were  omitted  from  the  power  series  approximation  given 
in  Eq.  (4.10).  The  actual  formula  used  during  calculation  of  trajectories  was 


y  =  r 


R0  +  R0  fp  sin2*  +  (fr/r  +  fc  +  ff(sin2* 
+  Rq  {sin2  (0-  02j2)  -  cos2* 


|) )  sin2*cos2* 


(4.12) 
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where 


or 


R0  =  6.378196  x  108cm 
R0  =  20,925,8^0  feet* 


fp  =  1/298.28 
f  =  10'5  0 


2,2 


=  -  20  degrees 


(4.13) 

(4.14) 

(4.15) 


were 

taken  as  the  best  average  of  all  the  values  suggested  by  Refs. 

24  to  29. 

Also 

fr  -  s(Eo  fp)2 

(4.16) 

f0  -  vy2'8 

(4.17) 

And 

ff  =  -3.328 

(4.18) 

was  optimized  for  the  adopted  values  of  RQ  and  fp. 

In  practice,  Eq.  (4.12)  evaluated  the  last  three  terms  in  single 
precision  arithmetic  and  then  used  double  precision  addition  when  adding  all 
five  values  together  to  achieve  the  final  result.  Comparison  with  the  exact 
iterative  procedure  sh 'wed  that  the  accuracy  achieved  with  Eq.  (4.12)  was  always 
better  than  0.05  feet  over  the  tested  altitude  range  of  zero  to  one  million  feet. 


4.5  Mean  Standard  Atmosphere 

In  this  investigation,  the  basic  model  used  to  define  the  varia¬ 
tion  of  density  with  altitude  was  the  U.S.  Standard  Atmosphere,  1962  (Ref.  31). 
As  far  as  corrections  for  oblateness  were  concerned,  the  true  perpendicular 
altitude  above  the  surface  was  used  to  define  the  height  as  per  Eq.  (4.12).  In 
essence,  the  atmosphere  was  assumed  to  be  a  uniformly  thick  blanket  enveloping 
the  Earth.  Figure  17  from  a  report  by  Cole  and  Kantor  (Ref.  32)  gives  ample 
evidence  to  support  his  assumption. 

Additional  information  regarding  standard  atmospheres  for 
different  latitudes  has  recently  appeared  in  the  literature  (Ref.  33).  However, 
in  view  of  the  lack  of  detail  concerning  the  space  between  each  mean  set  of 
values,  the  additional  accuracy  afforded  by  the  extra  referen'  .  curves  would 
not  have  been  worth  the  additional  programming  effort  necessary  to  make  them 
available  to  the  program  calculating  the  flight  trajectories.  In  addition, 
there  must  certainly  be  longitudinal  variations,  as  well  as  those  associated 
with  latitude,  and  the  current  absence  of  information  about  these  is  further 
justification  for  ignoring  both  effects.  Perhaps  the  most  important  point  in 
favour  of  using  a  single  standard  atmospheric  model  over  the  whole  Earth  is 
the  necessity  of  determining  the  effects  of  departures  from  a  non-rotating, 
circular  Earth  in  a  manner  which  distinguishes  them  from  those  due  to  atmos¬ 
pheric  density  variations. 
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4.6  Atmospheric  Density  Variations 

The  presence  of  spatial  density  variations  about  some  mean  values 
in  the  atmosphere  poses  a  problem  with  which  any  control  system  must  cope. 
Whether  the  purpose  of  the  control  system  is  to  achieve  lower  temperatures,  a 
terminal  position,  or  any  other  of  several  possibilities,  the  most  difficult 
part  of  the  problem  is  the  fact  that  the  density  variations  about  the  mean  are 
random  and  can  only  be  specified  by  way  of  statistical  data.  Since  this  paper 
deals  with  the  feasibility  of  the  COIfCAC  system  rather  than  with  the  exact  de¬ 
tails  of  the  hardware,  a  single  mean  curve,  the  U.S.  Atmosphere,  19&2,  was 
deemed  sufficiently  representative  for  all  but  a  few  reentries.  Then^  the  pro¬ 
blem  was  to  find  a  model  that  would  permit  a  reasonably  accurate  picture  of  the 
spatial  density  variations  to  be  found  in  the  atmosphere  about  the  mean  curve 
of  density  versus  altitude. 

One  question  that  was  deemed  important,  however,  was  whether  the 
shape  of  the  mean  curve  of  density  versus  altitude  was  itself  of  any  great 
importance.  In  order  to  answer  it,  a  number  of  quite  arbitrary  exponential 
curves  were  used  with  widely  varying  values  of  the  scale  height.  Actually, 
each  new  reference  curve  was  composed  of  two  exponentials.  The  higher  altitude 
portion  was  defined  by  the  scale  height  and  the  point  at  which  it  intersected 
the  U.S.  Standard  Atmosphere,  1962.  The  lower  curve  was  defined  by  the  cut-off 
point  for  the  upper  curve  and  a  second  point  of  intersection  with  the  U.S. 
Standard  Atmosphere,  1962  which  was  below  the  cut-off  point.  The  results  are 
given  in  Sec.  8.2. 

Two  types  of  spatial  density  variations  about  the  mean  were  used. 
One  was  based  on  altitude  and  the  other  was  generated  as  a  set  of  Gaussian 
random  numbers  which  were  spaced  at  equal  intervals  along  the  flight  path  length. 
Considerable  attention  was  paid  to  the  choice  of  these  two  types  due  to  the 
almost  complete  lack  of  data  regarding  spatial  density  variations  in  the  atmos¬ 
phere.  A  great  many  papers  could  be  found  which  gave  day  to  day  variations  at 
a  given  point.  Other  investigators  gave  the  results  of  calculations  which 
defined  new  mean  atmospheres  at  various  latitudes.  However,  Cole  and  Cantor 
(Ref.  32)  did  provide  some  information  regarding  the  maximum  spatial  density 
variations  about  the  mean  which  might  be  expected  as  well  as  the  distances  over 
which  the  departures  might  occur. 

For  the  density  variations  based  on  altitude,  an  actual  experi¬ 
mental  result  shown  in  Ref.  32  confirmed  a  suggestion  by  Etkin  (Ref.  34)  that 
a  simple  sinusoidal  variation  about  the  mean  of  density  versus  altitude  would 
be  realistic.  The  necessary  parameters  were  the  amplitude  of  the  variation 
about  the  mean,  the  wave  length,  and  the  point  of  zero  variation  or  the  nodal 
position.  A  range  of  values  were  used  and  the  results  can  be  found  in  Sec.  8.3. 

In  the  case  of  random  departures  from  the  basic  curve,  a 
statistical  approach  was  required.  The  method  chosen  to  simulate  the  variations 
about  the  mean  was  a  series  of  random  numbers  representing  fractional  deviations 
from  the  mean,  with  a  Gaussian  distribution.  After  specifying  the  variance  of 
the  series,  the  numbers  were  then  assigned  to  equidistance  positions  along  the 
flight  path.  In  the  absence  of  any  other  data,  linear  interpolation  was  used 
between  the  points  to  find  the  fractional  departure  from  the  mean  value  for  the 
density.  The  fraction  could  be  positive  or  negative,  but  was  limited  to  an 
absolute  magnitude  of  three  times  the  variance  (which  it  would  be  less  than  99$ 
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of  the  time  in  any  case)  in  order  to  prevent  bizarre  results  such  as  a  negative 
density.  In  order  to  provide  a  statistical  result,  a  number  of  trajectories 
was  computed,  one  for  each  of  several  different  series  of  random  numbers  having 
the  same  variance  and  correlation  interval.  (The  effect  of  the  random  density 
model  is  presented  in  Sec.  8.4.) 

While  the  above  method  of  simulating  the  random  variation  to  be 
expected  in  the  atmosphere  is  certainly  not  the  best  possible  model,  it  should 
at  the  very  least  be  adequate.  In  any  case,  the  lack  of  any  concrete  data 
upon  which  to  base  a  more  realistic  model  does  not  warrant  a  more  sophisticated 
approach.  It  is,  of  course,  possible  that  the  information  exists  in  classified 
reports.  If  so,  the  results  presented  in  Sec.  8.4  will  provide  an  immediate 
answer  as  to  whether  the  CONTAC  system  will  perform  in  a  satisfactory  manner 
under  all  expected  atmospheric  variations.  On  the  other  hand,  if  such  data 
actually  does  not  exist,  then  an  evaluation  of  the  CONTAC  system  must  wait 
until  it  has  been  determined. 


V.  EQUATIONS  OF  MOTION 
5.1  Basic  Reference  System 

The  choice  of  a  basic  reference  system  from  among  the  multitude 
of  possibilities  is  usually  dictated  by  the  problem  which  is  to  be  solved  as 
well  as  by  the  methods  of  solution  which  are  available.  In  many  situations,  as 
in  the  present  case,  the  choice  is  easily  made  due  to  the  many  restrictions  im¬ 
posed  by  ur favour able  reference  systems  which  might  require  an  unnatural  set  of 
equations  or  perhaps  extra  terms  in  the  equations.  Initially,  this  investigation 
was  begun  using  an  analogue  computer  of  limited  capacity.  When  the  problem  out¬ 
grew  the  available  equipment,  a  large  scale  digital  computer  was  substituted. 

The  original  reference  system  was  retained,  however,  due  to  the  background  of 
experience  and  confidence. 

The  basic  reference  system  chosen  is  a  set  of  axes  (l,J,K)  fixed 
at  the  mass  centre  of  a  rotating  Earth  (Fig.  18).  Thus,  the  position  of  t.  _ 
vehicle  is  naturally  expressed  in  terms  of  the  geocentric  latitude,  A,  the  geo¬ 
centric  longitude,  8(positive  to  the  East),  and  the  geocentric  radius,  r  (Fig. 
19);  the  velocity  of  the  vehicle  is  then  given  by  the  geocentric  azimuth,  ip 
(heading  with  respect  to  the  North-positive  to  the  East),  the  glide  angle,  7 
(angle  between  the  local  geocentric  horizontal  and  the  velocity  vector-positive 
downward),  and  the  magnitude  of  the  velocity  vector,  V  (Fig.  19).  However,  since 
all  aerodynamic  data  is  always  referenced  to  a  set  of  axes  which  are  aligned  in 
some'  way  with  the  vehicle,  a  second  reference  system  which  will  be  termed  flight 
axes  (i,j,k)  must  be  introduced  (Fig.  20).  It  should  be  noted  that  these  flight 
axes  are  not  fixed  to  the  body;  on  the  contrary,  the  (j)  axis  always  remains 
in  the  local  geocentric  horizontal  plane,  thereby  forcing  the  (T,kj  axes  to  re¬ 
main  in  the  local  geocentric  vertical  plane.  In  order  to  complete  the  orienta¬ 
tion  requirements,  the  [i]  axis  is  directed  along  the  velocity  vector,  V  ,  and 
the  (k]  axis  is  directed  downward  (in  this  study  I7I  <  90  degrees).  As  a  con¬ 
sequence  of  this  definition  for  the  flight  axes,  the  side  force  equation  will 
include  the  component  of  the  lift  in  that  direction,  namely:  L  sin  <p  .  For 
those  readers  who  have  not  encountered  this  vehicle  reference  system  of  flight 
axes  (for  example,  Refs.  35  and  36),  it  should  again  be  mentioned  that  this 
investigation  is  concerned  with  a  three  degree  of  freedom  system  of  equations 
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for  a  point  mass  vehicle.  Since  the  pitch,  roll,  and  yaw  angles  which  are  de¬ 
pendent  variables  in  a  six  degree  of  freedom  system  are  now  control  variables, 
it  is  rather  useful  to  simplify  the  three  force  equations  by  using  the  required 
components  of  the  lift  and  drag.  The  point  mass  vehicle  is  then  rotated  in  pitch 
and  roll  (sideslip  is  assumed  negligible  and  the  yaw  angle  is,  therefore,  set 
identically  equal  to  zero)  so  as  to  develop  these  required  lift  and  drag  com¬ 
ponents.  The  validity  of  the  assumption  regarding  the  representation  of  aero¬ 
dynamic  forces  has  been  examined  in  Section  III. 


5.2  Conversion  Between  Reference  Systems 

Since  the  equations  of  motion  are  written  in  the  flight  axes 
reference  system,  it  is  necessary  to  find  the  components  of  all  forces  in  this 
final  reference  system.  (Fig.  21).  Thus  the  conversion  matrices  from  any  of 
the  basic  or  intermediate  reference  systems  to  the  flight  axes  reference  system 
must  be  found._  A  rotation  about  the  _£^}_axis  (Fig.  22)  gives  the  transformation 
from  the  (l,J,K}  set  of  axes  to  the  (l',J',K}  set  of  axes.  The_sj;cond  rotation 
about  thefj'Jaxis  (Fig.  23)  gives  the  transformation  to  the  {n,e,v}  set  of  axes. 
The  third  rotation  about  the  {v}  axis  (Fig.  24)  gives  the  transformation  to  the 
set  of  axes.  Finally,  the  fourth  rotation  about  the  (j)  axis  gives  the 
transformation  to  the  (T,j,k)  reference  system  of  flight  axes  (Fig.  25).  These 
four  rotations  naturally  require  the  four  simple  conversion  matrices  given  with 
each  figure.  The  overall  conversion  matrices  from  any  of  the  intermediate  re¬ 
ference  systems  sire  then  easily  found  to  be 
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The  conversion  matrix  between  the' basic  reference  system  (I,J,K) 
and  the  flight  axes  reference  system  (i,j,kj  could  be  found  and  expressed  in 
the  same  manner;  however,  it  is  unnecessary  since  it  is  not  used. 

5.3  Vector  Equation  of  Motion 


An  evaluation  was  made  of  the  relative  magnitudes  of  all  forces 
which  may  act  upon  a  spacecraft  during  a  reentry  into  the  Earth's  atmosphere. 
Forces  due  to  other  planetary  bodies  including  the  moon,  the  solar  wind,  and 
magnetic  fields  often  have  an  appreciable  effect  on  the  orbits  of  particular 
satellites  when  these  forces  are  experienced  over  long  periods  of  time.  However 


for  a  reentry  trajectory,  it  can  easily  be  shown  (Ref.  38)  that  only  aerodynamic 
forces  and  the  Earth's  gravitational  field  have  a  non~negligible  effect  upon  the 
flight  path.  In  addition,  based  on  the  information  in  Ref.  23,  it  may  be 
assumed  that  the  basic  reference  system,  (l,J,K)  (see  Fig.  18),  is  rotating  in 
inertial  space  at  a  constant  angular  velocity. 

The  single  vector  equation  of  motion  which  is  required  for  a 
point  mass  moving  in  such  a  rotating  reference  system  is  given  by  (Ref.  39) 

=  a  +  2(wxV)  +  w  x(wxr)  (5.*+) 

where  a  is  the  acceleration  relative  to  the  (l,J,K)  reference  system,  i.e. 
relative  to  the  Earth. 

The  only  external  forces  are  either  aerodynamic  or  gravitational. 
It  is  appropriate,  therefore,  to  separate  the  left  hand  side  of  Eq.  (5.4)  and 
write  _  __  _ 

F  =  FA  +  mg  (5.5) 

The_relative  acceleration,  a,  in  the  Cl,J,K)  system  has  compon¬ 
ents  along  the  (i,j,k)  axes  that  are  obtained  by  the  rule  for  transforming  a 
derivative  between  rotating  frames  of  reference  (Ref.  40),  i.e. 

a  =  V  +  (ftxV)  (5.6) 

where  ft  is  the  angular  velocity  of  the{i, j ,k]  system  relative  to  the  (T, J,K)  and 
the  components  of  V  are  the  rates  of  change  of  the  components  of  V  along  the 
the  {i,j,k)  axes. 

Equations  (5*4)  to  (5-6)  can  be  combined  to  give 
F,  .. 

—  +  g  =  V  +  (ftxV)  +  2(wxV)  +  wx(wxV)  (5.  7) 


5.4  Expansion  of  the  Single  Vector  Equation  of  Motion 

The  expansion  of  each  of  the  terms  in  Eq.  (5.7)  in  the  {i,j,k} 
flight  axes  system  is  not  difficult  t.lthough  it  is  somewhat  involved  and 
lengthy.  The  three  components  of  each  term  will  first  be  found  in  whichever 
axes  system  is  most  convenient  then,  if  necessary,  transferred  to  the  {i,j,k} 
flight  axes  system,  and  finally  collected  to  give  the  three  complete  force 
equations.  Starting  with  the  first  term,  the  results  are  as  follows. 

The  aerodynamic  forces  are  given  by 
FA  =  -  I  CD(a)  |p(y)  V2  S 

+  I  CL(a)  |p(y)  V2  S  sin  9  *5-8) 

-k  CL(a)  gp(y)  V2  s  cos  9 

CL(a)  and  CD(oO  are  given  by  Eqs.  (3*3)  and  (3*4).  The  expression  for  the 
density,  p(y),  has  been  discussed  in  Sec.  4.5  and  the  evaluation  of  the  height 
y,  above  a  triaxial  ellipse  is  given  by  Eq.  (4.12). 
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_  __  The  gravitational  force  vector  is  given  by  Eq.  (4.3)  in  the 

5,e,v}  reference  system._  Transforming  the  components  from  this  intermediate 
reference  system  to  the  {i,j,k}  flight  axes  system  via  equation  (5*2)  gives 
the  gravity  vector  as 

g  =  T  (-g  siny  +  g.  cosycos^  +  g  cosysirr ip) 

+  1  (-  sin^  +  gecos^) 

+  k  (-g  cosy  -  g.  sinycos if/  -  gQsinysin^)  (5*9) 

r  a  u 

The  expressions  for  gf,  g^,  and  gg  are  developed  and  stated  in  Sec.  4.3. 
Although,  as  noted  in  Sec.  4.3,  g@  is  too  small  to  affect  the  reentry  trajectory 
of  a  spacecraft,  the  expressions  have  purposely  been  left  in  a  generalized 
form  so  that  the  precise  effect  of  any  term  in  the  series  may  be  ascertained 
should  it  be  desirable  at  some  future,  date.  In  addition,  the  values  of  gr  and 
g^  will  normally  include  only  the  major  effect  of  the  Earth's  oblate  gravity 
potential,  namely  C2  0  term.  However,  here  again,  the  expressions  are  left  in 
their  generalized  form.  The  truncated  expressions  which  were  actually  used  in 
the  calculations  may  be  found  by  a  comparison  of  Eqs.  (4.3)  and  (4.6)  to  be 
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[l  +  (gf  |  (3  sin2X-l)C2_0] 

(5.10) 

= 

(24 

2 

r 

[(j)  3  slnXeosM2j0j 

(5.U) 

ge  = 

0 

(5.12) 

The  terms  on  the  right  hand  side  of  Eq.  (5-7)  are  either  very 
simple  or  quite  complex.  The  first  term  is  simply 

•  • 

V  =  I  V  (5.13) 

Before  proceeding  with  the  next  term,  the  expression  for  ft 
should  be  found.  With  reference  to  Figs .22  to  25,  it  can  easily  be  seen  that  _ 
the  sum  total  of  all  the  rotations  which  are  required  to  rotate  from  the  {i, J,K) 
basic  reference  system  to  the  {i, j ,k}  flight  axes  reference  system  yields  the 
relative  rotational  velocity  which  is 

ft  =  K6+J'A  +  v  ^  -  j  y  (5.14) 

The  negative  sign  in  the  last  term  is  required  due  to  the  rotation  of  h  into  v 
which  is  in  the  opposite  sense  from  the  normal  right  hand  rule. 

At  this  point,  the  three  kinematic  relationships  for  the  position 
of  the  vehicle  r,  A  and  0  should  be  defined  since  the  last  two  are  required 
immediately.  They  are 


• 

-V  sin  y 

(5.15) 

• 

A  = 

X  cosycos^ 

(5.16) 

• 

e  = 

V  cosysin^ 
r  cosA 

(5.17) 

Then,  in  order  to  transform  ft  into  the  (i,j,k)  flight  axes  system,  the  trans- 
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formation  given  in  Eqs.  (5*2)  and  (5.3)  oust  be  used,  along  with  the  substitu¬ 
tion  of  Eqs.  (5.l6)  and  (5.17).  The  result  is 


_  _  •  y 

ft  =  i(^sin7  -  —  CqS^  sinycosysin^sin^) 

+  j  (-7  -  ~  COS7) 

_  •  y  2 

+  k(^cos7  -  — ^^cos  7sinXsin^) 


(5.18) 


Finally,  combining  the  expression  for  ft  with  V  =  i  V  yields  the  desired  cross 
product  for  the  second  term  on  the  right  hand  side  of  Eq.  (5.7).  The  result  is 

(ftxV)  =  I  (0) 

—  *  y2  2 

+  j  tyV  COS7  -  r  cosX  COS  7sin^s^nV') 

-  .  v2 

+  k  (7  V  +  -j- cos  7)  (5.19) 

The  next  term  is  somewhat  easier  to  express.  It  is  easily  seen 

(Fig.  23)  that 

CJ  =  K  w  (5.20) 

which  upon  substitution  of  Eq.  (5*3)  and  combination  with  V  =  i  v  results  in 
2(wx7)  =  i  (0) 

-  j(2uV[sin7cos?\cos^+  cosysinA]) 

+  k  (2uV  cosAsin^)  (5* 21) 

Finally,  the  last  term  is  easiest  to  expand  in  the  reverse 
direction.  It  is  also  easy  to  see  (Fig.  23)  that 

r  =  - v  r  (5.22) 

or  ______ 

r  =  I'r  cosA  -  K  r  sin?v  (5.23) 

which  upon  combination  with  Eq.  (5.20)  gives 

w  x  (u  x  r)  =  I’  u2  r  cosX  (5»2^) 

and  after  substitution  of  Eq.  (5.3)  results  in 

w  x  (uxr)  =  i  (w  rcosA[sin7cosA  +  cos7sin^cos^]  ) 

-  j  (u  rcos?\sinXsin^) 

+  k  (u  rcosX[cos7sinA  -  sin7sinAcos^]  )  (5*25) 

_  All  of  the  terras  in  Eq.  (5«7)  have  now  been  expanded  into  the 

(i,j,k)  flight  axes  reference  system.  Since  each  of  the  three  components  are 
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independent  of  the  other,  three  first  order,  non- linear, differential  equations 
may  be  written  by  collection  of  all  the  terms  in  Eqs.  (5.8),  (5.9)  (5  r?) 

(5.19),  (5.21),  and  (5.25).  The  result  is  ’  ’  ^  h 


•  2 

^  =  “  Cjj(ot)gp(y)V  S/m  -  g^siny  +  g^cosycosij/  +  ggCosysin^ 
2 

-w  rcosA(sin7cosA  +  cosysinAcos^) 


7V  =  -  CL(a)^p(y)V  —  cos<p  -  C0S7  -  g^osy  -  g^sinycos'^ 

2 

“  ggSinysin^  -  w  rcosA(cos7CosA  -  sinysinycos^/) 

-  2u>VcosAsin^ 

*v°°s?  ■  f  si*p  +  y^Z^sisH  .  vi#  +  e$cosy, 

2 

+  (o  rcosAsinAsinf  +  2<*)V(sin7cosAcos^+  cos7sinA)  (5.28) 

..  _  ,  In  addition  to  the  above  three  differential  equations,  and 

the  positional  differential  equations  (5.15),  (5.16),  and  (5. 17),  several  other 
relationships  are  required  which  have  already  been  stated  or  developed  previous- 
^'o\F0r1c?”vf!?ience»  th*y  wil1  be  mentioned  again  at  this  point.  Equations 
(3.3)  and  (3 ^0  give  the  relationship  between  the  lift  and  drag  coefficients, 
CL(a)  and  CD(a),  and  the  angle  of  attack, a.  The  relationship  between  the 
density,  p(y),  and  the  perpendicular  height  above  the  Earth's  surface,  y,  has 
been  discussed  in  Sections  4.5  and  4.6.  The  perpendicular  height  above  the 
sm-face  of  a  triaxial  Earth,  y,  is  given  by  Eqs.  (4.12)  to  (4.18).  Finally, 

Si*  f°r  gA’  and  g0  are  «iven  by  Eqs.  (5.10),  (5.11),  and  (5.12). 

Sbo^d  b®  desirable  to  include  additional  terms  in  the  oblate  gravity  field 

°f  equations  (4.1)  to  (4.5)  will  serve  to  illustrate  the  manner  by 

which  it  may  be  done. 


(5.26) 


(5.27) 


At  this  point  in  the  development  of  the  equations  of  motion, 
it  was  found  to  be  convenient  to  make  a  change  in  the  independent  variable  from 
time,  t,  to  path  length  aiong  the  trajectory,  s.  The  original  reasons  for  mak¬ 
ing  this  change  are  no  longer  a  justification  for  it.  However,  it  was  con¬ 
venient  to  maintain  the  new  independent  variable  and  the  following  first  order, 
non-linear,  differential  equations  were  the  ones  which  were  actually  integrated. 
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The  three  force  equations 
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The  three  positional  equations  become 
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Also,  the  differential  equation  for  the  time  is 

dt  1 
ds  "  y 


5-5  Numerical  Procedures 


(5.35) 


A  major  problem  in  the  development  of  digital  computer  programs 
is  the  problem  of  retaining  accuracy  during  the  course  of  the  computation  while 
keeping  the  number  of  points  (and  therefore  the  computation  tioe)  as  low  as 
possible.  In  particular,  the  calculation  of  flight  trajectories  of  reentry 
vehicles  often  requires  an  accuracy  of  six  or  eight  figures  in  the  altitude  or 
geocentric  radius,  respectively.  The  latter  requirement  certainly  exceeds  the 
precision  of  numbers  in  most  scientific  computers  when  round-off  error*  inherent 
in  all  integration  schemes  is  taken  into  account.  Truncation-error  is  easily 
reduced  by  using  one  of  the  standard,  high-accuracy  integration  techniques  such 
as  the  fourth  order  Runge-Xutta  method  and  the  proper  size  of  interval.  Eovever, 
round-off  error  proves  much  more  difficult  to  eliminate,  especially  when  a 
small  enough  interval  size  (and  therefore  a  large  number  of  points)  is  used  to 
reduce  the  truncation- error.  Consequently,  a  practical  method  oust  be  found  for 
reducing  round-off  error  and  also  of  reducing  the  number  of  integration  points 
while  simultaneously  keeping  the  truncation-error  to  a  minimum. 


*  Error  is  introduced  from  two  sources: 

(a)  Truncation-error  cay  be  defined  as  the  difference  between  the  true  value 
and  the  value  obtained  from  the  numerical  integration  formula  which  makes 
use  of  a  finite  number  of  teres  from  an  infini*-*  Kaclaurin’s  or  Taylor's 
series . 

(b)  Round-off  error  is  the  result  of  the  limited  accuracy  (usually  8  figures) 
available  within  the  computer. 
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The  solution  to  the  problem  of  round-off  error  took  advantage 
of  the  basic  characteristic  of  digital  calculations.  All  numerical  integration 
schemes  involve  adding  small  increments  to  the  variables  in  order  to  find  the 
next  value  after  each  step.  Accordingly,  these  increments  possess  far  more 
absolute  accuracy  than  their  respective  variables.  Thus,  if  the  variables 
could  take  advantage  of  this  increased  absolute  accuracy  when  incrementation 
took  place,  the  round-off  error  would  be  greatly  reduced.  Notice  that  the 
incrementation  stage  would  be  the  only  one  at  which  the  increased  accuracy 
(and  therefore  greater  computation  time)  would  be  used  since  the  calculation 
of  the  increments  would  use  the  variables  only  to  their  original  accuracy. 

In  practical  computational  terms  (using  popular  languages  such 
as  FORTRAN),  the  scheme  merely  involves  the  use  of  double-precision  addition 
whenever  the  increments  are  added  to  the  variables.  Normal  single-precision 
computation  of  the  increments  makes  use  only  of  the  most  significant  part  of 
the  variables.  In  addition,  an  analysis  of  the  whole  program  should  be  made 
as  other  sections  may  benefit  from  this  technique.  An  example  is  the  calcu¬ 
lation  of  the  altitude  over  an  oblate  Earth  when  only  the  geocentric  latitude 
and  the  geocentric  radius  are  available.  Having  calculated  the  "corrections" 
in  single  precision,  they  should  be  added  in  double  precision. 

A  major  reduction  in  the  number  of  steps  is  also  possible,  by 
using  a  continuously  varying  interval  step  size  dictated  by  truncation -error. 

It  has  been  found  that  the  truncation- err or  depends  on  how  quickly  the  variables 
are  changing  as  well  as  the  step  size.  In  the  calculation  of  flight  trajectories 
of  entry  vehicles  using  geographical  coordinates,  the  key  variables  are  the  glide 
angle,  the  azimuth,  and  the  latitude.  Experience  has  shown  that  whenever  all 
of  these  variables  are  reasonably  stationary,  the  interval  step  size  can  be 
rather  large.  The  criterion  which  was  used  to  determine  the  interval  size 
(within  a  certain  preset  range)  was  that  the  change  in  all  these  angles  per 
step  should  be  kept  below  a  certain  minimum  value.  By  varying  the  interval 
size  parameter,  the  value  could  be  f-~  .1  for  which  the  truncation-error  was 
accepta>"  ..e,  since  the  double-precision  incrementation  always  reduced  the  total 
round-off  error  to  a  negligible  value. 

Two  basic  methods  were  used  to  check  the  results  of  the  computa¬ 
tions.  First,  the  atmosphere  was  removed  (by  setting  the  density  to  zero)  and 
orbits  were  calculated  for  all  the  situations  for  wbi'h  simple  analytical  solu¬ 
tions  were  known.  Circular  and  elliptical  orbits  were  compared  over  non¬ 
rotating  and  rotating  Earth  reference  systems.  In  every  case,  the  additional 
terms  involving  the  Earth's  rotation  produced  the  proper  result  when  changes 
were  expected.  On  the  other  hand,  parameters  such  as  the  orbital  period  and 
the  major  axis  (for  elliptical  orbits)  remained  the  same.  Typical  accuracies 
were  six  or  seven  figures  over  several  orbits.  These  latter  results  also 
gave  confidence  in  the  numerical  integration  procedure  referred  to  above.  The 
check  on  the  accuracy  of  Eq„  (4.12)  was  made  via  the  test  program  developed  to 
check  the  original  validity  of  its  results. 

In  order  to  check  the  accuracy  of  the  digital  program  when  the 
atmosphere  was  present,  the  terms  due  to  rotation  and  oblateness  were  removed. 

The  equations  were  then  further  simplified  (the  atmosphere  was  made  exponential, 
only  zero  roll  angles  were  considered,  and  gravity  was  held  constant  at  its 
average  value)  and  the  results  compared  with  the  output  of  the  identical  equa¬ 
tions  which  had  been  wired  into  an  analogue  computer.  Over  the  range  for  which 
the  density  approximation  on  the  analogue  computer  was  reasonably  accurate  (from 
150,000  feet  up  to  300/500  feet)  there  was  no  discernable  difference  between  the 
two  results. 
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VI.  RESULTS  AND  DISCUSSION  FOR  A  CIRCULAR,  NON-ROTATING  EARTH 
6.1  Introduction 

The  next  three  sections  of  this  paper  present  an  analysis  of  the 
perl  fiance  of  the  CONTAC  system  with  results  obtained  from  the  digital  computer 
program  listed  in  Appendix  C.  Section  VI  is  devoted  to  presenting  the  overall 
capabilities  of  the  CONTAC  system  as  well  as  its  limitations. 

The  method  of  presenting  the  results  has  been  discussed  in  Sec. 
1.3.  Apart  from  the  peak  temperature,  the  two  main  quantities  of  interest  are 
the  equivalent  down  range  and  cross  range.  These  values  are  illustrated  and  de¬ 
fined  in  Fig.  3. 


In  order  to  give  the  reader  a  better  understanding  of  the  way 
in  which  the  CONTAC  system  works,  a  detailed  history  of  all  important  variables 
is  given  in  Figs,  lib  and  11c.  Figure  11a  is  a  similar  history  of  an  uncontroll¬ 
ed  reentry  and  may  be  used  as  a  comparison.  It  is  quite  c2ear  that  the  CONTAC 
system  is  easily  able  to  control  the  temperature  and  prevent  the  additional  rise 
of  65OF0  which  occurred  for  the  uncontrolled  reentry. 

In  order  to  provide  the  details  of  the  way  in  which  the  altitude 
angles  are  controlled,  an  expanded  scale  has  been  used  for  a  portion  of  the 
reentry  in  Fig.  lid.  Notice  first  that  the  line  of  time  vs.  down  range  is 
straight  over  this  small  portion  of  the  reentry.  Consequently,  the  abcissa  of 
the  graph  has  been  labelled  as  time  in  addition  to  down  range  in  order  to  pro¬ 
vide  a  better  feeling  for  the  time  scale  of  the  cycling  of  the  angle  of  attack 
by  the  CONTAC  system.  The  most  important  information  to  be  gained  from  Fig.  lid 
is  the  rates  involved  in  the  process.  Notice  that  the  fuel  economy  will  be 
quite  high  due  to  infrequent  use  of  the  reaction  jets. 

6.2  Footprint  Optimization 

The  standard  reentry  conditions  at  350,000  feet  altitude  were 
chosen  to  be  26,000  feet  per  second  for  the  speed  and  1.0  degree,  for  the  glide 
angle.  It  might  be  noted  that  the  speed  is  between  250  and  300  feet  per  second 
greater  than  the  orbital  velocity  at  this  height,  which  therefore  allows  the 
spacecraft  the  possibility  of  having  reached  this  position  from  an  altitude  of 
between  250  and  300  miles.  A  slight  increase  in  the  initial  speed  or  glide 
angle  would  allow  higher  initial  orbital  altitudes  up  to  500  miles  if  that  con¬ 
dition  were  desirable. 

Two  criteria  were  used  in  determining  the  best  flight  path  which 
could  be  obtained  from  the  initial  conditions.  The  most  important  factor  was 
to  prevent  throughout  the  reentry,  the  occurrence  of  a  temperature  greater  than 
the  peak  value  experienced  during  an  equilibrium  glide.  In  general,  this  re¬ 
quirement  was  slightly  relaxed  to  allow  temperature  peaks  which  were  greater  by 
20F°  or  30F°.  The  second  principal  requirement  was  to  obtain  as  large  a  cross 
range  as  possible  throughout  the  whole  extent  of  the  down  range  capability. 

In  addition,  some  effort  was  nded  to  increase  the  available 
down  range  increment.  It  is  here  that  a  further  reason  for  the  choice  of 
26,000  feet  per  second  as  the  initial  velocity  was  found.  On  Fig. 9,  it  is 
shown  that  the  speed  lost  up  to  the  point  of  the  first  temperature  maximum  can 
easily  be  controlled  to  a  value  between  150  and  /'50  feet  per  second  by  the  use 
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of  the  proper  entry  angle  of  attack.  As  noted  in  Sec.  1.9>  this  reduction  in 
the  vehicle's  speed  at  the  end  of  Phase  II  can  be  produced  without  causing  a 
significant  difference  in  altitude,  glide  angle,  or  down  range  distance  at 
that  position.  This  condition  results  from  the  variation  of  the  drag  coeffic¬ 
ient  with  the  angle  of  attack  while  the  lift  coefficient  remains  essentially 
constant.  Consequently,  due  to  the  choice  of  speed  at  the  initial  position, 
the  speed  at  the  end  of  Phase  II  is  either  very  slightly  above  (when  low  angles 
of  attack  with  relatively  low  drag  coefficients  were  used)  or  just  below  (when 
high  angles  of  attack  with  much  larger  drag  coefficients  were  used)  circular 
orbital  velocity.  In  a  vacuum,  the  former  speed  results  in  a  perigee  condition 
which  would  produce  an  additional  revolution  around  the  Earth  yielding  360 
degrees  of  down  range  increment  before  the  spacecraft  descends  to  the  same 
altitude.  On  the  other  hand,  even  in  a  vacuum,  the  latter  speed  would  produce 
an  apogee  condition  which  brings  about  an  immediate  decrease  in  altitude. 
However,  it  is  quite  obvious  that  the  presence  of  atmospheric  drag  will  tend  to 
alter  the  situation.  Except  when  complicated  by  Earth  oblateness,  even  the 
use  of  a  very  small  angle  of  attack  (at  5°>  Cj)  ~  Cp  ^n)  for  the  former  situa¬ 
tion;  results  only  in  an  extended  glide  at  practically  constant  altitude  (the 
rise  is  limited  to  30,000  feet)  until  the  speed  falls  below  the  circular  tailue 
by  about  1000  feet  per  second.  On  the  other  hand,  the  continued  use  of  a 
high  angle  of  attack  for  the  latter  situation  reduces  the  speed  very  quickly 
resulting  in  very  little  incremental  down  range.  In  both  extreme  situations 
as  well  as  any  between,  the  incremental  down  range  which  is  available  by  means 
of  subsequent  attitude  control,  is  very  nearly  the  same  once  the  speed  falls 
below  the  circular  value  by  about  1500  feet  per  second.  Therefore,  any  small 
speed  variations  at  the  end  of  Phase  II  on  the  order  of  one  hundred  or  even 
fifty  feet  per  second  will  result  in  large  variations  in  incremental  down  range 
during  Phase  III  of  the  reentry.  As  a  consequence,  the  situation  provides  the 
guidance  system  with  the  capability  of  large  down  range  increments,  without 
interference  from  the  CONTAC  system,  due  to  the  possibility  of  an  extended 
glide  nt  very  low  angles  of  attack  during  which  the  speed  is  lost  very  slowly. 

It  should  be  noted  that  the  above  discussion  is  not  altered  by 
the  inclusion  of  Earth  rotation.  The  immunity  is  due  to  the  fact  that  although 
the  relative  speed  can  change  by  up  to  -  1500  feet  per  second  at  the  equator 
for  initially  due  West  or  due  East  reentries,  respectively,  the  relative  cir¬ 
cular  orbital  speed  also  changes  by  exactly  the  same  value.  For  other  initial 
azimuths , "the  same  relative  change  also  occurs.  Consequently,  the  same  perigee 
and  apogee  conditions  result  upon  the  use  of  low  and  high  entry  (i.e.,  during 
Phase  II)  angles  of  attack,  respectively. 

With  regard  to  satisfying  the  first  criterion,  the  CONTAC  system 
behaved  very  well.  The  operational  limits  were  difficult  to  obtain  since  im¬ 
practical  situations  had  to  be  imagined  before  any  degredation  was  found  in  the 
peak  temperature  performance.  The  discussion  of  the  parameters  which  were 
varied  during  the  reentry  is  confined  to  the  attitude  at  this  particular  point; 
the  result  of  the  differences  or  errors  in  initial  conditions  is  presented  in 
Sec.  6.3. 


The  entry  angle  of  attack  was  limited  to  a  maximum  of  65  degrees 
in  order  to  prevent  excessive  heating  rates  at  the  rear  of  the  spacecraft.  It 
may  be  questionable  whether  larger  pitch  angles  could  be  maintained  due  to 
problems  of  stability  and  control.  However,  assuming  that  these  problems  can 
be  solved,  there  is  no  reason  why  angles  as  high  as  75  or  even  80  degrees  cannot 
be  used. 
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The  lower  limit  for  the  entry  angle  of  attack  was  approximately 
35  degrees.  When  smaller  values  were  used,  the  first  maximum  temperature  was 
too  large,  since  too  little  lift  was  developed  to  prevent  the  vehicle  from  des¬ 
cending  too  low  on  the  initial  penetration  into  the  atmosphere. 

The  entry  roll  angle  was  always  held  at  zero  during  this  investi¬ 
gation.  It  is  probable  that  a  moderate  roll  angle  during  Phase  II  could  in¬ 
crease  the  cross  range  without  affecting  the  peak  temperature  control.  Further 
investigation  into  this  question  might  prove  useful. 

The  skip  angle  of  attack  (which  is  controlled  by  the  guidance 
system)  is  employed  as  the  lower  limit  during  Riase  III.  The  largest  value  was 
the  lesser  of  55  degrees  and  the  entry  angle  of  attack.  The  lowest  value  used 
was  5  degrees.  The  guidance  system  might  employ  even  smaller  values  down  to 
zero  if  the  upper  surface  of  the  spacecraft  is  not  subject  to  or  can  cope  with 
large  heating  rates.  It  should  be  noted  that  when  the  large  entry  angles  of 
attack  of  55  or  65  degrees  were  used,  the  CONTAC  system  often  prevented  the 
vehicle  from  assuming  the  skip  angle  of  attack.  As  soon  as  it  commences  active 
operation  at  the  beginning  of  Phase  III,  the  CONTAC  system  disallows  lower 
limits  if  they  would  permit  too  high  a  peak  temperature.  Again,  it  is  noted 
that  this  operation  is  completely  automatic. 

The  skip  angle  of  roll  (which  is  also  controlled  by  the  guidance 
system)  is  employed  as  the  lower  limit  during  Phase  III.  The  "smallest"  value 
was  a  full  minus  ninety  degrees  to  the  left  (or  plus  ninety  to  the  right). 

Smaller  values  were  not  used  since  they  resulted  in  less  cross  range  which  was 
not  in  accordance  with  the  second  criterion.  The  "largest"  value  used  was 
zero  degrees.  The  CONTAC  system  often  prevented  the  vehicle  from  assuming  the 
skip  bank  angle  when  it  was  near  ninety  degrees;  the  conditions  which  led  to 
this  situation  were  the  same  as  those  which  prevented  very  small  skip  angles  of 
attack  from  being  used  and  again  the  CONTAC  system  automatically  disallowed  the 
value  and  prevented  the  temperature  from  achieving  too  high  a  value. 

The  upper  limit  on  the  angle  of  attack  during  both  Phases  III 
and  IV  was  55  degrees.  This  value  was  fixed  at  the  attitude  which  yielded  the 
maximum  vertical  force.  The  upper  limit  for  the  bank  angle  was  also  the  value 
which  yielded  maximum  vertical  force,  namely  zero  degrees.  The  reasons  for 
this  choice  have  already  been  given  in  Sec.  1.10. 

The  lower  limit  for  both  the  angle  of  attack  and  the  bank  angle 
did  not  affect  the  operation  of  the  COKTAC  system.  Therefore,  the  choice  of 
the  values  used  will  be  discussed  under  the  second  criterion  of  cross  range 
optimization. 

While  operating  within  the  above  values  for  the  attitude,  the 
CONTAC  system  performed  the  task  of  restraining  the  peak  temperature  experienced 
during  the  reentry  wit1*  rut  any  degredation  in  performance  for  any  combination  of 
these  values.  As  shown  i  Fig.  26,  the  peak  temperature  experienced  during  the 
reentry  varied  only  -  20  Fahrenheit  degrees  about  the  peak  value  experienced 
during  an  equilibrium  glide.  At  the  same  time,  cross  range  was  practically 
constant  (between  9  degrees  and  11  degrees)  throughout  most  of  the  available 
down  range.  Therefore,  it  is  possible  to  obtain  an  "isothermal  footprint" 

(with  regard  to  the  peak  temperature  experienced  during  the  reentry)  of  consider¬ 
able  area  when  the  CONTAC  system  is  employed  to  constrain  the  temperature  during 
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the  flight.  The  dimensions  are  over  one  hundred  degrees  equivalent  longitude 
in  down  range  and  about  twenty  degrees  equivalent  latitude  in  cross  range. 
Considering  the  fact  that  the  continental  United  States  has  a  "length"  of  only 
forty  degrees  in  equivalent  longitude  and  is  less  than  eighteen  degrees  of  equi¬ 
valent  latitude  wide,  this  "isothermal  footprint"  is  about  three  times  the  area 
of  that  country.  Consequently,  during  any  given  12  (or  at  most  24)  hour  period, 
enough  range  control  is  available  when  using  the  CONTAC  system  to  land  at  any 
point  within  or  near  the  United  States,  as  long  as  the  terminal  position  is 
underneath  or  at  least  near  the  orbital  envelope. 

The  attainment  of  the  second  criterion,  the  maximum  possible 
cross  range,  was  attempted  by  programming  the  skip  bank  angle  (Phase  III)  and 
the  lower  limit  for  the  angle  of  roll  (Fnase  IV)  within  the  restraints  imposed 
by  the  first  criterion.  During  each  Phase  of  the  reentry,  several  possibilites 
were  investigated  until  the  best  overall  performance  was  attained  throughout 
the  flight. 


Note  that,  the  final  results  for  cross  range  presented  in  this 
paper  should  not  be  considered  the  maximum  possible  values  since  it  was  not 
feasible  to  consider  every  parameter  combination.  In  addition,  the  calculation 
was  terminated  at  approximately  100,000  feet  where  the  speed  was  about  3000  feet 
per  second.  As  a  consequence,  about  3  degrees  of  equivalent  latitude  or  about 
200  miles  of  additional  cross  range  is  available  before  the  final  touchdown. 
However,  as  pointed  out  in  Section  1.12,  this  final  portion  of  the  flight  has 
not  been  .included  in  the  investigation.  Therefore,  no  details  for  Phase  V 
can  be  given. 


On  the  other  hand,  density  variations  in  the  atmosphere  and  other 
sources  of  error  may  reduce  the  anticipated  cross  range.  It  was  invariably 
found  that  when  the  relationships  which  yielded  maximum  cross  range  were  in¬ 
correctly  adjusted,  the  result  was  not  only  a  decrease  in  cross  range,  but 
also  either  increased  or  decreased  down  range  depending  on  the  direction  of  the 
mismatch.  Therefore,  it  may  be  necessary  to  sacrifice  a  portion  of  the  cross 
range  capability  in  order  to  retain  a  degree  of  down  range  mobility  in  the 
latter  phases  of  the  flight. 

At  this  point,  a  detailed  description  of  the  methods  used  to 
maximize  the  cross  range  will  be  given  for  Phases  I  to  IV. 

During  Phase  I,  no  cross  range  control  is  available  by  means  of 
attitude  adjustments  since  atmospheric  forces  are  negligible.  In  conjunction 
with  the  rotation  of  the  Earth,  it  might  seem  possible  to  provide  additional 
cross  range  by  appropriately  delaying  or  advancing  the  retro- fire  time.  For 
instance,  in  the  case  of  a  polar  orbit  with  the  landing  zone  near  the  equator, 
a  retro- fire  delay  of  20  minutes  would  appear  to  provide  5  additional  degrees 
of  cross  range.  However,  the  results  in  Sec.  7.2  are  contrary  to  such  a 
possibility  since  the  midline  of  the  footprint  closely  follows  the  vacuum 
orbital  track  over  the  Earth's  surface. 

During  Phase  II,  no  cross  range  control  was  attempted.  For 
entry  angles  of  attack  in  the  vicinity  of  ;5  degrees,  a  very  slight  increase 
in  cross  range  may  be  possible  by  trading  off  some  vertical  force  for  small 
amounts  of  side  force.  However,  the  speed  is  so  great  that  only  large  forces 
will  produce  a  turning  radius  with  any  appreciable  effect.  Thus,  since  large 
bank  angles  cannot  be  used  without  degrading  the  vehicle's  peak  temperature 
performance,  there  is  too  little  to  be  gained  from  complicating  the  reentry  with 
an  additional  parameter. 
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During  Phase  III,  about  ten  to  twenty  percent  of  the  cross  range 
may  be  realized  by  using  the  optimum  bank  angle.  Figure  27  shows  that  the 
best  value  for  the  skip  angle  of  roll  is  about  75  degrees  and  this  value  was 
adopted  for  all  subsequent  trajectories.  Another  possibility  for  producing  the 
most  cross  range  was  that  of  using  a  bank  angle  based  on  speed  or  perhaps  alti¬ 
tude  as  suggested  by  Wagner  (Ref.  6).  However,  the  result  was  usually  so  close 
to  75  degrees  that  the  additional  complexity  was  not  considered  worth  while. 

It  is  during  Phase  IV  that  most  of  the  cross  range  capability 
is  realized.  The  lower  limits  for  the  angle  of  attack  and  the  bank  angle  are 
not,  of  course,  used  until  well  after  the  beginning  of  Phase  IV.  This  is  be¬ 
cause  the  CONTAC  system  retains  control  of  the  attitude  until  so  much  energy 
has  been  dissipated  that  the  convective  heating  rate  will  continue  to  decline 
regardless  of  the  attitude  assumed.  Once  this  condition  has  been  reached,  the 
guidance  system  will  assume  full  control. 

The  basic  concept  behind  the  attempt  to  maximize  the  cross 
range  is  the  relationship  proposed  by  Wagner  (Ref.  6) .  The  problem  in  this 
investigation  was  how  to  implement  the  procedure  on  a  practical  basis.  The 
method  which  was  adopted  is  to  base  the  maximum  allowable  bank  angle  (either 
left  or  right)  on  a  variable  connected  with  the  shape  of  the  trajectory.  The 
three  variables  considered  were  the  speed,  the  altitude,  and  the  glide  angle. 

The  energy  was  also  considered  in  order  to  ccibine  both  speed  and  altitude  in 
some  manner.  However,  the  glide  angle  was  rejected  on  the  grounds  that  its 
variation  during  this  part  of  the  flight  was  an  increase  while  a  decrease  in 
the  bank  angle  v/as  obviously  required.  This  left  the  speed,  altitude,  and 
energy  of  the  vehicle  as  a  basis  upon  which  to  proceed. 

The  next  step  was  to  find  the  best  angle  of  attack  for  the 
lower  limit  during  Phase  IV.  The  value  found  to  yield  the  largest  cross  range 
was  the  angle  for  (L/D)max  which  occurs  at  22.7  degrees  for  the  vehicle  under 
study.  This  value  could  be  anticipated  due  to  the  fact  that  the  maximum  down 
range  also  results  from  the  use  of  this  angle. 

At  this  point,  a  number  of  simultaneous  tradeoffs  were  investi¬ 
gated.  No  attempt  is  /lade  to  justify  the  procedure  nor  even  the  final  result. 
However,  very  briefly,  it  involved  the  determination  of  which  power  relationship 
(square  root,  linear,  or  quadratic)  between  the  bank  angle  and  one  of  the  speed, 
the  altitude,  or  the  energy  would  yield  the  most  cross  range.  The  result  that 
was  consistantly  best  for  all  the  reentries  studied  is  the  relationship  shown 
in  Fig.  28  in  which  the  bank  angle  was  a  linear  function  of  the  speed  over  the 
interval  shown.  A  few  results  were  found  for  which  a  slightly  greater  cross 
range  was  achieved  by  using  the  square  root  of  the  speed  over  the  specified 
interval;  however,  there  was  no  way  of  predicting  when  the  result  would  be 
larger.  Therefore,  the  linear  variation  was  adopted  as  the  best  solution. 

In  general,  the  above  combination  of  the  angle  of  attack  for 
(L/D)max  and  the  relationship  shown  in  Fig.  28  for  the  bank  angle  as  the 
lower  limits  during  Phase  IV  proved  a  proper  choice  in  one  other  way.  In  most 
situations,  the  azimuth  change  at  100,000  foot  altitude  was  about  70  degrees. 

A  slight  additional  increase  until  a  speed  of  1000  feet  per  second  is  reached 
during  Fhase  V  would  yield  a  final,  azimuth  change  of  about  80  to  90  degrees. 

It  is  certainly  convenient  to  be  heading  in  the  proper  direction  to  achieve 
the  maximum  possible  additional  cross  range  if  it  is  required  during  the  last 
stages  of  the  reentry. 
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However,  no  calculations  were  made  for  Phase  V  of  the  reentry. 
But,  it  is  worth  mentioning  that  any  procedures  which  are  used  to  optimize  the 
cross  range  must  take  into  account  the  final  approach  and  landing.  Certainly, 
some  cross  range  capability  must  be  sacrificed  in  order  to  obtain  a  sensible 
landing  pattern  and  a  touchdown  with  a  margin  of  safety  and  an  allowance  for 
adverse  conditions  at  the  landing  site. 

6.3  Effect  of  Non-Standard  Initial  Conditions 


The  results  presented  in  this  section  deal  with  the  effects  pro¬ 
duced  by  a  difference  between  the  desired  and  the  actual  initial  values  of  the 
speed  and  the  glide  angle.  The  cause  of  the  difference  has  been  dealt  with  in 
some  detail  in  Section  1.8.  However,  the  investigation  should  not  be  confined 
only  to  that  limited  range  of  errors  in  attitude  misalignment  and  incorrect 
retro-fire  velocity  increment  that  have  occurred  in  recent  flights.  Therefore, 
the  results  of  large  deviations  in  the  speed  and  the  glide  angle  will  be 
examined. 


Before  any  results  are  presented,  it  should  be  noted  that  there 
are  actually  two  problems  to  consider.  The  more  important  is  whether  the 
COKTAC  system  can  cope  with  the  deviations  from  the  expected  initial  values. 

This  question  is  easily  answered.  If  the  peak  temperature  during  the  reentry 
does  not  exceed  by  more  than  a  small  value  the  peak  temperature  which  is  ex¬ 
perienced  during  an  equilibrium  glide,  then  this  test  is  successful.  If 
unsuccessful,  the  first  temperature  maximum  is  always  higher  than  the  standard 
peak  temperature .  On  the  other  hand,  when  the  test  is  successful,  the  peak  which 
occurs  during  the  flight  (whether  it  be  the  first  or  a  subsequent  maximum)  will 
be  equal  to  the  standard  peak  value. 

The  second  question  is  whether  or  not  the  vehicle  may  still  land 
at  its  original  destination.  Due  to  changes  in  the  initial  conditions,  the 
shape,  size,  and  position  of  the  footprint  is  altered.  The  alterations  are 
analyzed  and  the  basic  differences  in  the  footprint  are  given.  It  can  then  be 
determined  whether  or  not  the  footprint  still  contains  the  desired  landing  site. 

Figure  29  shows  the  results  of  deviations  in  the  initial  speed 
which  range  from  25,700  to  26,300  feet  per  second.  Wi-ch  respect  to  the  first 
problem,  it  c an  be  concluded  that  positive  errors  in  the  initial  speed  will  not 
cause  any  degredation  in  the  performance  of  the  CONTAC  system.  However, 
peculiarly  enough,  it  is  the  negative  errors  (due  to  a  lower  than  normal  initial 
speed)  which  cause  some  difficulty.  For  small  errors,  between  zero  and  -100 
feet  per  second,  the  increase  in  the  peak  temperature  is  also  small,  amounting 
to  only  50F°  for  most  points  on  the  footprint  and  only  90F°  for  the  extreme 
down  range  positions.  However,  for  errors  in  the  vicinity  of  -300  feet  per 
second,  the  performance  definitely  begins  to  decline;  peak  temperatures  are 
200F°  greater  for  most  points  on  the  footprint  and  at  the  extreme  down  range 
position,  the  peak  is  another  sixty  degrees  higher  still.  Calculations  could 
be  made  for  even  greater  errors  in  the  initial  speed,  but  a  simple  check  on  the 
initial  altitude  from  which  such  a  reentry  could  be  initiated  reveals  that  even 
for  an  error  of  -300  feet  per  second  (i.e.  an  initial  speed  at  350,000  feet  of 
25,700  feet  per  second),  the  maximum  height  from  which  the  vehicle  could  have 
descended  is  only  130  miles.  In  addition,  the  errors  in  retro-fire  alignment 
attitude  would  have  to  be  about  30  degrees  in  alignment  and  about  250  feet  per 
second  in  velocity  increment.  Such  large  errors  would  lie  outside  the  normal 
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operating  conditions.  Finally,  still  larger  errors  (on  the  negative  side  at 
least)  have  no  chance  at  all  of  occurring  unless  some  deliberate  action  to  de¬ 
crease  the  speed  is  taken  at  an  altitude  below  100  miles .  Since  this  height 
is  below  normal  orbiting  altitudes,  the  action  would  have  to  take  place  subse¬ 
quent  to  retro-fire  and  is  not  considered  a  feasible  possibility  in  “the  pro¬ 
duction  of  errors. 

With  regard  to  the  position  of  the  altered  footprint,  no  problems 
are  apparent.  The  maximum  down  range  possible  with  the  -300  feet  per  second 
error  is  still  60  degrees  of  equivalent  longitude  greedier  "than  the  minimum  down 
range  with  the  +100  feet  per  second  error.  And  with  the  latter  error,  global 
down  range  is  already  available.  Therefore,  with  the  +300  feet  per  second 
error,  the  guidance  system  will  only  have  bo  take  "the  vehicle  around  the  Earth 
one  additional  orbit  if  the  minimum  down  range  position  is  past  the  landing 
site.  With  regard  to  cross  range,  rt  is  quite  clear  that  no  problems  exist  here 
either.  The  minimum  footprint  width  is  still  almost  equal  to  the  equivalent 
latitude  span  of  the  United  States.  Also,  note  that  the  minimum  down  range 
increment  or  footprint  length  is  still  greater  than  the  equivalent  longitude 
span  of  that  country. 

Figure  30  shows  the  results  of  deviations  in  the  initial  glide 
angle  which  range  from  0.65  to  1.5  degrees.  Reference  5  indicates  that  such 
an  interval  is  far  larger  than  would  be  expected  due  to  normal  errors  in 
attitude  alignment  and  velocity  increment.  However,  Fig.  8  shows  that  this 
interval  is  to  be  expected  if  large  retro- fire  errors  are  present.  Negative 
errors  in  the  initial  glide  angle  are  easily  handled  by  the  CONTAC  system 
and  without  any  degradation  in  performance.  But  Tor  positive  errors  even  as 
small  as  0.25  degrees,  the  CONTAC  system  fails  to  prevent  the  peak  tempera¬ 
ture  from  rising  beyond  the  standard  peak  value.  An  analysis  of  the  actual 
trajectories  indicates  that  in  every  case  the  first  temperature  maximum  is 
higher  than  the  standard  peak  value  due  to  the  inability  of  the  spacecraft  to 
produce  enough  vertical  force.  As  a  result,  the  vehicle  descends  to  a  lower 
than  standard  height  in  the  atmosphere  at  that  position.  The  decrease  in 
height  results  in  a ^higher .density  without  a  compensating  decrease  in  speed. 
Therefore,  although  the  CONTAC  system  is  not  directly  responsible  for  the  tempera¬ 
ture  increase  (recall  that  it  does  not  go  into  active  operation  until  the  be¬ 
ginning  of  Phase  III  and  these  temperature  peaks  occur  at  the  end  of  Phase  II), 
nevertheless,  the  vehicle  must  be  designed  to  withstand  the  additional  rise  in 
the  peak  tempera ture  of  20oF°  if  the  initial  glide  angle  may  be  as  high  as  1.25 
degrees.  However,  it  should  be  noted  that  the  additional  rise  for  an  initial 
glide  angle  of  1.5  degrees  is  only  about  300F°  over  the  standard  peak  tempera¬ 
ture.  Consequently,  the  relationship  is  highly  non-linear. 

With  regard  to  the  position  of  the  altered  footprint  there  may 
be  a  problem  for  initial  glide  angles  greater  than  about  1.4  degrees.  Beyond 
this  value,  the  maximum  down  range  may  be  less  than  the  minimum  down  range  for 
initial  glide  angles  less  than  C.75  degrees.  Whether  the  additional  down  range 
(obtainable  from  permitting  the  vehicle  to  travel  once  around  tne  Earth  in  order 
to  overlap  the  original  down  range  position)  can  be  usefully  employed  by  the 
down  range  guidance  scheme  is  a  question  which  must  be  answered  by  simulated 
real  time  runs.  The  cross  range  is  not  a  problem;  there  is  more  than  the  usual 
value . 


A  few  additional  points  which  apply  to  both  Figs.  29  and  jO 
should  be  mentioned.  First,  there  has  been  no  attempt  to  optimize  the  upper 
and  lower  limits  for  the  attitude  when  incorrect  initial  conditions  in  the 
speed  or  glide  angle  were  present.  Second,  the  results  of  an  incorrect  initial 
value  for  the  down  range  position  has  not  been  taken  into  account.  A  one 
degree  error  in  the  initial  down  range  would  require  about  a  twelve  second 
error  in  retro-fire  time  assuming  that  the  initial  conditions  of  speed  and 
glide  angle  were  correct.  However,  the  initial  down  range  is  also  affected 
by  the  conditions  which  produce  errors  in  the  initial  reentry  values.  For 
instance,  an  increased  initial,  glide  angle  of  1.35  degrees  (which  results  in 
a  decreased  down  range  over  the  whole  footprint)  is  accompanied  by  a  decreased 
down  range  at  the  initial  position  of  eleven  degrees  of  equivalent  longitude  when 
the  retro-fire  takes  place  at  250  miles.  Consequently,  it  can  be  seen  that  the 
range  control  problem  may  be  quite  difficult  to  solve,  but  a  solution  should 
be  available  which  is  also  compatible  with  the  concept  and  operation  of  the 
CONTAC  system. 


VII.  RESULTS  AND  DISCUSSION  -  OBLATE,  ROTATING  EARTH 
7-1  Introduction 


Section  VI  presented  the  problems  associated  with  variations  in 
the  initial  reentry  conditions  of  speed  and  glide  angle.  However,  it  will 
certainly  be  possible  to  design  for  the  spacecraft  the  equipment  necessary  to 
compute  the  proper  attitude  alignment  and  retro-fire  velocity  increment  at 
the  given  orbital  position.  Therefore,  these  problems  should  not  be  considered 
to  be  very  serious.  On  the  other  hand,  the  effects  of  Earth  rotation  and 
oblateness  must  be  compensated  for  this  taskjmay  not  oe  simple  in  all  situations. 

In  order  to  understand  the  individual  effects  of  the  Earth's 
rotation  and  oblateness,  calculations  were  performed  in  which  only  one  of 
these  was  present  at  a  time.  All  the  effects  were  then  combined  and  the  over¬ 
all  effect  noted. 

7.2  Effect  of  the  Earth's  Rotation 


Figure  31  gives  the  results  of  calculations  for  a  number  of 
different  reentry  positions  and  directions  when  it  was  assumed  that  the  Earth 
was  circular,  but  rotating. 

An  analysis  of  the  peak  temperatures  shows  that  the  CONTAC 
system  encounters  a  little  difficulty  in  controlling  the  temperature  when 
rotation  and  an  initial  Westerly  heading  along  the  equator  combine  to  in¬ 
crease  the  initial  relative  speed  between  the  vehicle  and  the  atmosphere. 

On  the  other  hand,  this  effect  is  reversed  for  an  initial  Easterly  heading  and 
lower  than  normal  peak  temperatures  result  .  However,  since  the  maximum 
difference  in  each  case  amounts  to  only  -  125  F°  for  due  West  and  due  East 
equatorial  reentries,  respectively,  the  effect  of  rotation  on  temperature 
peaks  cannot  be  considered  a  major  problem.  In  most  situations,  the  effect 
will  probably  be  beneficial  since  the  loss  in  inertial  velocity  during  launch 
does  not  favour  the  establishment  of  orbits  which  travel  in  a  Westerly  di¬ 
rection  against  the  rotation  of  the .Earth. 
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With  respect  to  the  shape  of  the  footprint,  the  most  obvious 
part  of  the  result  is  the  shift  towards  negative  values  of  cross  range  without 
any  substantial  change  in  down  range.  A  comparison  should  be  made(Fig.  29).  In 
those  cases,  for  a  non- rotating  Earth,  errors  in  the  initial  speed  of  -  300 
feet  per  second  produced  large  shifts  in  final  down  range.  For  these  cases  in 
which  there  is  a  ±1500  feet  per  second  change  in  speed  relative  to  the  Earth 
due  to  an  initial  Westerly  or  Easterly  heading,  respectively,  along  the  equator, 
the  down  range  shift  averages  only  about  -  8  degrees  of  equivalent  longitude 
which  is  only  about  four  percent  of  the  total  down  range  distance.  It  may  be 
concluded,  therefore,  that  for  a  given  set  of  inertiuo.  initial  conditions,  the 
Earth's  rotation  is  not  a  significant  factor  with  respect  to  the  final  down 
range  position  of  the  vehicle.  The  physical  explanation  has  already  been  pre¬ 
sented  in  Sec.  6.2. 


However,  with  regard  to  the  cross  range,  an  important  shift 
occurs  in  the  footprint.  Fortunately,  the  extent  of  the  shift  is  quite  pre¬ 
dictable.  A  comparison  has  been  made  with  orbital  tracks  for  vacuum  tra¬ 
jectories.  In  almost  every  case  studied,  the  extent  of  the  shift  is  either 
the  same  as  or  algebraically  less  (i.e.  whether  the  actual  shift  is  positive 
or  negative)  by  up  to  1  degree  than  the  value  which  would  be  found  due  to  the 
rotation  of  the  Earth  under  the  equivalent  circular  orbit  (Fig.  3).  In  most 
cases,  the  agreement  was  better  than  §■  degree.  The  example  easiest  to  visualize 
is  the  case  for  which  the  vehicle  is  initially  heading  due  north  from  the 
equator.  For  a  circular  vacuum  orbit  through  the  initial  position,  the  Earth 
will  have  rotated  -10.8  degrees  by  the  time  the  spacecraft  reaches  the 
equator  again  at  180  degrees  down  range.  The  actual  shift  in  this  case  also 
amounts  to  -10.8  degrees.  A  check  upon  all  other  cases  yields  similar  results. 

The  observed  result  may  be  explained  as  follows.  In  symmetric 
fligl£,  with  zero  bank  and  sideslip,  neither  the  lift  nor  the  drag  have  a 
horizontal  component  perpendicular  to  V.  Thus  the  only  force  available  to 
produce  horizontal  curvature  of  the  path  (in  an  Earth- fixed  reference  frame) 
is  the  Coriolis  force.  This  is  the  same  whether  the  atmosphere  is  present  or 
not,  (except  for  the  gradual  reduction  of  V  by  the  drag).  Hence  the  horizontal 
curvature  of  the  path  is  largely  unaffected  by  the  atmosphere. 

7.3  Effect  of  the  Earth's  Oblate  Shape 

Figure  52  gives  the  results  of  calculations  for  a  number  of 
different  reentry  positions  and  directions  when  it  was  assumed  that  the  Earth's 
poles  were  depressed  to  give  the  correct  oblate  shape,  but  that  the  gravitation¬ 
al  field  was  purely  a  central  force  and  the  Earth  was  not  rotating. 

An  analysis  of  the  peak  temperatures  shows  that  the  C0NTAC  sys¬ 
tem  also  encountered  a  little  difficulty  when  the  oblate  shape  of  the  Earth 
caused  the  effective  initial  glide  angle  (see  Appendix  B)  to  be  larger  than 
normal.  For  situations  in  which  the  effective  initial  glide  angle  is  smaller 
than  normal,  the  C0NTAC  system  was  able  to  maintain  the  peak  temperature  at 
the  standard  value. 

The  situation  may  be  compared  directly  with  the  results  of  Fig. 

30  in  which  both  smaller  and  larger  than  normal  values  were  used  for  the  geo¬ 
centric  initial  glide  angle.  The  extreme  negative  difference  between  the 
effective  and  the  geocentric  initial  glide  angle  is  -0.18  degrees,  and  occurred 
for  an  initial  latitude  of  30  degrees  North  and  an  initial  heading  of  due 
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North.  The  peak  temperature  for  both  this  footprint  and  the  one  where  the 
geocentric  initial  glide  single  was  0.75  degrees  is  the  normal  peak  value.  In 
addition,  the  cross  range  is  the  same  for  each.  However,  down  range  achieved 
in  the  oblate  Earth  case  is  a  relatively  small  but  significant  amount  less  at 
both  ends  of  the  footprint  them  for  the  case  of  the  geocentric  initial  glide 
smgle  of  0.75  degrees.  This  difference  is  in  good  agreement  with  the  expected 
result . 


The  extreme  positive  difference  between  the  effective  and  the 
actual  initial  glide  smgle  is  +0.18  degrees,  and  occurred  for  an  initial 
latitude  of  60  degrees  South  smd  sm  initial  heading  of  due  North.  For  this 
reentry,  the  pesik  temperature  for  both  this  footprint  smd  the  one  where  the 
initial  glide  angle  was  1.25  degrees  is  somewhat  greater  than  the  standard 
peak  value,  but  the  difference  is,  as  expected,  less  for  the  case  of  the 
oblate  Esirth  reentry.  The  peaks  were  about  l60F°  -  190Fo  greater.  On  the 
other  hsrnd,  both  the  down  range  and  the  cross  range  achieved  under  both  situa¬ 
tions  is  substantially  the  same,  which  indicates  that  the  continued  decrease  in 
the  height  due  to  the  oblate  shape  has  a  continued  effect  on  range  while  for 
peak  temperatures,  the  effect  ceases  at  the  first  temperature  maximum  (which 
occurs  at  about  18.5  degrees  down  range  from  the  initial  position). 

The  use  of  other  geographical  locations  for  the  initial  position 
produced  similar  results.  A  slight  rise  in  the  peak  temperature  was  noted  for 
the  South-to-North  trajectory,  which  produces  an  increase  in  the  effective 
initial  glide  angle  of  only  0.07  degrees.  The  equator- to-pole  reentry,  which 
produces  a  decrease  of  0.07  degrees,  experiences  a  slight  increase  in  down 
range. 


Ii,  general,  the  results  due  to  the  shape  oblateness  of  the  Earth 
often  showed  that  much  of  the  effect  cam  be  accounted  for  by  the  difference  be¬ 
tween  the  effective  and  the  geocentric  initial  glide  angle.  In  addition, 
although  there  is  a  substantial  difference  in  many  situations  between  results 
for  the  same  effective  (over  an  oblate  Earth)  and  geocentric  (over  a  spherical 
Earth)  initial  glide  angles,  there  is  at  least  qualitative  agreement  in  all 
cases.  However,  it  should  not  be  considered  serious  if  the  correction  factors 
for  non-normal  values  of  the  effective  and  the  geocentric  initial  glide  angles 
do  not  agree.  The  guidance  system  must  reduce  to  zero  either  type  of  error. 
Since  they  are  equal  in  magnitude,  errors  due  to  non- normal  values  in  the 
actual  initial  glide  angle  will  provide  the  necessary  information  on  the  extent 
of  the  errors  which  may  be  expected  on  account  of  the  shape  oblateness  of  the 
Earth. 


There  is  another  possible  approach  to  the  problem  of  shape 
oblateness.  Rather  than  relying  on  the  CONEAC  -  guidance  system  to  correct 
for  the  errors  introduced  by  incorrect  values  in  the  effective  initial  glide 
angle,  it  should  be  possible  to  produce  the  correct  down  range  by  an  appro¬ 
priate  adjustment  of  the  actual  geocentric  initial  glide  angle.  The  adjustment 
would  have  to  take  account  of  the  geographical  location  of  the  initial  position. 
The  change  in  the  actual  initial  glide  angle  could  be  produced  by  correcting 
the  retro- fire  attituie  alignment  and  velocity  increment. 

7.4  Effect  of  the  Earth’s  Oblate  Gravity  Field 

The  results  of  including  the  oblateness  term  in  the  gravitational 
field  for  the  Earth  showed  a  negligible  difference  between  the  various  foot- 
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prints  which  result  from  different  initial  positions  and  initial  heading  di¬ 
rections.  Consequently,  it  may  be  concluded  that  even  the  C2  0  coefficient 
has  a  negligible  effect  on  the  reentry  of  a  spacecraft  and  that  all  other  terms 
may  be  completely  ignored. 

One  oinor  effect  was  found  which  tends  to  cancel  the  shape 
oblateness  contraction  in  down  range  which  resulted  when  proceeding  from  polar 
to  equatorial  regions.  For  these  cases,  the  small  reduction  in  gravitational 
force  on  the  vehicle  due  to  the  additional  distance  between  the  vehicle  and 
the  equatorial  bulge  resulted  in  a  small  increase  in  down  range  of  the  order 
of  5  degrees  at  the  minimum  down  range  position  and  10  degrees;  at  the  maximum 
down  range  position.  Since  this  increase  amounted  tc  less  than  5%  of  the  total, 
and  no  other  results  were  affected,  it  was  not  considered  significant.  Note 
that  in  the  case  of  shape  oblateness,  the  decrease  in  down  range  was  double 
this  value  when  the  initial  position  was  at  the  pole. 

7-5  Effect  of  the  Earth's  Rotation  and  Oblateness 


Figure  33  gives  the  results  of  calculations  for  a  nusber  of 
different  reentry  positions  and  directions  for  a  rotating  oblate  Earth. 

An  analysis  of  the  peak,  temperatures  shows  that  the  COlflAC 
system  again  encountered  the  same  expected  difficulty  in  controlling  the  peak 
temperature  when  rotation  and  the  initial  heading  combine  to  increase  the 
relative  speed  between  the  vehicle  and  the  atmosphere.  In  these  situations, 
an  initial  heading  which  is  due  West  together  with  an  initial  equatorial  posi¬ 
tion  results  in  an  increase  in  the  peak  temperature  of  about  10CF°  above  the 
standard  peak  temperature.  Almost  exactly  the  same  result  occurred  in  the 
case  of  a  circular,  but  rotating  Earth.  Also,  the  decrease  of  10QF°  due  to 
an  initial  Easterly  heading  is  again  present. 

In  addition,  the  expected  increase  in  the  peak  temperature  also 
occurs  when  the  effective  initial  glide  angle  is  larger  than  normal  although 
in  the  rotating  Earth  environment  the  increase  is  relatively  small.  The  worst 
situation  studied  occurred  for  an  initial  latitude  of  60  degrees  South  with  an 
initial  heading  of  due  North;  for  this  case,  the  rise  was  only  a  modest  100F° 
over  all  but  the  extreme  down  range  portion  of  the  footprint  for  which  a  rise 
of  20GF°  was  incurred. 

•  • 

It  is  conceivable  that  a  higher  temperature  increase  would  re¬ 
sult.  if  the  initial  heading  and  initial  latitude  combined  to  produce  both  a 
higher  than  normal  effective  initial  glide  angle  together  with  a  higher  than 
normal  effective  initial  speed.  Hovever,  in  order  to  produce  the  former, 
almost  due  North  initial  headings  are  required  while  the  latter  requires  almost 
due  West  values.  The  latitude  must  also  be  large  if  the  former  is  to  be  signi¬ 
ficant,  which  in  turn  reduces  the  latter.  Therefore,  it  is  deemed  quite  unlikely 
that  peak  temperature .  will  be  more  than  100F°  above  the  standard  peak  during 
reentries  for  which  the  inertial  initial  conditions  are  &  speed  of  26,000  feet 
per  second  and  a  flight  path  angle  of  1.0  degrees.  Finally,  even  if  the 
effective  initial  glide  angle  is  greater  than  1.0  degrees  and  increases  in 
the  peak  temperature  higher  than  100F°  occur,  it  would  be  possible  to  compen¬ 
sate  for  this  in  advance  of  the  reentry  by  the  appropriate  adjustment  in  the 
retro-fire  attitude  alignment  and  velocity  increment. 
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With  regard  to  the  shape  and  position  of  the  footprint,  the  com¬ 
bined  effects  of  rotation  and  oblateness  upon  both  the  down  range  and  the  cross 
range  do  produce  a  situation  which  is  difficult  to  analyze.  However,  the  re¬ 
sults  may  be  considerably  simplified  with  the  aid  of  Figs.  31  and  32  and  the 
analysis  which  applies  to  each.  Very  broadly  speaking,  the  two  effects  seem 
.to  combine  in  an  almost  linear  fashion  and  produce  the  range  distortions 
associated  with  shape  oblateness  together  with  the  cross  range  shift  due  to 
rotation.  For  example,  the  rotational  case  used  in  Sec.  7*2  i :.i  which  the  vehicle 
is  initially  heading  due  North  from  the  equator  results  in  the  identical  cross 
range  shift  of  10.8  degrees  of  equivalent  latitude  at  180  degrees  of  equivalent 
longitude  down  range;  in  addition,  this  same  case  now  manages  to  have  an  extend¬ 
ed  footprint  with  almost  the  same  maximum  down  range  as  the  oblate  shape  case. 

-Continuing  with  this  same  footprint,  the  cross  range  shift  at  360  degrees  of 
equivalent  longitude  down  range  is  21.4  degrees  as  compared  to  21.7  degrees 
for  the  circular  orbit  in  vacuum  which  goes  through  the  initial  position. 

The  above  results  are  illustrated  in  Fig.  33a.  Many  more  examples  (there  were 
40  pairs  of  point  to  choose  from  as  well  as  the  general  outline  of  each  foot¬ 
print)  -far  too  numerous  to  mention  bear  out  the  conclusions that  the  effects  of 
rotation  and  oblateneso  combine  with  very  little  interaction. 

In  conclusion,  it  should  be  noted  that  the  effects  of  both  ro¬ 
tation  and  shape  oblateness  upon  the  shape  and  size  of  the  footprint  are  large. 
The  effect  of  gravitational  oblateness  is  usually  negligible  and  is  so  completely 
masked  by  shape  oblateness  that,. it  is  never  significant.  With  regard  to  shape 
oblateness,  th.  primary  effect  is  to  change  the  down  range  position  and  span  of 
the  footprint  without  producing  any  major  change  in  the  cross  range .  On  the 
othei  hand,  rotation  produces  a  dramatic  cross  range  shift  without  any  major 
effect  upon  the  down  range.  However,  the  cross  range  shift  due  to  rotation  is 
accurately  predicted  by  the  vacuum  orbital  track  over  the  Earth's  surface  of  an 
object;  which  passes  through  the  initial  position  of  the  vehicle.  It  must  have 
the  same  azimuth  and  also  possess  the  correct  speed  (the  glide  angle  should  be 
zero)  relative  to  the  Earth's  surface  to  place  the  object  in  a  circular  orbit. 
Consequently,  although  the  cross  range  shift  is  quite  large  and  cannot  be  ignored 
. if  accurate  navigation  is  to  be  achieved,  the  guidance  system  will  be  able  to 
predict  the  extent  of  the  shift  as  a  function  of  the  down  range. 


VIII.  RESULTS  AN})  DISCUSSION  -  DENSITY  PERTURBATIONS 

8.1  Introduction 

As  noted  in  Sec.  4.6,  there  is  very  little  data  available  on 
spatial  density  variations  in  the  Earth's  atmosphere.  As  a  result,  several 
different  types  of  variations  were  proposed  and  calculations  carried  out  to 
determine  if  they  caused  any  degradation  in  the  CDNTAC  system  performance. 

Each  sub-section  gives  a  detailed  account  of  the  extent  of  the  variation  for 
which  the  CONEAC  system  could  correct,  so  that  its  limitations  are  established. 

8.2  Effect  of  Different  Basic  Models 

An  attempt  was  made  to  determine  just  how  small  the  atmospheric 
scale  height  had  to  be  before  the  CONTAC  system  would  be  unable  to  correct  for 
larger  than  normal  density  gradients.  The  details  of  tlie  density  versus 
altitude  relationship  are  shown  in  Fig.  34.  Scale  heights  between  14,000  feet 
and  26,000  feet  were  used  during  the  initial  portion  of  the  reentry.  It  can  be 
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seen  that  the  increase  in  density  at  230,000  feet  due  to  the  use  of  the  smallest 
scale  height  produced  a  density  increase  of  about  250^  which  is  an  extreme 
value  in  view  of  the  available  data  (Ref.  33)* 

The  results  were  very  satisfactory.  The  use  of  the  C0NTAC  system 
resulted  in  no  increase  in  the  peak  temperature  larger  than  10P°  above  the 
standard  peak  value  and  the  range  results  were  within  2 %  of  the  values  found 
when  using  the  U.S.  Standard  Atmosphere,  1962.  The  1959  ARDC  Standard. Atmos¬ 
phere  was  also  used.  Again,  as  expected,  there  was  no  effect  noted. 

8.3  Effect  of  Altitude- Dependent  Perturbations  in  the  Density 

In  order  to  simulate  variations  dependent  upon  the  altitude, 
information  found  in  Ref.  32  was  adopted.  The  resenblence  to  a  sine  wave 
suggested  the  curve  shown  in  Fig.  35.  The  variations  are  all  about  the  U.S. 
Standard  Atmosphere,  1962.  After  the  reentry  was  calculated  for  this  shape, 
a  number  of  others  were  investigated.  The  wave  length  was  varied  from  25,000 
feet  to  100,000  feet.  The  amplitude  was  increased  to  give  a  fractional  varia¬ 
tion  of  ±0.2  and  the  nodal  point  was  shifted  to  test  the  sensitivity  to  all 
phase  relationships. 

The  most  severe  case  used  a  wave  length  of  25,000  feet  and  an 
amplitude  of  0.2  for  the  fractional  variation.  The  above  values  together  with 
a  scale  height  of  18,000  feet  for  parts  of  the  U.S.  Standard  Atmosphere,  1962 
yield  a  density  increase  of  in  an  altitude  range  of  4,250  feet  in  the  region 
one -half  of  a  wave  length  from  the  nodal  point.  Since  Ref.  32  indicates  that 
the"density  may  vary  by  as  much  as  50  percent  when  descending  through  an 
altitude  increment  of  4,250  feet"  and  that  this  is  the  maximum  increase  to  be 
expected,  it  was  felt  that  further  calculations  were  not  necessary. 

Again,  the  results  were  completely  satisfactory  in  that  the 
C0NTAC  system  was  able  to  maintain  control  of  the  temperature.  The  down  range, 
however,  did  vary  by  as  much  as  5$  from  the  values  obtained  by  using  the  U.S. 
Standard  Atmosphere,  1962. 

8.4  Effect  of  Random  Density  Perturbations 

The  simulation  of  random  spatial  density  variations  presented 
an  unavoidable  difficulty.  Data  is  not  available  from  which  to  set  up  a.  real¬ 
istic  model  of  the  variations  which  might  be  expected.  A  little  information 
combined  with  intuitive  judgements  was  used  to  postulate  a  random  variability 
model  based  upon  the  distance  travelled. by  the  spacecraft  during  the  reentry. 

The  detailed  description  of  the  model  is  given  in  Sec.  2.6. 

The  results  may  best  be  described  by  Reference  to  Fig.  36.  The 
combination  of  variance  and  correlation  interval  above  the  hatched  zone  usually 
resulted  in  situations  with  which  the  CONIAC  system  could  not  cope.  Increases 
in  the  peak  temperature  of  about  300F°  usually  resulted  for  variances  which 
were  025  to  0.1  greater  than  the  variance  shown  acceptable  for  a  given  corre¬ 
lation  interval.  For  combinations  in  the  hatched  zone,  the  peak  temperature 
was  usually  less  than  25F°  greater  than  the  standard  peak  value;  however,  in 
isolated  examples,  an  increase  of  between  50F°  and  100F°  was  experienced.  For 
variances  which  were  0.05  to  0.1  less  than  the  variance  immediately  under  the 
hatched  region,  the  peak  temperature  was  always  equal  to  the  standard  peak  value. 
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An  exhaustive  study  could  not  be  made  due  to  the  statistical 
nature  of  the  problem.  A  complete  set  of  results  would  have  required  over  a 
thousand- individual  runs,  each  requiring  up  to  three  or  four  minutes  of  com¬ 
puter  time. 

A  large  number  of  runs  is  needed  because  a  probability  analysis 
is  required.  .For  each  combination  of  variance  and  corrflation  interval,  a  num¬ 
ber  of  different  random  density  distributions  would  have  to  be  used,  and  a 
graph  deduced  of  the  probability  of  exceeding  a  given  temperature  for  that 
combination.  Repeating  for  a  number  of  variance-correlation  pairs  would  pro¬ 
vide  the  data  for  a  Figure  like  36,  in  which  the  hatched  area  would  be  bounded 
by  curves  of  constant  probability  of  exceeding  a  set  temperature,  e.g.,  1$ 
and  10$  probability  of  exceeding  a  300F°  rise.  In  order  to  have  a  reasonable 
degree. of  confidence  in  the  results,  a  great  many  individual  runs  must  be  made 
for  each  pair  of  values  for  the  variance  and  correlation  interval.  Since  cer¬ 
tain  results  would  occur  infrequently,  the  suggested  figure  for  the  number  of 

runs  in  each  set  is  one  hundred.  It  seems  likely  that  al  least  ten  pairs  or 

sets  would  be  required.  For  the  graph  presented,  about  2§  hours  of  time  was 

used  to  make  about  50  individual  runs  for  7  pairs  of  values  for  the  varia<  cc 

and  correlation  interval. 

One  additional  factor  of  importance  must  be  noted.  In  addition 
to  the  radiation-equilibrium  temperature,  the  actual  temperature  based  on  Eq. 
(3.6)  was  calculated  for  the  nose  stagnation  point.  In  general,  the  equilibrium 
temperature  is  held  at  a  constant  value  by  the  COflTAC  system  for  a  lengtly period. 
For  this  situation,  the  difference  between  the  real  peak  temperature  (in  which 
the  skin  possesses  some  heat  capacity)  and  the  equilibrium  peak  temperature  is 
usually  less  than  5F°.  For  this  reason,  the  real  peak  temperature  has  hereto¬ 
fore  been  ignored.  However,  when  random  sj;atial  density  variations  were  en¬ 
countered  which  had  a  correlation  interval  less  than  200  miles  (and  especially 
when  this  distance  was  less  than  100  miles  -  at  20,000  feet/sec,  10-0  miles  are 
traversed  in  only  26  seconds),  there  was  often  a  large  difference  between  the 
two  peak  values.  In  one  instance  when  the  equilibrium  peak  temperature  was 
about  300F°  greater  than  the  standard  velue,  a  heat  capacity  of  2.C  BTU/ft2 
(which  might  be  provided  by  a  0.03  inch  thick  skin  of  Molybedenum  weighing  1.6 
lb/ft2)  reduced  the  excess  temperature  by  20$  to  about  240F°.  Only  a  few 
of  the  50  runs  provided  a  situation  for  which  this  result  could  occur.  Con¬ 
sequently,  no  definitive  conclusions  could  be  drawn.  But  it  will  certainly  be 
very  fruitful  to  determines  the  extent  to  which  even  a  limited  heat  capacity 
for  the  vehicle’s  skin  can  reduce  excess  peak  temperatures  caused  by  random 
spatial  density  variations  in  the  Earth's  atmosphere. 


IX.  CONCLUSIONS 

9.1  Validity  of  the  CONTAC  System  Concept 

The  results  of  this  investigation  indicate  that  the  CONTAC 
system  which  has  been  proposed  is  a  highly  effective  method  for  reducing  the 
peak  temperature  experienced  during  the  reentry  of  a  manned,  lifting  spacecraft. 
The  maximum  temperature  which  a  vehicle  experiences  during  a  fixed  atttitude 
reentry  is  lowest  for  an  equilibrium  glide.  This  maximum  is  taken  as  the 
standard  peak  value.  The  CONTAC  system  prevents  the  actual  peak  temperature 
experienced  during  the  reentry  from  exceeding  that  value  as  long  as  the  initial 
conditions  are  not  such  as  to  produce  a  value  for  the  first  temperature  maximum 


50 


that  is  higher  than  the  peak  value  for  the  equilibrium  glide.  In  addition,  if 
the  first  peak  temperature  is  higher,  then  the  CONTAC  system  is  able  to  prevent 
all  subsequent  peaks  from  exceeding  the  first . 

With  regard  to  range  control,  the  CONTAC  "yetem  is  intended  for 
use  with  any  type  of  guidance  scheme  which  uses  closed  .  >op  control.  The  avail¬ 
able  down  range  increment  for  the  vehicle  in  this  study  (which  is  similar  to  the 
M2-F2  now  currently  under  investigation)  is  in  excess  of  100  degrees  of  equivalent 
longitude.  The  combined  cross  range  for  turns  both  to  the  left  and  to  the 
right  is  greater  them  20  degrees  of  equivalent  latitude  throughout  all  of  the 
above  value  for  down  range.  Within  these  limits,  the  footprint  is  essentially 
isothermal,  i.e.  the  temperature  experienced  during  the  reentry  is  the  same  re¬ 
gardless  of  the  terminal  position  of  the  vehicle  within  the  footprint. 

9.2  Effects  of  Adverse  Conditions  on  the  Ke.utry 

The  performance  of  the  CONI  AC  system  has  been, investigated  under 
several  types  of  adverse  conditions.  The  first  type  is  errors  in  the  initial 
values  of  speed  and  glide  angle.  The  second  consists  of  those  effects  which 
can  be  predicted  in  advance,  such  as  the  Earth's  rotation  and  oblateness.  The 
third  type  consists  of  variations  from  standard  values  of  the  atmospheric  dens- 
sity  which  are  both  altitude  dependent  and  random. 

Large  errors  in  the  initial  speed  and  glide  an^le  are  considered 
unlikely  since  they  cam  be  produced  only  by  extremely  large  retro-fire  errors. 
(Past  experience  shows  that  retro- fire  errors  can  be  kept  quite  small.)  At 
the  deboost  position,  attitude  alignment  must  be  in  error  by  about  5  degrees  or 
the  velocity  increment  must  he  in  error  by  about  50  fps  before  any  significant 
changes  in  the  initial  speed  and  glide  angle  are  noted  at  the  reentry  altitude 
of  350,000  feet.  However,  if  large  retro- fire  errors  are  present,  the  following 
results  are  noted,  •  Negative  errors  in  the  initial  speed  up  to  100  feet  per 
second  (larger  errors  are  improbable)  could  cause  an  increase  in  the  peak  tem¬ 
perature  up  to  100  F°.  However,  an  increase  in  the  initial  glide  angle  of  0.5 
degrees  causes  an  increase  of  300  F°  in  the  peak  temperature. 

The  Earth’s  rotation  causes  a  maximum  increase  of  about  100  F° 
in  the  peak  temperature  when  the  Earth's  rotation  and  the  initial  heading  com¬ 
bine  to  produce  an  effective  initial  speed  1500  feet  per  second  above  the 
standard  value  of  26,000  feet  per  second.  When  the  Earth’s  oblateness  and  the 
idLhis.i.  i^ding  and  latitude  combine  to  produce  an  effective  initial  glide  angle 
of  0,i8  degress  greater  them  the  standard  value  of  1.0  degree,  the  maximum  in¬ 
crease  in  the  peak  temperature  is  again  about  100  F°.  However,  the  rotational 
condition  is  unlikely  to  occur  in  practice  since  satellite  orbits  areuusually 
in  the  same  ser.je  as  Earth's  rotation,  and  the  oblateness  effect  may  be 
compensated  for  in  advance  of  the  reentry  by  an  appropriate  adjustment  of  the 
retro-fire  attitude  alignment  and  velocity  increment. 

Atmospheric  density  perturbations  based  on  altitude  did  not,  with¬ 
in  reasonable  limits,  produce  any  situations  for  which  the  CONIAC  system  could 
not  correct.  The  maximum  likely  vertical  density  gradients  had  no  effect  on 
the  CONTAC  system  performance.  Ontheopther  hand,  random  spatial  (roughly 
horizontal)  density  variations  produced  a  very  different  result.  The  investi¬ 
gation  described  a  basis  for  judging  the  adequacy  of  the  CONTAC  system  when 
sufficient  information  about  random  density  variations  in  the  Earth's  atmos¬ 
phere  become  available.  Insufficient  data  could  be  found  for  the  altitudes  of 
concern,  namely  between  220,000  feet  to  290,000  feet,  in  order  to  arrive  at  a 
firm  conclusion  at  this  time. 


51 


REFERENCES 


1.  Wingrove,  R.C. 


2.  Bryson,  A.E. 
Denham,  W.F. 
Carrol,  F.J. 
Mikami ,Kinya 


A  Survey  of  Atmospheric  Reentry  Guidance  and  Control 
Methods,  IAS  31«t  Annual  Meeting,  January  21-23,  1963, 
New  York,  N.Y.,  IAS  Paper  No.  63-86. 

Determination  of  Lift  or  Drag  Programs  to  Minimize 
Re-Entry  Heating,  Journal  of  the  Aerospace  Sciences, 
Vol.  29,  No.  4,  April,  i.962. 


3.  Levinsky,  E.  S.  Application  of  Inequality  Constraints  to  Variational 

Problems  of  Lifting  Re-Entry,  IAS  29th  Annual  Meeting, 
January  23  -  25,  196l,.New  York,  New  York,  IAS  Paper 
No.  61-21. 


4.  Stalony-Dobrzan8ki,pJ.j,n^.ry  g^^ance  Control  Using  Temperature  Rate 

Flight  Control  System,  AIAA  2nd  Aerospace  Sciences 
Meeting,  January  25-27,  1965,  New  York,  N.Y., 

AIAA  Paper  No.  65-47. 


5.  Collins,  R.  R. 
Dotts,  H.W. 
Hoyler,  W.F. 
Hecht,  K.F. 

6.  Wagner,  W.E. 


7.  Smith,  R.H. 
Menard,  J.A. 


8.  Gaines,  L.M. 
Suroer,  T.E. 


9.  Thompson,  M.O. 
Weil,  J. 
Holleman,  E.C. 


10.  Etkin,  B. 


Guidance,  Control,  and  Propulsion  Systems,  Gemini  Mid¬ 
program  Conference  Including  Experimental  Results, 

February  23-25,  1966,  Manned  Spacecraft  Center,.  Houston, 
Texas,  NASA  S P-121.  '  _ 

Roll  Modulation  for  Maximum  Re-Entry  Lateral.  Range, 

Journal  of  Spacecraft,  Vol’.  2,  No.  5,  September-October 
1965. 

Supercircular  Entry  and  Recovery  with  Maneuverable  Manned 
Vehicles,  National  IAS,  ARS  Joint  Meeting,  June  13-16, 

.  1961,  Los  Angeles,-  California,  IAS  Paper  No.6l-ll4-l8o8. 

Prediction  of  Optimum  Approach  and  Landing  Techniques  for 
Manned  Re-Entry  Gliders,  National  IAS,  ARS  Joint  Meeting, 

June  13-16,  1961,  Los.  Angeles,  Calif.  IAS  Paper,-  No. 

6I-II&I809. 

An  Assessment  of  Lifting  Reentry  Flight  Control  Require¬ 
ments  during  Abort,  Terminal  Glide,  and  Approach  and 
Landing  Situations,  AGARD  Specialists'  Meeting  on 
Stability  and  Control,  Sept.  20-23,  1966,  Cambridge,  England. 

Longitudinal.  Dynamics  of  a  Lifting  Vehicle  in  a  Circular 
Orbit,  UTIAS  Report  No.  65,  AF06R  TN  63-191,  Feb.  i960. 


11.  Drummond,  A.M 


Performance  and  Stability  of  Hypervelocity  Aircraft  on 
a  Minor  Circle,  UTIAS  Report  to  bq.. published. 


12.  Technology- Week  Including  Missiles  and  Rockets,  American  Aviation 

Publications,.  Inc.,  Vol.  19>  No.  20,  July  25,  1966. 

13.  Geiger,  R.E.  Experimental  Lift  and  Drag  of  a  Series  of  Glide  Con¬ 

figurations  at  Mach  numbers  12.4  and  17.5,  IAS  National 
Summer  Meeting,  June  28-July  1,  i960,  Los  Angeles, 
California,  IAS  Paper  No.  60-92. 

14.  .-Feltermaa,  JD.E.  An  Analysis  of  the  Delta-Wing  Hypersonic  Stability  and 

lteal,.  L,  Jr.  Control  Behaviour  at  Angles  of  Attack  Between  30°  and 
90°,  NASA  TN  D-1602,  March  1963. 

15.  Wells,  W.R.  Tables  of  Aerodynamic  Coefficients  Obtained  from  Developed 

Armstrong,,  W.O.  Newtonian  Expressions  for  Complete  and  Partial  Conic  and 

Spheric  Bodies  at  Combined  Angles  of  Attack  and  Sideslip 
with  Some  Comparisons  with  Hypersonic  Experimental  Data, 
NASA  TR  R-127,  1962. 

16.  •  Smith,  R.Hc  Supercircular  Entry  and  Recovery  with  Maneuverable  Manned 

Menard,  J.A.  Space  Vehicles,  LAS-ARS  Joint  Meeting,  June  13-16,  1961, 

Los  Angeles,  Calif.,  IAS  Paper  No.  6l-ll4-l8o8. 

17.  Chapman,  D.R.  An  Analysis  of  the  Corridor  and  Guidance  Requirements 

for  Supercircular  Entry  into  Planetary  Atmospheres, 

NASA  TR  R-55,  I960. 

18.  Mugler,  J.P.  Jr.  Static  Longitudinal  Aerodynamic  Characteristics  at  Tran- 

Olstad,  W.B.  sonic  Speeds  of  a  Blunted  Right  Triangular  Pyramidal 

Lifting  Reentry  Configuration  for  Angles  of  Attack  up 
to  110°,  NASA  TN-  D-797,  June  1961. 

19.  Fournier,  P.G.  Wind-Tunnel  Investigation  at  High  Subsonic  Speed  of  the 

Static  Longitudinal  Stability  Characteristics  of  a  Winged 
Reentry  Vehicle  Having  a  Large  Negatively  Deflected  Flap- 
Type  Control  Surf ace, NASA  TN  D-2030,  October,  1963. 

20.  -Ware,  G.  M.  Low-Subsonic-Specd  Static  Stability  of  Right-Triangular - 

Pyramid  and  Half-Cone  Lifting  Reentry  Configurations, 

NASA  TN  D-646,. February  1961. 

21.  -Crabtree,  L.F.  Estimation  of  Heat  Transfer  to  Flat  Plates,  Cones  and 
Doomett,  R.L.  Blunt  Bodies,  Royal  Aircraft  Establishment  Technical 
Woodley,  J.G.  Report  No.  65137 >  July  1965* 

22.  Leo  Rute  A  Study  of  Aerodynamic  Effects  of  Isothermal  and  Tempera¬ 

ture  Gradient  Atmospheres  on  Reentry  Trajectories, 
Polytechnic ■ Institute  of  Brooklyn  Department  of  Aero¬ 
space  Engineering  and  Applied  Mechanics,  PIBAL  Report 
No.  738,  March  1962. 

23.  Parvin,  R.H.  The  Earth  ana  Inertial  Space,  Part  I  -  Motions  of  the 

Earth,  Aerospace  Engineering,  Vol.  18,  No.  4,  April  1959* 


53 


24..  .Koz°,i,  Yoshihide  Numerical  Results  on  the  Gravitational  Potential  of  the 

Earth  from  Orbits,  Proceedings  of  the  First  International 
Symposium  on  the  Use  of  Artificial  Satellites  for 
Geodesy',  April  26  -  28,  1962,  Washington,  D.C. 

25.  Kaula,  W.M.  Tesseral  Harmonics  of  the  Gravitational  Field  and 

Geodetic  Datum  Shifts  Derived  from  Camera  Observa¬ 
tions  of  Satellites,  RASA  TN  D-1848,  June  1963. 

26.  Kaula,  W.M.  Improved  Geodetic  Results  from  Camera  Observations  of 

Satellites,  Journal  of  Geophysical  Research,  Yol.  68, 

No.  18,  Sept.  15,  1963. 

27.  Wagner,  C.A.  The  Drift  of  a  24-Hour  Equatorial.  Satellite  Due  to  an 

Earth  Gravity  Field  Through  4th  Order,  NASA  TN  D-2103, 
February  1964. 

28.  Wagner,  C.A.  The  Equatorial  Ellipticity  of  the  Earth  as  Seen  from 

Four  Months  of  Syncoo  II  Drift  Over  the  Western  Pacific, 
NASA  TN  D-3313,  July  1966. 

29.  Wagner,  C.A.  The  Equatorial  Ellipticity  of  the  Earth  from  Two 

Months  of  Syncoo  II  Drift  Over  the  Central  Pacific, 

NASA  TN  D-3315,  July  1966. 

30.  Baker,  R.M.  Jr.  Three-Dimensional  Drag  Perturbation  Technique,  APS 

Journal,  Yol.  30*  Ho.  8,  August  i960,  AFOSR  TN  59-767, 
July  1,  1959. 

31.  Staff  U.S.  Standard  Atmosphere,  1962. 

32.  Cole,  A.  E.  Horizontal  and  Vertical  Distribution  of  Atmospheric 

Kantor,  A.J.  Density  Up  to  90  Km,  Air  Force  Surveys  in  Geophysics, 

No.  157,  AFRCRL-64-483,  June  1964. 

33.  Cole,  A.E.  Model  Atmospheres  (Chapter  2,  Handbook  of  Geophysics 

Court,  A.  and  Space  Environment),  Air  Force  Cambridge  Research 

Kantor,  A.J.  Laboratories,  S.G.- Hancom  Field,  Bedford,  Massachusetts . 

34.  Etkin,  B.  Private  Communications,  1966. 

35.  Sommer,  S.C.  Point  Return  from  a  Lunar  Mission  for  a  Vehicle  that 

Short,  B.J.  Maneuvers  Within  the  Earth's  Atmosphere,  NASA  TN  D-1142, 

Novembex*  1961. 

36.  Nielsen,  J.N.  Three-Dimensional  Orbits  pf  Earth  Satellites,.  Including 

Goodwin,  F.K.  Effects  of  Earth  Oblateness  and  Atmospheric  Rotation, 

Mersam,  W.A.  NASA  MEMO  12-4-58A, . December  1958. 

37-  Fine,  J.H.  Stability  of  Flight  Paths  of  Lifting  Vehicles  During 

Entry  into  Planetary  Atmospheres,  ITT IAS  Technical  Note 
No.  48,  July  1961,  AFOSR  1488. 


54 


38.  Brown,  R.C. 
Bciille,  R.V, 
Giffin,  (r.Di 

39 »  Rtkin,  B. 


1K).  .Etkin,  B. 


S ix-Degree -of -Freedom  Flight  Path  Study  Generalized 
Computer  Program,  Bart  I,  Problem  Formulation,  WADD 
Technical  Report  Bo.  60-781,  Part  I,  May  1961. 

Dynamics  of  Flight,  J,  Wiley  &  Sons,  Inc.,  New  York, 

1  OCA  ' 


Dynamics  of  Flight,  Appendix  A-6,  -  Derivative  of  a 
Vector. 


PROTECTIVE  LIFE  -  HOURS 


1200  1600  Z'JOO  2400  2800  3200 
TEMPERATURE  ~  *F 


nwsiPTs 


FIGURE  1  ENDURANCE  TIME  vs  TEMPERATURE 
FOR  COATED  REFRACTORY  METALS 

(From  ‘Refractory  Metals  and  their  Protection’, 
Thomas  D.  Cooper  and  Oscar  0.  Sharp,  Aerospace 
Englneerini,  Vol.  22,  No.  1,  January  1963.) 


COMPUTER  PREDICTED 


Rodnoy  C.  Hinpovo,  IAS  31st  Annual  Weetinc,  January  21-23, 1963, 
Now  York,  N.Y.,  IAS  Papor  No.  63-86. 


Hl-10  Artiitic  liptHiMiH 
Mmir,  HtwHwnw.  California 
Jaxwry,  I  Hi 


n«««<  typical  VCHICLC  CONFIGURATION 


MAXIMUM  LATITUDE  OF 
ORBITAL  ENVELOPE  EQUAL 


SINGLE  ORBITAL 
TRACK  OVER 


ORBITAL 

INCLINATION 

ANGLE 

RELATIVE  TO 
THE  EARTH 


EQUATOR 


ORBITAL  ENVELOPE 
Earth* 8  rotation, 
under  the  envelo- 
least  once  even' 


(By  reason  of  the 
complete  surface 
.  :i  covered  at 
24  hours) 


FIGURE  5  ORBITAL  ENVELOPE 


FIGURE  6  THE  FIVE  REENTRY  PHASES 
(Hot  To  Scale) 


FIGURE  7  PHASE  I  -  PRE-ATMOSPHERIC  FLIGHT 


INITIAL  POSITION  POR  START  OP  CALCULATION 
fbese  II  starts  at  390,000  faat  with  an 
inertial  speed  of  24,000  ft/aac  and  an 
inertial  glide  angle  of  1.00  degrees. 


TYPICAL  RESULTS  AT  FIRST 


ANCLE  OP 

SPEED 

ALTITUDE 

DOWN  RANGE 

CLIDE  ANCLE 

ATTACK 

(fpe) 

(faat) 

(degrees) 

(degrees) 

39® 

29S70 

271*00 

11.92 

0.022 

*0« 

29S20 

273900 

IS. 02 

0.025 

*90 

29770 

279100 

17.69 

0.031 

990 

29630 

275*00 

17.59 

0.0*7 

*9® 

29630 

272*00 

IS. 2* 

0.07* 

/jUyli  of  attack  coaatant  at  a 
.  value  between  39°  and  65>  o, 
nodulated  between  tbaao  llnlta 
to  control  tba  Sown  rang a. 


ROTS:  Orbital  velocity  at  279,000 
foot  la  29,770  ft/aac. 


■  FIRST  Vax  SICKALS  END  OP  PHASE  II 


FIGURE  9  PHASE  II  -  OYIXALL  DOW  SAKE  CONTROL 


»  G0VTAC  SYSTEM  W  PASSIVE  OPERATION 
(Attitude  Controlled  by  Culdance  Syataa) 

-  OOKTAC  SYSTEM  III  ACTIVE  OPERATION 


■  skip  PotnoN  op  ns  reentry 

(Nay  ba  extended  or  reduced  to  provide  aadlua  down  fangs  control  by 
tba  appreciate  choice  of  tba  attitude.  Sene  croaa  range  control  la 
aleo  available  during  tbla  portion  by  lining  large  bank  anglea.) 


.SECOND  da,,  AGAIN  MINCS  C0NTAC 
SYSTEM  INTO  ACTIVE  OPERATION 


n 

MIU 

r~ 

SI 

ICINNIK 

.  SIGNALS 
OP  PRASE  IK 


Skip  angle  of  attack  and  aklp  bank  aagla 
•  are  bold  conatant  or  nodulated  by  the 

ayatan  to  provide  range  control. 


OONTAC  ayaton  no molly  oporatoa 
-for  only  a  few  eeconde  and  brings 
attitude  to  tba  aklp  position. 


PEAT  TEMPERATURE  USUALLY  OCCURS 
AT  A' OUT  249,000  PEET  AM)  SICNALS 
THE  EhD  or  PHASE  III  I 


OPTIMUM 

-ATTITUDE 

RZACHED 


C3NTAC  SYSTEM  SL0NLY  SRINCS  1NE  ATTITUDE 
TO  TEE  OPTINM  POSITION  (99°  for  the  angle 
of  attack  and  0°  for  tba  baok  angle)  MILE 
KEEPING  TEE  TS7MRATURI  CONSTANT 


length  of  oklp  portion  la  dictated  prlaarlly  by  spend  inherited  Iron 
phono  II  of  tba  reentry.  When  c Oakland  vltfc  a  low  angle  of  attack, 

•  Rp,—  slightly  higher  than  orbital  velocity  produces  extended  skips 
up  to  120  degrees  of  ogulvnlant  longitude.  Mien  large  angles  of  attack 
and  largo  bank  anglea  are  c Oakland  with  a  spaed  slightly  lose  then 
orbital  velocity,  the  skip  portion  of  pheee  III  can  be  elialoeted 
cnpliuli  and  the  total  deem  range  Increase!  during  phase  III  can 
be  reduced  to  about  10  degrees  of  soul valent  longitude. 


Jptinun  attitude 
held  constant 


PICUIE  10  PRASE  III  -  SQUILIVRIUN  CLIDE  ATTAINMENT 


FIGURE  U«  DETAILED  HISTORY  OF  ANGLE  OF  ATTACK  FOR  A  PORTIOR  OF  A  COKTAC  CONTROLLED 


ALERT  CALCULATION  SECTION 
TO  TRY  ACAIN  WITH  A  SMALLER 
SIZE  PATH  STEP  AND  SET  UP 
PATH  STEP  TO  BE  USED 


ALERT  CALCULATION 
SECTION  TO  COMPUTE 
REN  VALUES  OF 
(M,t«)  *nd  (L,tL) 


Kin  1mm 

Ttftt 


is  (T-TU.l.ln) 

*  *d«»db»nd 


ftiHtlmi 


fluiUva.  I, 


FIGURE  12  CONTAC  STSTM  PROGRAM  FLOW  DIAGRAM 


•It  ORICIMAL  POINTS  IN  DICITAL  CALCULATION 


Td«*dt>»nd 


TEMPERATURE 


FIRST  POINT  AT  WHICH, 
>Td««db*nd  " 


DISTANCE  ALONG  FLIGHT  PATH 


FIGURE  13  TEMPERATURE  v*  DISTANCE  ALONG  FLIGHT  PATH 
DIGITAL  POINTS  AND  Tdt,db,nd  DIFFICULTt 


ACTUAL  POSITION  AT  WHICH 


FICUM  1A  TEMPERATURE  v«  DISTANCE  ALONG  ELICKT  PATH 
DICITAL  POINTS  AND  Td.*<ib»iui  CIPPICULTT 


VEHICLE 

DYNAMICS 


FIGURE  15  ONBOARD  DICITAL  CONTAC  SYSTM  PLOW  DIACRAM 


HORIZ 


FI CUKE  20  FLIGHT  AXES  COORDINATE  REFERENCE  SYSTEM 


00  axle  la  pointing  out  of  page  fro*  Indicated  poaltlon 


(v)  axla  1*  pointing  into  page  at  lndicatad  poaltlon. 


I 


i 


e 

o 


m  a 

o  o 

O. 


*8 

*8  • 
8  *± 

m 

u 

8  -8 
M  8 


M  8 
« 

e.  « 


o  e. 

u  O 
P  *i 
O  8 

$  - 

•H  8 
4J  M 

8  W 
«H  e 
O  *H 

0.  o 


IM 

iH 


i 


FICURE  24  CONVERSION  BETWEEN  SECOND  INTERMEDIATE  REFERENCE  SYSTEM  {n,«,v)  FIGURE  23  CONVERSION  BETWEEN  FIRST  INTERMEDIATE  REFERENCE  S.oTEN  {I* 

AND  THIRD  INTERMEDIATE  REFERENCE  SYSTEM  (h.J.v)  "TO  SECOND  INTERMEDIATE  REFERENCE  SYSTEM  <n,«,v} 


1.  {j>  axis  Is  pointing  out  of  page  from  indicated  position 

LOCAL 


FIGURE  25  CONVERSION  BETWEEN  THIRD  INTERMEDIATE  REFERENCE  SYSTEM  {h , j ,7} 
AND  FLIGHT  AXES  REFERENCE  SYSTEM 


3300 


DOW  RAMCC  ACKIIVtD  AT  100,000  PUT  (4*tr«««) 


rtCUM  2*  VARIATION  OF  LOVE*  LIMIT  FOR  TRI  BARK  ANCLE  vs 
SPEED  DURING  FNASE  IV  OF  THE  RIENIRT 


FWURE  »  CROSS  NAME  n  IXMK  RAN8E  FOR  VARIOUS  RMTIAL  SPEEDS 
WER  A  CIRCULAR,  NON-ROTATWQ  EARTH 


FIGURE  31e  CROC  RANGE  n  DOWN  RANGE  FOR  TiO  INITIAL  POSITIONS 
OVER  A  CIRCULAR,  ROTATING  EARTH 


FItfUW  36  OPERATIONAL  REGION  OF  CONTAC  SYST*M  FIGURE  35  FRACTIONAL  VARIATION  OF  DENSITY  WITH  ALTITUDE 


APiENDIX  A 


The  Equilibrium  Glide  Trajectory 

An  equilibrium  glide  trajectory  is  &  coplan&r  reentry  during 
which  the  lift  is  always  equal  to  the  net  difference  between  the  gravitational 
and  the  centrifugal  forces  on  the  vehicle,  Since  drag  is  not  compensated  by 
thrust,  the  speed  is  gradually  reduced  thereby  requiring  an  increase  in  lift 
to  offset  the  loss  in  centrifugal  force.  Since  the  angle  of  attack  is  fixed, 
this  requirement  is  met  by  continually  descending  to  a  lower  altitude,  but  a 
controlled  rate  so  that  the  lift  never  becomes  too  large.  As  a  consequence, 
the  reentry  is  characterized  by  the  absence  of  any  phugoid  or  long  period 
oscillations  involving  the  exchange  of  speed  for  altitude  and  vice  versa. 


For  a  simple  coplanar  reentry  over  a  circular  non-rotating  Earth, 
the  equations  of  motion  are  given  by 


m  V  +  D  -  mgsin7  =  0 


(A.l) 


m  V  7  +  L  +  m 


C0S7  =  0 


(A.2) 


For  an  equilibrium  glide,  all  the  important  effects  of  interest 
occur  between  100,000  feet  and  350,000  feet.  Consequently,  the  values  of  r  and 
g  will  vary  by  less  than  -1^  from  the  average  values  over  the  above  altitude 
range.  Therefore,  it  is  reasonable  to  make  r  and  g  constant. 


During  an  equilibrium  glide,  it  is  assumed  that  the  m  V  y  term 
in  Eq.  (A.2)  is  negligible  compared  to  the  lift.  Note  that  this  term  cannot 
always  be  equal  to  zero  since  this  would  require  a  constant  glide  angle. 


Finally,  the  analysis  win  be  restricted  to  the  region  for  which 
7  «  1.  (For  the  vehicle  used  in  this  investigation,  7  <  5°  until  the  speed 
drops  to  4300  feet/sec  at  an  altitude  of  140,000  feet.)  .Also,  the  lift  and 
drag  are  expressed  in  the  usual  manner,  in^terms  of  Cp  and  Cp.  Equations  (A.l) 
and  (A.2)  become  ' 

?  +  Cj  £pV2  S/m  -  g7  =  0  (A.3) 

CL  Jpvfym  =  (g  -  ^  )  (A-4) 

Equation  (A. 4)  is  then  rewritten  to  give 


2  m  1  -  s2 
p  -  CL  S  r  s2 


(A.5) 


where  the  speed  ratio  is 


s  =  V/^rg  (A. 6) 

Equation  (3.7)  can  be  used  to  find  the  temperature  as  a  function  of  speed. 

The  result  is 


T4  =  Ks2,15  (l-s?)^ 


(A.7) 


where 


(A.8) 


K  = 


17,600  fr/rg) 


3.15 


€  CT 


N^n(  26 , 000 )  3  * ' 15n6 . 002378 


2m 


at 


It  is  straightforward  to  show  that  the  peak  temperature  occurs 


s*  =  *hjl  =  O.8165 


(A. 9) 


The  use  of  GM  and  R0  from  equations  (4.7)  and  (4.8),  respectively, 
and  using  g  and  r  for  the  average  height  of  250,000  feet  gives  the  values  at 
the  position  of  the  peak  temperature  as 


/  *  V-/8 

**  ■ 3035  Gw)  * 

(A.10) 

V*  =  21,050  feet/sec 

(A. 11) 

p*  =  4.722  x  10"8  (“cV"-)  slu«s/ft3 

(A. 12) 

Equations  (A. 10)  to  (A. 12)  show  that  the  peak  temperature 
always  occurs  at  the  same  speed.  Also,  the  height,  y*  (actually  the  density) 
is  dependent  only  on  the  lift  coefficient  and  the  wing  loading.  From  Eq. 

(A. 10),  the  smallest  value  for  T*  is  obtained  for  the  maximum  value  of  vertical 
force,  i.e.  of  C^S/m.  Consequently,  the  smallest  T*  is  obtained  by  holding  the 
vehicle  at  the  attitude  for  CLmax«  The  expression  for  CL  in  Eq.  (3.3)  is  a 
maximum  for  an  angle  of  attack  of  55  degrees,  which  is,  therefore,  the  optimum 
angle. 


At  the  optimum  angle  of  55  degrees,  Cip/m  is  0.5417  sq.ft/slug 
for  the  subject  vehicle.  Substitution  of  this  value  together  with  e«/Rn  =1.0 
into  Eqs.  (A.IO)  and  (A. 12)  yields  the  values  in  Table  A.l. 

An  equilibrium  glide  can  be  realized  very  closely  in  practice  if 
the  proper  initial  conditions  are  achieved  at  the  initial  position.  Shown  on 
Figure  26  are  the  values  of  speed  and  glide  angle  which  are  required  at  350,000 
feet  altitude  above  a  circular,  non-rotating  Earth.  Initially,  the  speed  and 
glide  angle  are  adjusted  to  produce  a  negligible  value  for  mV>  in  Eq.  (A. 2). 
it  is  found  that  the  required  condition  of  a  negligible  mV7  (relative  to  L) 
then  prevails  throughout  the  flight  without  any  additional  control.  The  values 
of  T*,  V*,  p*  are  also  shown  in  Table  A.l  for  the  trajectory  computed  from 
the  exact  equations.  The  valtie  of  y*  is  also  included  for  this  case  and  may 
be  comp  i  to  the  value  of  y*(p*)  which  results  from  Eq.  (A. 12)  when  the  U.S. 
Standard  cmosphere,  1962  is  the  curve  of  density  versus  altitude.  The  same 
curve  was  used  for  the  numerically  computed  trajectory. 

The  agreement  between  the  two  sets  of  data  shows  that  equations 
(A. 10)  to  (A. 12)  are  quite  accurate.  Consequently,  the  standard  peak  temperature, 
T*,  as  given  by  Eq.  (A. 10)  is  shown  to  be  independent  of  (L/D).  Also,  the 
smallest  value  of  T*  is  achieved  for  the  smallest  value  of  (m/e2R  CjrS) .  In 
addition,  because  of  the  1/8  power,  £he  gains  to  be  realized  fromneven  a  50$ 
reduction  in  this  parameter  are  not  large;  when  T*  equals  3000°F,  the  reduction 
in  T*  is  only  320F0  or  only  a  little  more  than  10$. 


A2 


TABLE  A1 


Flight  Conditions  at  the  Position  of  Peak  Temperature  During  an  Equilibrium  Glide 


Value 

Approximate  Analytical 
Results* 

Numerical  Integration  of 
Exact  Equations* 

28l7°F 

28l0°F 

V* 

21,050  feet/sec 

20,800  feet/sec 

P* 

8.72  x  10“ 8s lugs/ft 3 

9.27  x  10“8  slugs/ft3 

y* 

245,500  feet*** 

243,800  feet*** 

+  equations  (A.IO)  to  (A. 12) 

T  initial  conditions  of  Big.  26 
***  corresponding  to  p*  for  the  U.S.  Standard  Atmosphere,  1962 


APPENDIX  B 


The  Effective  Initial  Glide  Angle 

The  definition  of  the  effective  initial  glide  angle  over  an 
oblate  Earth  is  quite  straightforward.  It  simply  takes  account  of  the  effective 
change  in  height  of  the  vehicle  which  occurs  during  Phase  II  of  the  reentry 
due  to  the  change  in  the  Earth's  radius.  This  portion  of  the  flight  has  been 
chosen  since  it  is  during  this  period  that  both  the  peak  temperature  and  the 
overall  down  range  are  determined.  As  shown  in  Pig.  37,  the  change  in  the 
Earth's  radius  may  be  used  to  give 

sinA?  =  (B.l) 

where  Ay  is  defined  in  Fig.  37. 

For  this  investigation,  both  y  and  A7  are  always  very  small  during  Phase  II. 
Consequently,  the  small  angle  approximation  may  beiused  to- give:  v.  c  ■  V  v  .Ue 


7eff 


+  & 
'  e  s 


(B.2) 


In  practice,  Ay  is  always  very  easy  to  find.  It  is  the  difference 
in  the  height  over  the  Earth's  surface  when  the  same  change  that  takes  place 
during  Phase  II  in  geocentric  latitude  also  occurs  at  the  average  geocentric 
radius  which  is  used  during  Phase  II.  The  value  of  s  is  the  flight  distance 
travelled  during  this  portion  of  the  reentry. 
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40ft  •  0.9 

tot*  •  9.0 

aotv  •  o.o 

CO  TO  IT* 

im  toon  coot  >  i  m.  in,  m 

IM  *90f  I  IM.  lit.  ITT 

MTi  •  ((  IIRICIMII  •  COIIOIII  •  IhIHttll 

i  -  colUHMil  •  coiiohtai  n  • 

9  (Ml  »  COttMtMl  /  ft 

I  -  II 91 91 1 1  "  C01I0HTAI  •  T9  /  V 

•9t9  •  ti9t9iti*ti9iouTai/eoti«*M*t  •  iimuiimm/ii 

I  »  II  IIMH9UI  /  CHICMMI  »  C01IMI  I  t  COlloUTAI 

1  *  llRIHltll  II  •  TM  /  V 

«0fV  •  II  COItUMII  •  CHI91II  •  llOIOflTAI 
I  •  1I9ICA9MI  #  C01IMLT4I  ||  • 

>  mj  ♦  a  •  cmionfai  •  I  •  i.«  I 


If*  MM*D  •  0.9 

Mlrt  •  0.9 
Kin  • 

*9*  TV  •  0.0 


9*0  •  9.9 

9  TO  •  0.* 

ITT  |M  KMT  »  4.0  I  in.  110.  191 
•  14  MMT  •  -  1.9 

IT*  091090  •  9.0 


It*  IM  TM9TT  I  Ilf.  119,  Ilf 
Ilf  rfM  -  1*99.9  •  TOfofv  t  T.O 

CHI  M«r*l 
90  TO  ITt 

II*  9*4*09  •  10*0900.9 


*4*0  III  09  OAT  A  TH  Iff  1  MMKM  CM  *9  MT«9f  •  M9H  0*9TM  CMT»K 
III  1*1  )  I  |99.  14).  Iff 

I9»  CHI  M4T9C 


MUTTIM  ICIMl  *.*  9*«(  9 


tUOAOVf  |Nf  ICtfll  9.9 


m 

149 

IT 

It* 

m 


*f  A*  I*  09  TUIMKI  M*  MWIWIHI  M9*  MCI *1*9 T 
IM  A9MI  9imt.  1.0  I  T  199.  191.  199 
MM  14.111  VI*tf*»Tt9fT*.9*««|T.9*f*(f.*M»ir.T*|*|f 
909OATI  Wt.9.  MI9.4  I 

9*91 1C 10*94 1  |  19*.  19*.  19* 

. .  »  144.  19*.  149 

tri*if* 


Itt  AO  • 

*n  *0  190 
19*  A*  • 

90  TO  149 
199  CHI  9MIM4 
*0  TO  199 
I9A  *0  a  91 

14* 


■MM9  •  f  |*|  T9 


TUMT* 


«*f*IT  •  9*909* 

9*10  *  9019 IT  •  9  99**4 

991  TM  •  MI*:T  •  9*9*** 

TMTM  •  T0I9IT  •  MfltH 

A*Tf  IT  •  It  Ml  A*ff4T.  -  MITA  I 

MTItT  •  f.O 

o**tr  *  i«mii 
40MCT  •  0,0 

A9T4ST  •  tlUH  AMAH Mf MU-9.0»m«Tl.|  .91.  -  4 
M9T09  •  *199*1 

A9M9CT  •  *.* 


1*1  191  *l«|tl  I  |*1.  1*9.  194 

1*9  V*  •  *99*  T I  9**V9T  *  4119  /  M  I 

*99 

CA  TO  1*9 

1*9  V*  •  *99*TI  4M  /  9*  I 

V|*  •  M  *  V* 

M  TO  |»9 

14*  V*  •  Vf*ff* 

Tl*  *  M  9  V* 

1*4  191  *9*T  I  1*4.  I9f.  IM 

1*4  VM  •  M  *  *CMI  I  MOM  I 

*•  •  f*  •  *9191  IM9M  I 

Ml  •  M  9  M 

T9  •  91  *  9C0II  9«|*  I 

Tf  •  9»  •  *9t*f  941*  I  -  M  9  9*  9  *CM I  MITA0  • 

r»9  »  W  9  N  ♦  T9  •  Tf 

CAM**  •  *419*11  94TMI  M  /  449*  TIMM  I.  MI9IT  | 

T*  •  *99911  9}  •  Wt  I 

T9*  •  T*  •  T9 

191  99.9  -  M4I9*|*|TI  I  |9f»  14*.  1*9 

1*9  9*1#  •  *|f*9l  9|f*.  9*|«IT  |  »  NTUH  Tf  /  99  I 

0*  TO  I9| 

t*9  9||*  •  *91*91  41401*.  9*1919  I 

M  T*  IT| 

1*9  Ml*  •  MTAM  Tf  I  T9  I 


|9IT||tf*ft09  09  f*TM 9  HKfUMM  •MKftfttf 
OKT9KT  09  94*01*91  »0*  TM  <MT9Vf  MMtltltt 
CHI  *49*41 
CHI  C94TCI 
I  104  IT  •  1. 


0.* 

M9I  01**1.  09T1  I 
0*1  4901  ACtCt*.  091994  I. 
r*l  *994  T4ICI9.  0*199*  I. 
I  4TICI9  /  OOtAtf  ) 

0*1  *901  *T9tK.  OMIT*  I. 

I  ATICI9  t  9MI4AT  I 

(VM  I9M  4T9MIC.  0*1  TT*  I. 
I  ATICI9  t  KHIf  I 

041  *901  A* H 1C.  091TT*  I. 

•MCI* 

04*414 


*T*9I |ff 
•40**4  • 
motor  • 

TI4I9L  • 

AT4HC  • 

At94|C  . 

AT40IC  • 

*TM|C  • 

*T*IIC  • 

•T*tlC  • 

•1944T  • 

4or»*r  • 

•**iit  •  i 
9  *0 

ll«H  •  0 

IIAlf  •  0 
1**11  »  * 
tom  •  -  t 

AfV9T9  •  9 

■H*H  •  I 
9Wt4  •  r> 

9tf  •  •  4 

■ITItt  •  0 
9991 TT  •  9 

IM  *001IT9f*TTi  I  194.  19*.  199 

MITfl*.  *1  MWf.1iMK.(*l«ri, 

I  *r90*t.4T*lTT.T909Tl.TMir|.r4l09T,T*l09».Kt*T». 

f  AtflTf.ACflf 4«*9l0vT.ACttft»9fI**9.T«99tS.T9f alt. 

t  T9flTf.TCflft 

MITfl*.  91  T*.*0 

IM  fOOMCITTlIt  I  10*.  199.  10* 

14*  CITTV*  I  H».  I0T.  IM 

*9‘  Tf  I*.  Tl 

OO  TO  199 

MITMA.  01 

CHI  MHIT 


191  191  MUTT  l  1 991. 1999. 1991 
1991  141  19ITC91I99I  I  |99l*|90t.|90f 
19*9  A*  •  0*1 
1901  cm  I  touf 


C  909941  94TM9  99*9  1M99WT  I  Of  191  Til 
I9|  MTH*9 


MVOVtlM  KH«C  l.l  *AM 


7 


lll»rc  OCK. 

Ururutiar  40*  t« 

C  W«i(W*HK  COTIac  cool  KMH  CCMftC 


fMit  u/*#rvTiar  u  oiiimd  »o  c«4Mj  thoic  ttn  >uni  mmicm  ca» 
•  f  Hfllio  IMMI  (Mllix  or  CQ*Mt 

cnwvn  ittffNMTi  #p*  urim  muuit 

Cmaoa  /C0|/  MKiMIMiTAHiH.liltr.rRlHiMllt 

»  »M'NN.&-*e*N 

cm*om  /coi/  *T»i>n.iM*cr.riMCf.tMOut.iMCfCH..i«MCTr.iMCr*, 
»M«#l.orMCf.**6T0l.0*eTCN.0f«X. MACON 
1C°^*I  /CfU/  oI!iMt>WM,W<ll,W<Mt*W,l>*T*|ll0i*|l*W'W>1*1* 
CO-«o^  /coie/  i(nM,ilin<t,iC(ilt,UCIi,K(IU,mi,Mli2. 
t  4*4 1#.  4*4104. 4a|*.ia{a2,4*|NX.iAiar.  mu, 

{  Tflt4;,rT4)4l.4CCl*l'f  l.iMXOf.ilXM.iCCftM. 

’  «tJ1  .ACC)*'.  4(0*#. (afoul. CaMIf.fNMtt.  (*•*}. 

*  (Wtl.fKMt.dKMt.INMr.llMU.IIAni.lMlItf, 

’  IWl 

£0MM  /to/**/  «  l*4I.Tl«4ff  .t|M4ll,ria*X.0r|INa,»Nil.MUr. 

I  M<*Oi.0Mt«i,\MiovT,t(WA«.}p<utr.t»wn.SMaei. 

*  4Cir<X,4C*<CI.4e>CC0N 

ccmo*  /coo/  Mtort  .na*(.*mt  #.#*4114  .min.mi  aoa.  Mia/,  crm*. 
C«MV«  /coi 00/  f ia*i t. *•/«**. r»{ at t, ret ii t.t note, ttami. Trio* r. 
!  TMttit)  .*«Miii»i(MKtr(ti,rn»i.r2i4i.M<*i. 

/  r  i  cal*.  »2cait.M<*.iT.»r»Mim.c»<wu  i.nai. 

«  r**Ktliti. rt'M*. 144*04. TaaaaM.taaaaL.raNXi. 

%  *«<•(.  i*na*# ,  r«acar 

crtMrx.  /Cn/o/i/  tHOiM.tlMr 

o/v«if  memo*  irATfacars  #o*  corn*  vaaaimci 
oo/ai*  mkiiioh  iwiii9*i«»,t3un,u<ni.Mi(»t 
cxvfi#  Miciitoa  tOMiaut 


tMDMIK  KMK  t.l  MM 

ttwukv  »  mvn  nmicn  mv  m  ouacia  wn«  coma*.  or  n*m 
«•»»  miiataf  immi  mi  or  muff 

I  III  I*  vaiao*  Mia#  mmh  aval* 

*f«"»  »HM  HMM 

l  ll*  I*  reran  tacrx  iwiti  f«cu 

•ffroi  NK«  iaUM 

*  !!•  It  M#t«t  mu  tMiH 

Miaoi  *  a  i  *#  v*iai  uw 

*•  ,M  W  MIMt  U1M  MHH  U|«l  A*  I  a# 

MOM  11*101  ACC  If  It  «l 
(Mill  1*0*11  l«K)  IM4I  IM*t 

„  ««•*»  l**M»  I  MIX  IMtt  (MUr 

II*  l»*  f *  UMtf  iiaiu  eiatx  ari«* 

MAM  MilX  M*ll  imwi  IN  M«A III) 

1M4I*  1MAH  IMOl  l«NI  IVUIl 

iiimi  Aotftr  uim  AsiMi  «mcm 

M  »*♦  f*  Mioa»  Mil#  ru*x  Mtaat  ni«#  i*h« 
inoar  f*4i#  (UIK  taiax  rata#  Ifua 
WTI  -  roa  MtTI  •  -|«  *  muff  AM 
to  ff  Aina  If  Meant. 

-  rot  M| T|  •  r|  muff  AAA  to 
AC  Aiwa  in  aMtlHtf  or  ru  * 

m  5*1  4i  *r*Aosui  c*otui  iiai  t carat »i 

MltX  ffarosni  tMMflfl  fMrOflll 

IM  I**  *«  *NOJ*0  AOt.lN  Mfroo 

AMW|  Mint  (0*41* 

CONftf  Mlttl  Ml  41*  Ml  AM  rvfXT  lift 

***  »*•  coatn  ffAtja  vAiNCf  Meat  a  nous  tiaci 

ANAMa  ANTflN  Mtfia  rtNAM  IK#*I 

cram  ArriMtlM rat i 
COMM  rwtreiaiCNii 

IM1  10*  If  tcida  Alicia  fMCta  laiwr  imai 

fafatf  MMM  MtUt  IMAM  INOtCt 

rviMI  mow  OfAAtf  Airttf  ACMTf 

Aticuv  ac i it  I  itfutf  tcartf  rrfitf 

turn  (uim 

ifua  lot  rircai  tM**r  m  exeat  cr  cmbaa 

Ml  Ml  *44*11  44  Mil  M|*4t  NHIII 

fix  CM*  04  rum  IM  I4CMI 

#••1  MfXI  rtrxi  »HCM  CVCM 

•  tfAtj  AMIN  ##MI  ******  pmin 

rvMt  #nur  Ail  ax  mmi  ittau 

»*«**•  AtM“.  H#NU  Atroai  44  MAI 

••Mr  llama  wofuf  aoaAi  Mm* 

n|«r*c  *4i»K  tUM*  rwifua  (MX 

«**«#*  lyfui  iiimi  iitnii 

«w*f  -  fltcrr  ror  MfMI*  TWM  muff 
“»  C44AaA«A  XI  atrala  tuoiwutiaf  PitCM 
*aoo  Iff  riftai  M#Mt  NM  44  Mil  USH  MM  fa* 

At  ►Ml  N0AI  1C*  IX 


lUXH/IIM  term  A.f 

(C*V*  f»*»e*rN»i  roa  tioari  vmiaxis 

erm+om  /coarfX/  vmar.rxait.cortx.oovMiif  i 

cnaaoa  / Cvaa/  it#.»,if,C4*M»,»«t,MitAt,ltAlAl,fl«||,»|tAI, 

•  l'W*tll|*#t.t4 

CIWVM  .*(*■( f/  *t((.#0(*MI(MMN(>#9u*|||| 

wxr  memo*  ftAT(*f*rs  #m  ci^au  vaatAXIf 

zz-  ”u:;:s  ran xisz?™-”™-™".'™.*,* 

onuut  MI4UIM  Mfo 

r«/:v*4»acr  tutiwau  *0*  umh  vaAiauei 

#  *oi*44  tacr  t’H.Hroi 

ft AfiafNtf  roa  SutMWtlM  vxiAXIf 
COaana  /mart/  **»#W|}f  j.ori *0A>*(IA 

("*y*1  U*®!4'  MCt.MOt.rt.thlMtj.iflfi.ijts.Kitj.jm,, 

coaaw  /cc*av/  Man 1 1 1 .(Mni •(# uir •  ti*r«i 

*****  rarcillpa  icai|*c«ts  fM  tlMOwtlK  AMlAXff 
ooni  ax(«oA.*rio 

OOWXf  Mtcif  10a  M 

MfCMrra  cum 

rouimracr  l#»tfw«ri  roa  tvuouti**  vaaiAXff 
IWIIHItCI  |*,f4| 

IIICIUWIM  Tt»«  aao  fiacatioa  imiMm  rot  ir.H  ruiam 
IINMIM  **111.  toil*  J«Mf  ll.iatli 

oauxr  MKltiiw  foivta 

4fUU44fXt  flATIWNIf  roa  4X4#  V 4* 1 4*4(1 
raulviiiacf  i i* • out i 

IXAl  *IT4  m«t( 

aria  aatvia.M.iaf.M.ta/a.aoo.  #*«.a.t*a.r*a.a.T*a^/ 

#414  «J»o/*M.«iaioi%xa.oiao,oiM.o*aa/ai.of.o*.xa.ofa.iHa* 
i  oiac.oxaa.otao/ 


t  acrvac  *.*  #im  * 


c...  !!!B  I?  °!  *•**  ••  Hiww  mmp  aorvx  cw#vax 

II*  IM  AMI  (M«.  #»|  |  |  it*. 

['*  t2.*l-i.*«l»*.wi,^.MN«*i.MMx.rxMr.rMMi.xiAC*a» 

n*i  rata  •  **i ax 

ratal  •  Mia  •  mi  a 

•««MI  •  MMOt 

•AW  «  rtlMM  •  #44 Maw 

IM  tri  AMI  CO* rx.  01  »  I  lit.  1X41. If* 

carter  •  racrx 

rraaev  •  f.«  •  faraav 
MMAV  •  KCVX 

XCrCa  •  orccx 

ura  rr#i»t  •  .  iiiiMfim 

lit  in  *aoi  term,  x  |  s  lit.  3fai.ua 

IM  *f*o  if. let  M#iaT.T*itf ,r>4ax.*atax.TAiar.r»i*t.XM£ 

IMI  MAt  •  -  »m*t  •  van#  ♦  fl.fr  #Wi»f  »  r  tnabx 

v*l*  •  -  #v*iac  r  r«(i»  *  Cl.fr  #m»»  »  •  minx 

Ilia  4*IWI  •  -  401(4*10*11 

Ml  ere  AMI  corex.  «ta  I  I  Iff.  Iffl.llf 

in  acA#»  if. in  Ariowt.ANAir.AMsx.ANiNX.AAiar.xxxrTAfx.Acaif. 

n  •.  **•**»•  WO**.  MAT  ?.  fomx  .INW  .xxrv 

MM  INI  •  -  rriort  •  n.r  *  imi  mint  »  •  anaox 

A«ACf  •  4*41  •  ANAI 

A»IN  •  ATime  •  A.»*  re  1,0  *  A#t«or  I  •  AN  I  NX 

•■INI  •  A*IN  •  an; a 

flfCAf  •  CCAfX 

NWW  »  t.l 

AIXUI  •  IMMt.1 

xcion  •  KCtati*i 

KCJM  •  KClIim 

ACCINt  •  AN*  III  .  ACCtftCII.  -  ACCltMU.  -  ACClfllll.  MIN) 
rNNNt*  •  eaarx 
IWX  •  M 


fa  IDM4II  aria. a.  art.#  » 

I*  rrraiTi  rria.o  t 

I«  loiwti  rrio.a  i 

ii  roTNAti  trrf.a  /  un  /  mi. a  i 

si::; '“"•*•  '■••• ' «»•  -•••  »•••  • 

::  n:-.:  ”■•••  •  *c  • 


•WMOWIIM  WW  i,»  mm  > 


lit  |f|  AM|  «WIM,  )  |  |M|  IMiiiu 

l»4  MM  If. Ml  tMAM.il 

1 _  _ WTtM,»MAl>,|M«AAf|tTtMtMA»MtAM«< 


•  Mill  •  AM  MW 

•MM  •  Mill  *«  I  MM  A 
«"!*  *  •**«*  9  ^INIII 

•  !WI1  /  I.MIMIII 

*«MM  •  M*V1  /  II 

•  MMt  -  NMIt  /  /.AMIMIAI 

....  •  *V*A»|  /  MMIHIfl 

tlA|  TIAAf  •  TI«Mi 

tl*»»  •  (,« 


NUI  •  •.• 

•MAI  •  MM 

•mai  •  nun 
•Siam  «  nun 
M  •  *.« 

m*  ••leut  •  -  aas i mown 
Wi%t  •  -  MtirriMtt 

lit  in  AMI  M»m.  AM  I  I  IM,  IM. 

1 H  MM  Ifc.MI  M,  tW  *  * 

Ml  OUT  •  Mill 
VflMt  f  twill 
Ml  Min  |  lit,  IM,  IM 

lit  MU«  •  M(||  «  | - 

t»Al#  •  full  A  I 

•  fin  •  a 

•  Mil  • 

IM  MfWH  •  AMI)  A  | 
mim  •  mm  a  i 

M|W  • 

MIA#  • 


MAI 

mu 

tM* 


mil 

•  MMOA 

t*mi 

tin 

•  AMIAA 

mn 

•  AAAWM 

win 

w  hum 

*wm 

♦  oi* 

Mill 

•  on 

mm 

*  Ml 

Mill 

•  *1* 

wm 

•  Oil 

Mill 

•  »l* 

TWI1I 

•  Oil 

t*|A| 

•  Olf 

Mill 

r  tit 

tWIAl 

*  *M 

tlli 

•  MAN* 

-  nt rut  «  A|| 

• 

♦ 

•now?  *  AM 

• 

1  1.0  • 

• 

i  i.o  a  mowr 

wtum  amt*  a.0 

iai  m  ami  io*mi.  pi m  i  i  i«z,  iin.iu 
i'.'Ji 

»KMI  •  W.A 

'MU  • 

tUWOl  •  I.NM 

t«M*A  •  TMMSIAI 
tAWMl  •  fwwMSIfl 

TAM»«  •  AMAttl  -  tUMilll,  -  fMMtltl,  -  rmilll, 

•  M*Ml^tMCMIM,  -  TMUlltl.  -  IUUIIII)  Ml<*« 


******  . 

M*  1*1  Ct  •  I.  * 
.Mwtnin  •  i.i 

arntniti  .  •.« 


miti 

•mt» 

♦Olttl 
Ml  ttltl 
mm  Mnn 
m»  tMAiiin 
iai  mwiim 


•I  Cult 
1  Of Cult 
•  t.e 
i  im.iAii.nit 


"•ttiiin  •  cwMiiri 

unuti 


,  «  *ma#mi»i  •  tiiti  iMM,m 

MAIM  .UMtilM,  •  IIMIIII.  -  TPIMtll>,  M  » 


IA«  Ml  AMI  U»m,  im  I  |  |44,  IAAI.IAA 

,4%  US 


Min 
i*ii  im  • 
ivm  • 


i.twiimm 

•M#M  •  ■( MM 
W  •  •*•*•! 
t.M  •  A|M 
t.«  •  « 


t*l  •  t.t  •  » 

...  it*  •**•!•  I  IM,  I  A*,  uv 
l*»  MtAt  •  MN  I  .CM  /  MTt|«  | 

***  Ml  AMI  AMtM.  MM  I  I  |AA.  |||.  |M 
!*♦  CAU  CMS  It 


HtMW  tO  MSItHIIAMM  NNIM  IMMWTIM  IAI  Ml 


WtfOVTi  M  i.l  f*M 
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MHVtIM  KNM  «.! 


UNt'UflM  M/Hf  cm*  «WMM  *CXI 


IMII  M«'M  II  MlfCHO  »a  CtlClOt'f  *X  t«ca«*«tl  10  M  UMt 
to  r«*  ixo  vtiufi  o*  t**  *uimn  at  «m  at 

iwimiu  f*«*  m»  mi 

(owk  Italian  »o«  weaai  waatttm 
co-m*.  tenant  yiiui.(nmiui 

(«w>i  /even/  i<'fciL(i«i/coi*ii/iot*iio«i/cef#n/iBr»2UM> 

toaaot  /<»*•/  v/0.O.*'».6x*Na.ttl0.8ti.l4*.tiNTM.Tlft*.0M**» 

t  M0VI0.l|»l.l» 

COMCM  /CM* */ 

I  *(0.f»X*l1.«f**lT 

tiw  /ca*M/  4*7f.tto.«co.t>«.4*o.*t0.if l0.4to.4O*.4T»t.otO 
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formance  of  the  attitude  control  system  were  found  to  be  moderate  and  pre¬ 
dictable.  The  maximum  temperature  increment  associated  with  them  was  found 
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of  the  footprint  due  to  rotation  was  found  to  be  within  70  miles  of  the,  value 
that  would  occur  for  the  corresponding  orbit  in  vacuum.  Oblateness  could 
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central  coordinate  system.  The  results  of  density  variation  in  the  Earth's 
atmosphere  were  not  serious.  Large  increases  in  the  maximum  skin  temperature 
occurred  only  when  extremely  large  spatially  random  density  disturbances  were 
encountered  by  the  vehicle.  . 


DD 


FOR* 

1  JAN  «4 


1473 


Security  Classification 


Security  Clgtiificrtion 


KEY  WORDS 

Reentry  -  heating 

Reentry  -  lifting  body 

Reentry  -  temperature  control 

Reentry  -  effect  of  earth  oblateness 

Reentry  -  effect  of  density  perturbations 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  name  end  address 
of  tho  contractor,  subcontractor,  grantee,  Department  of  De¬ 
fense  activity  or  other  organisation  ( corporate  miihot)  issuing 
thr  report, 

la.  REPORT  SECUMTY  CLASSIFICATION:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
“Restricted  Data"  is  included.  Marking  la  to  be  in  accord¬ 
ance  with  appropriate  security  regulations. 

2b.  GROUP:  Automatic  downgrading  la  specified  in  DoD  Di¬ 
rective  5200. 10  and  Armed  Forcer  Industrial  Manual.  Enter 
the  group  number.  Also,  when  apt  Ucable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ised. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  In  all  cases  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES  If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progresa,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  name(s)  of  autboifs)  as  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  initi'’U 
If  military,  show  rank  end  branch  of  service.  The  name  ot 
the  principal  author  is  an  absolute  minimum  requirement 

6.  REPORT  DATE:  Enter  the  date  of  the  report  es  dey, 
month,  year,  or  month,  year.  If  more  than  one  date  appears 
cn  the  report,  use  date  of  publication. 

7a.  TOTAL  NUMBER  OP  PAGES:  The  total  page  count 
ahouid  follow  normal  pagination  procedures,  be,,  enter  the 
number  of  pages  containing  information, 

7b.  K 'JMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report, 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  writteiv 

8b,  Sc,  A  Id.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

9b.  OTHER  F.EPORT  NUMBER(S):  If  the  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  numberfs). 

10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim¬ 
itations  on  further  dissemination  of  the  report,  other  than  those 

cro  m-55t 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  “Qualified  requesters  may  obtain  copies  of  thlc 
report  from  DDC.” 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized.  ” 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


"U.  S.  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC  Other  qualified  users 
shall  request  through 


(5)  "All  distribution  of  this  report  is  controlled.  Qual¬ 
ified  DDC  users  shall  request  through 


If  the  report  has  been  furnished  tc  the  Office  of  Technical 
Servicsa,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known.  J 

1L  SUPPLEMENTARY  NOTES:  Uee  for  additional  explana¬ 
tory  notes. 

1Z  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pry¬ 
ing  (or)  the  research  and  development-  Include  address. 

13-  ABSTRACT:  Enter  an  abstract  giving  s  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port.  If  additional  spsce  is  required,  a  continuation  sheet  shall 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  epd  with 
an  Indication  of  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  (TS).  (S).  (C),  or  (V). 

There  is  no  limitation  cn  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
indax  entries  for  cataloging  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  la  required.  Identi¬ 
fiers,  auch  as  equipment  model  designation,  trade  name,  military 
projact  coda  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rule?,  and  weights  is  optional 


Security  Classification 


